tigaR: integrative significance analysis of temporal differential gene expression induced by genomic abnormalities.

BACKGROUND: To determine which changes in the host cell genome are crucial for cervical carcinogenesis, a longitudinal in vitro model system of HPV-transformed keratinocytes was profiled in a genome-wide manner. Four cell lines affected with either HPV16 or HPV18 were assayed at 8 sequential time points for gene expression (mRNA) and gene copy number (DNA) using high-resolution microarrays. Available methods for temporal differential expression analysis are not designed for integrative genomic studies. RESULTS: Here, we present a method that allows for the identification of differential gene expression associated with DNA copy number changes over time. The temporal variation in gene expression is described by a generalized linear mixed model employing low-rank thin-plate splines. Model parameters are estimated with an empirical Bayes procedure, which exploits integrated nested Laplace approximation for fast computation. Iteratively, posteriors of hyperparameters and model parameters are estimated. The empirical Bayes procedure shrinks multiple dispersion-related parameters. Shrinkage leads to more stable estimates of the model parameters, better control of false positives and improvement of reproducibility. In addition, to make estimates of the DNA copy number more stable, model parameters are also estimated in a multivariate way using triplets of features, imposing a spatial prior for the copy number effect. CONCLUSION: With the proposed method for analysis of time-course multilevel molecular data, more profound insight may be gained through the identification of temporal differential expression induced by DNA copy number abnormalities. In particular, in the analysis of an integrative oncogenomics study with a time-course set-up our method finds genes previously reported to be involved in cervical carcinogenesis. Furthermore, the proposed method yields improvements in sensitivity, specificity and reproducibility compared to existing methods. Finally, the proposed method is able to handle count (RNAseq) data from time course experiments as is shown on a real data set.

The use of in vitro peptide binding profiles and in silico ligand-receptor interaction profiles to describe ligand-induced conformations of the retinoid X receptor alpha ligand-binding domain.

It is hypothesized that different ligand-induced conformational changes can explain the different interactions of nuclear receptors with regulatory proteins, resulting in specific biological activities. Understanding the mechanism of how ligands regulate cofactor interaction facilitates drug design. To investigate these ligand-induced conformational changes at the surface of proteins, we performed a time-resolved fluorescence resonance energy transfer assay with 52 different cofactor peptides measuring the ligand-induced cofactor recruitment to the retinoid X receptor-alpha (RXRalpha) in the presence of 11 compounds. Simultaneously we analyzed the binding modes of these compounds by molecular docking. An automated method converted the complex three-dimensional data of ligand-protein interactions into two-dimensional fingerprints, the so-called ligand-receptor interaction profiles. For a subset of compounds the conformational changes at the surface, as measured by peptide recruitment, correlate well with the calculated binding modes, suggesting that clustering of ligand-receptor interaction profiles is a very useful tool to discriminate compounds that may induce different conformations and possibly different effects in a cellular environment. In addition, we successfully combined ligand-receptor interaction profiles and peptide recruitment data to reveal structural elements that are possibly involved in the ligand-induced conformations. Interestingly, we could predict a possible binding mode of LG100754, a homodimer antagonist that showed no effect on peptide recruitment. Finally, the extensive analysis of the peptide recruitment profiles provided novel insight in the potential cellular effect of the compound; for the first time, we showed that in addition to the induction of coactivator peptide binding, all well-known RXRalpha agonists also induce binding of corepressor peptides to RXRalpha.

miRNA-520f Reverses Epithelial-to-Mesenchymal Transition by Targeting ADAM9 and TGFBR2.

Reversing epithelial-to-mesenchymal transition (EMT) in cancer cells has been widely considered as an approach to combat cancer progression and therapeutic resistance, but a limited number of broadly comprehensive investigations of miRNAs involved in this process have been conducted. In this study, we screened a library of 1120 miRNA for their ability to transcriptionally activate the E-cadherin gene CDH1 in a promoter reporter assay as a measure of EMT reversal. By this approach, we defined miR-520f as a novel EMT-reversing miRNA. miR-520f expression was sufficient to restore endogenous levels of E-cadherin in cancer cell lines exhibiting strong or intermediate mesenchymal phenotypes. In parallel, miR-520f inhibited invasive behavior in multiple cancer cell systems and reduced metastasis in an experimental mouse model of lung metastasis. Mechanistically, miR-520f inhibited tumor cell invasion by directly targeting ADAM9, the TGFß receptor TGFBR2 and the EMT inducers ZEB1, ZEB2, and the snail transcriptional repressor SNAI2, each crucial factors in mediating EMT. Collectively, our results show that miR-520f exerts anti-invasive and antimetastatic effects in vitro and in vivo, warranting further study in clinical settings. Cancer Res; 77(8); 2008-17. ©2017 AACR.

Crystal structure of the copper-containing quercetin 2,3-dioxygenase from Aspergillus japonicus.

Quercetin 2,3-dioxygenase is a copper-containing enzyme that catalyzes the insertion of molecular oxygen into polyphenolic flavonols. Dioxygenation catalyzed by iron-containing enzymes has been studied extensively, but dioxygenases employing other metal cofactors are poorly understood. We determined the crystal structure of quercetin 2,3-dioxygenase at 1.6 A resolution. The enzyme forms homodimers, which are stabilized by an N-linked heptasaccharide at the dimer interface. The mononuclear type 2 copper center displays two distinct geometries: a distorted tetrahedral coordination, formed by His66, His68, His112, and a water molecule, and a distorted trigonal bipyramidal environment, which additionally comprises Glu73. Manual docking of the substrate quercetin into the active site showed that the different geometries of the copper site might be of catalytic importance.

The nuclear receptor ligand-binding domain: a family-based structure analysis.

Nuclear receptors (NRs) are ligand-dependent transcription factors that play a central role in various physiological processes. The pharmaceutical industry has great interest in this gene-family for the discovery of novel or improved drugs for treatment of, for example, cancer, infertility, or diabetes. The usage of three-dimensional coordinates of protein structures to analyse and predict interactions with ligands is an important aspect of this process. All NR ligand-binding domains have a similar fold, which allows for comparison of the structures of their three main functional sites: the ligand-binding pocket, the cofactor-binding groove, and the dimerization interface. We performed an analysis of nearly one hundred NR ligand-binding domain structures, and identified the functionally important residues. The combined knowledge about the shape of the binding sites and the residues involved in the binding is important for drug design in two ways. First, knowledge about the location of residues that interact with a ligand in all crystal structures or in certain subfamilies assists in the design and docking of drugs. Second, similarities and differences in the residue types of the most frequent ligand- and cofactor-binding residues provide insight about potential cross-reactivity of ligands or cofactors.

Systemic miRNA-7 delivery inhibits tumor angiogenesis and growth in murine xenograft glioblastoma.

Tumor-angiogenesis is the multi-factorial process of sprouting of endothelial cells (EC) into micro-vessels to provide tumor cells with nutrients and oxygen. To explore miRNAs as therapeutic angiogenesis-inhibitors, we performed a functional screen to identify miRNAs that are able to decrease EC viability. We identified miRNA-7 (miR-7) as a potent negative regulator of angiogenesis. Introduction of miR-7 in EC resulted in strongly reduced cell viability, tube formation, sprouting and migration. Application of miR-7 in the chick chorioallantoic membrane assay led to a profound reduction of vascularization, similar to anti-angiogenic drug sunitinib. Local administration of miR-7 in an in vivo murine neuroblastoma tumor model significantly inhibited angiogenesis and tumor growth. Finally, systemic administration of miR-7 using a novel integrin-targeted biodegradable polymeric nanoparticles that targets both EC and tumor cells, strongly reduced angiogenesis and tumor proliferation in mice with human glioblastoma xenografts. Transcriptome analysis of miR-7 transfected EC in combination with in silico target prediction resulted in the identification of OGT as novel target gene of miR-7. Our study provides a comprehensive validation of miR-7 as novel anti-angiogenic therapeutic miRNA that can be systemically delivered to both EC and tumor cells and offers promise for miR-7 as novel anti-tumor therapeutic.

Correction: A Functional Screen Identifies Specific MicroRNAs Capable of Inhibiting Human Melanoma Cell Viability.

[This corrects the article DOI: 10.1371/journal.pone.0043569.].

A functional screen identifies specific microRNAs capable of inhibiting human melanoma cell viability.

Malignant melanoma is an aggressive form of skin cancer with poor prognosis. Despite improvements in awareness and prevention of this disease, its incidence is rapidly increasing. MicroRNAs (miRNAs) are a class of small RNA molecules that regulate cellular processes by repressing messenger RNAs (mRNAs) with partially complementary target sites. Several miRNAs have already been shown to attenuate cancer phenotypes, by limiting proliferation, invasiveness, tumor angiogenesis, and stemness. Here, we employed a genome-scale lentiviral human miRNA expression library to systematically survey which miRNAs are able to decrease A375 melanoma cell viability. We highlight the strongest inhibitors of melanoma cell proliferation, including the miR-15/16, miR-141/200a and miR-96/182 families of miRNAs and miR-203. Ectopic expression of these miRNAs resulted in long-term inhibition of melanoma cell expansion, both in vitro and in vivo. We show specifically miR-16, miR-497, miR-96 and miR-182 are efficient effectors when introduced as synthetic miRNAs in several melanoma cell lines. Our study provides a comprehensive interrogation of miRNAs that interfere with melanoma cell proliferation and viability, and offers a selection of miRNAs that are especially promising candidates for application in melanoma therapy.

Kisspeptin signaling is indispensable for neurokinin B, but not glutamate, stimulation of gonadotropin secretion in mice.

Kisspeptins (Kp), products of the Kiss1 gene that act via Gpr54 to potently stimulate GnRH secretion, operate as mediators of other regulatory signals of the gonadotropic axis. Mouse models of Gpr54 and/or Kiss1 inactivation have been used to address the contribution of Kp in the central control of gonadotropin secretion; yet, phenotypic and hormonal differences have been detected among the transgenic lines available. We report here a series of neuroendocrine analyses in male mice of a novel Gpr54 knockout (KO) model, generated by heterozygous crossing of a loxP-Gpr54/Protamine-Cre double mutant line. Gpr54-null males showed severe hypogonadotropic hypogonadism but retained robust responsiveness to GnRH. Gonadotropic responses to the agonist of ionotropic glutamate receptors, N-methyl-d-aspartate, were attenuated, but persisted, in Gpr54-null mice. In contrast, LH secretion after activation of metabotropic glutamate receptors was totally preserved in the absence of Gpr54 signaling. Detectable, albeit reduced, LH responses were also observed in Gpr54 KO mice after intracerebroventricular administration of galanin-like peptide or RF9, putative antagonist of neuropeptide FF receptors for the mammalian ortholog of gonadotropin-inhibiting hormone. In contrast, the stimulatory effect of senktide, agonist of neurokinin B (NKB; cotransmitter of Kiss1 neurons), was totally abrogated in Gpr54 KO males. Lack of Kp signaling also eliminated feedback LH responses to testosterone withdrawal. However, residual but sustained increases of FSH were detected in gonadectomized Gpr54 KO males, in which testosterone replacement failed to fully suppress circulating FSH levels. In sum, our study provides novel evidence for the relative importance of Kp-dependent vs. -independent actions of several key regulators of GnRH secretion, such as glutamate, galanin-like peptide, and testosterone. In addition, our data document for the first time the indispensable role of Kp signaling in mediating the stimulatory effects of NKB on LH secretion, thus supporting the hypothesis that NKB actions on GnRH neurons are indirectly mediated via its ability to regulate Kiss1 neuronal output.

EPR characterization of the mononuclear Cu-containing Aspergillus japonicus quercetin 2,3-dioxygenase reveals dramatic changes upon anaerobic binding of substrates.

Quercetin 2,3-dioxygenase (2,3QD) is a copper-containing dioxygenase that catalyses the oxidation of the flavonol quercetin to 2-protocatechuoylphloroglucinol carboxylic acid with concomitant production of carbon monoxide. In contrast to iron dioxygenases, very little is known about copper dioxygenases. We have characterized 2,3QD from the fungus Aspergillus japonicus by electron paramagnetic resonance spectroscopy (EPR). At pH 6.0, 2,3QD shows a mixture of two EPR species. The major form has parameters typical of type 2 Cu sites (g// = 2.330, A// = 13.7 mT), the minor one has a more distorted geometry (g// = 2.290, A// = 12.5 mT). Anaerobic addition of the substrate quercetin results in a different, single species EPR spectrum with g// = 2.336, A// = 11.4 mT, parameters, which are in-between those of the type 2 and type 1 Cu sites in the Peisach-Blumberg (g// vs. A//) plot. After turnover, a new EPR signal is observed, which is ascribed to the carboxylic acid ester product complex. This spectrum is similar to that of the native enzyme at pH 10.0 and has g-tensor parameters suggesting a trigonal bipyramidal site. Of a variety of flavonoids studied, only flavonols are able to bind to the copper centre of 2,3QD. Nine flavonols with different hydroxylation patterns at the A- and B-ring have been analysed. They cluster in two different regions of the Peisach-Blumberg plot and show that the presence of a 5-OH group has a large effect on the A// parameter. Several differences are noted between A. japonicus 2,3QD and the enzyme from A. niger German Collection of Microorganisms 821.