High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

Abstract Background The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery. Results Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer. Conclusions Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.

Quantum Mechanics Insight into the Microwave Nucleation of SrTiO3 Nanospheres

An extensive investigation of strontium titanate, SrTiO3 (STO), nanospheres synthesized via a microwave-assisted hydrothermal (MAH) method has been conducted to gain a better insight into thermodynamic, kinetic, and reaction phenomena involved in STO nucleation and crystal growth processes. To this end, quantum chemical modeling based on the density functional theory and periodic super cell models were done. Several experimental techniques were employed to get a deep characterization of structural and optical features of STO nanospheres. A possible formation mechanism was proposed, based on dehydration of titanium and strontium clusters followed by mesoscale transformation and a self-assembly process along an oriented attachment mechanism resulting in spherical like shape. Raman and XANES analysis renders a noncentrosymmetric environment for the octahedral titanium, while infrared and first order Raman modes reveal OH groups which are unsystematically incorporated into uncoordinated superficial sites. These results seem to indicate that the key component is the presence of distorted TiO6 clusters to engender a luminescence property. Analysis of band structure, density Of states, and charge map shows that there is a close relationship among local broken symmetry, polarization, and energy split of the 3d orbitals of titanium. The interplay among these electronic and structural features provides necessary conditions to evaluate its luminescent properties under two energy excitation.

Advanced QSAR Studies on PPAR delta Ligands Related to Metabolic Diseases

PPAR delta is a nuclear receptor that, when activated, regulates the metabolism of carbohydrates and lipids and is related to metabolic syndrome and type 2 diabetes. To understand the main interactions between ligands and PPAR delta, we have constructed 2D and 3D QSAR models and compared them with HOMO, LUMO and electrostatic potential maps of the compounds studied, as well as docking results. All QSAR models showed good statistical parameters and prediction outcomes. The QSAR models were used to predict the biological activity of an external test set, and the predicted values are in good agreement with the experimental results. Furthermore, we employed all maps to evaluate the possible interactions between the ligands and PPAR delta. These predictive QSAR models, along with the HOMO, LUMO and MEP maps, can provide insights into the structural and chemical properties that are needed in the design of new PPAR delta ligands that have improved biological activity and can be employed to treat metabolic diseases.

Dynamic Range of Vertebrate Retina Ganglion Cells: Importance of Active Dendrites and Coupling by Electrical Synapses

The vertebrate retina has a very high dynamic range. This is due to the concerted action of its diverse cell types. Ganglion cells, which are the output cells of the retina, have to preserve this high dynamic range to convey it to higher brain areas. Experimental evidence shows that the firing response of ganglion cells is strongly correlated with their total dendritic area and only weakly correlated with their dendritic branching complexity. On the other hand, theoretical studies with simple neuron models claim that active and large dendritic trees enhance the dynamic range of single neurons. Theoretical models also claim that electrical coupling between ganglion cells via gap junctions enhances their collective dynamic range. In this work we use morphologically reconstructed multi-compartmental ganglion cell models to perform two studies. In the first study we investigate the relationship between single ganglion cell dynamic range and number of dendritic branches/total dendritic area for both active and passive dendrites. Our results support the claim that large and active dendrites enhance the dynamic range of a single ganglion cell and show that total dendritic area has stronger correlation with dynamic range than with number of dendritic branches. In the second study we investigate the dynamic range of a square array of ganglion cells with passive or active dendritic trees coupled with each other via dendrodendritic gap junctions. Our results suggest that electrical coupling between active dendritic trees enhances the dynamic range of the ganglion cell array in comparison with both the uncoupled case and the coupled case with cells with passive dendrites. The results from our detailed computational modeling studies suggest that the key properties of the ganglion cells that endow them with a large dynamic range are large and active dendritic trees and electrical coupling via gap junctions.

Morphological Alterations in Newly Born Dentate Gyrus Granule Cells That Emerge after Status Epilepticus Contribute to Make Them Less Excitable

Computer simulations of external current stimulations of dentate gyrus granule cells of rats with Status Epilepticus induced by pilocarpine and control rats were used to evaluate whether morphological differences alone between these cells have an impact on their electrophysiological behavior. The cell models were constructed using morphological information from tridimensional reconstructions with Neurolucida software. To evaluate the effect of morphology differences alone, ion channel conductances, densities and distributions over the dendritic trees of dentate gyrus granule cells were the same for all models. External simulated currents were injected in randomly chosen dendrites belonging to one of three different areas of dentate gyrus granule cell molecular layer: inner molecular layer, medial molecular layer and outer molecular layer. Somatic membrane potentials were recorded to determine firing frequencies and inter-spike intervals. The results show that morphologically altered granule cells from pilocarpine-induced epileptic rats are less excitable than control cells, especially when they are stimulated in the inner molecular layer, which is the target area for mossy fibers that sprout after pilocarpine-induced cell degeneration. This suggests that morphological alterations may act as a protective mechanism to allow dentate gyrus granule cells to cope with the increase of stimulation caused by mossy fiber sprouting.

TGF-beta 1 modulates the homeostasis between MMPs and MMP inhibitors through p38 MAPK and ERK1/2 in highly invasive breast cancer cells

Background: Metastasis is the main factor responsible for death in breast cancer patients. Matrix metalloproteinases (MMPs) and their inhibitors, known as tissue inhibitors of MMPs (TIMPs), and the membrane-associated MMP inhibitor (RECK), are essential for the metastatic process. We have previously shown a positive correlation between MMPs and their inhibitors expression during breast cancer progression; however, the molecular mechanisms underlying this coordinate regulation remain unknown. In this report, we investigated whether TGF-beta 1 could be a common regulator for MMPs, TIMPs and RECK in human breast cancer cell models. Methods: The mRNA expression levels of TGF-beta isoforms and their receptors were analyzed by qRT-PCR in a panel of five human breast cancer cell lines displaying different degrees of invasiveness and metastatic potential. The highly invasive MDA-MB-231 cell line was treated with different concentrations of recombinant TGF-beta 1 and also with pharmacological inhibitors of p38 MAPK and ERK1/2. The migratory and invasive potential of these treated cells were examined in vitro by transwell assays. Results: In general, TGF-beta 2, T beta RI and T beta RII are over-expressed in more aggressive cells, except for T beta RI, which was also highly expressed in ZR-75-1 cells. In addition, TGF-beta 1-treated MDA-MB-231 cells presented significantly increased mRNA expression of MMP-2, MMP-9, MMP-14, TIMP-2 and RECK. TGF-beta 1 also increased TIMP-2, MMP-2 and MMP-9 protein levels but downregulated RECK expression. Furthermore, we analyzed the involvement of p38 MAPK and ERK1/2, representing two well established Smad-independent pathways, in the proposed mechanism. Inhibition of p38MAPK blocked TGF-beta 1-increased mRNA expression of all MMPs and MMP inhibitors analyzed, and prevented TGF-beta 1 upregulation of TIMP-2 and MMP-2 proteins. Moreover, ERK1/2 inhibition increased RECK and prevented the TGF-beta 1 induction of pro-MMP-9 and TIMP-2 proteins. TGF-beta 1-enhanced migration and invasion capacities were blocked by p38MAPK, ERK1/2 and MMP inhibitors. Conclusion: Altogether, our results support that TGF-beta 1 modulates the mRNA and protein levels of MMPs (MMP-2 and MMP-9) as much as their inhibitors (TIMP-2 and RECK). Therefore, this cytokine plays a crucial role in breast cancer progression by modulating key elements of ECM homeostasis control. Thus, although the complexity of this signaling network, TGF-beta 1 still remains a promising target for breast cancer treatment.

Regulation of neurogenesis and gliogenesis of retinoic acid-induced P19 embryonal carcinoma cells by P2X2 and P2X7 receptors studied by RNA interference

Embryonic carcinoma cells are widely used models for studying the mechanisms of proliferation and differentiation occurring during early embryogenesis. We have now investigated how down-regulation of P2X2 and P2X7 receptor expression by RNA interference (RNAi) affects neural differentiation and phenotype specification of P19 embryonal carcinoma cells. Wild-type P19 embryonal carcinoma cells or cells stably expressing shRNAs targeting P2X2 or P2X7 receptor expression were induced to differentiate into neurons and glial cells in the presence of retinoic acid. Silencing of P2X2 receptor expression along differentiation promoted cell proliferation and an increase in the percentage of cells expressing glial-specific GFAP, while the presence of beta-3 tubulin-positive cells diminished at the same time. Proliferation induction in the presence of stable anti-P2X2 receptor RNAi points at a mechanism where glial proliferation is favored over growth arrest of progenitor cells which would allow neuronal maturation. Differently from the P2X2 receptor, inhibition of P2X7 receptor expression during neural differentiation of P19 cells resulted in a decrease in cell proliferation and GFAP expression, suggesting the need of functional P2X7 receptors for the progress of gliogenesis. The results obtained in this study indicate the importance of purinergic signaling for cell fate determination during neural differentiation, with P2X2 and P2X7 receptors promoting neurogenesis and gliogenesis, respectively. The shRNAs down-regulating P2X2 or P2X7 receptor gene expression, developed during this work, present useful tools for studying mechanisms of neural differentiation in other stem cell models. (C) 2012 ISDN. Published by Elsevier Ltd. All rights reserved.

TGF-β1 modulates the homeostasis between MMPs and MMP inhibitors through p38 MAPK and ERK1/2 in highly invasive breast cancer cells

Abstract Background Metastasis is the main factor responsible for death in breast cancer patients. Matrix metalloproteinases (MMPs) and their inhibitors, known as tissue inhibitors of MMPs (TIMPs), and the membrane-associated MMP inhibitor (RECK), are essential for the metastatic process. We have previously shown a positive correlation between MMPs and their inhibitors expression during breast cancer progression; however, the molecular mechanisms underlying this coordinate regulation remain unknown. In this report, we investigated whether TGF-β1 could be a common regulator for MMPs, TIMPs and RECK in human breast cancer cell models. Methods The mRNA expression levels of TGF-β isoforms and their receptors were analyzed by qRT-PCR in a panel of five human breast cancer cell lines displaying different degrees of invasiveness and metastatic potential. The highly invasive MDA-MB-231 cell line was treated with different concentrations of recombinant TGF-β1 and also with pharmacological inhibitors of p38 MAPK and ERK1/2. The migratory and invasive potential of these treated cells were examined in vitro by transwell assays. Results In general, TGF-β2, TβRI and TβRII are over-expressed in more aggressive cells, except for TβRI, which was also highly expressed in ZR-75-1 cells. In addition, TGF-β1-treated MDA-MB-231 cells presented significantly increased mRNA expression of MMP-2, MMP-9, MMP-14, TIMP-2 and RECK. TGF-β1 also increased TIMP-2, MMP-2 and MMP-9 protein levels but downregulated RECK expression. Furthermore, we analyzed the involvement of p38 MAPK and ERK1/2, representing two well established Smad-independent pathways, in the proposed mechanism. Inhibition of p38MAPK blocked TGF-β1-increased mRNA expression of all MMPs and MMP inhibitors analyzed, and prevented TGF-β1 upregulation of TIMP-2 and MMP-2 proteins. Moreover, ERK1/2 inhibition increased RECK and prevented the TGF-β1 induction of pro-MMP-9 and TIMP-2 proteins. TGF-β1-enhanced migration and invasion capacities were blocked by p38MAPK, ERK1/2 and MMP inhibitors. Conclusion Altogether, our results support that TGF-β1 modulates the mRNA and protein levels of MMPs (MMP-2 and MMP-9) as much as their inhibitors (TIMP-2 and RECK). Therefore, this cytokine plays a crucial role in breast cancer progression by modulating key elements of ECM homeostasis control. Thus, although the complexity of this signaling network, TGF-β1 still remains a promising target for breast cancer treatment.

Structure-based drug design studies on a series of aldolase inhibitors

Human African trypanosomiasis, also known as sleeping sickness, is a major cause of death in Africa, and for which there are no safe and effective treatments available. The enzyme aldolase from Trypanosoma brucei is an attractive, validated target for drug development. A series of alkyl‑glycolamido and alkyl-monoglycolate derivatives was studied employing a combination of drug design approaches. Three-dimensional quantitative structure-activity relationships (3D QSAR) models were generated using the comparative molecular field analysis (CoMFA). Significant results were obtained for the best QSAR model (r2 = 0.95, non-cross-validated correlation coefficient, and q2 = 0.80, cross-validated correlation coefficient), indicating its predictive ability for untested compounds. The model was then used to predict values of the dependent variables (pKi) of an external test set,the predicted values were in good agreement with the experimental results. The integration of 3D QSAR, molecular docking and molecular dynamics simulations provided further insight into the structural basis for selective inhibition of the target enzyme.

Hausdorff Distance evaluation of orthodontic accessories' streaking artifacts in 3D model superimposition

The aim of this study was to determine whether image artifacts caused by orthodontic metal accessories interfere with the accuracy of 3D CBCT model superimposition. A human dry skull was subjected three times to a CBCT scan: at first without orthodontic brackets (T1), then with stainless steel brackets bonded without (T2) and with orthodontic arch wires (T3) inserted into the brackets' slots. The registration of image surfaces and the superimposition of 3D models were performed. Within-subject surface distances between T1-T2, T1-T3 and T2-T3 were computed and calculated for comparison among the three data sets. The minimum and maximum Hausdorff Distance units (HDu) computed between the corresponding data points of the T1 and T2 CBCT 3D surface images were 0.000000 and 0.049280 HDu, respectively, and the mean distance was 0.002497 HDu. The minimum and maximum Hausdorff Distances between T1 and T3 were 0.000000 and 0.047440 HDu, respectively, with a mean distance of 0.002585 HDu. In the comparison between T2 and T3, the minimum, maximum and mean Hausdorff Distances were 0.000000, 0.025616 and 0.000347 HDu, respectively. In the current study, the image artifacts caused by metal orthodontic accessories did not compromise the accuracy of the 3D model superimposition. Color-coded maps of overlaid structures complemented the computed Hausdorff Distances and demonstrated a precise fusion between the data sets.

Novel aryl beta-aminocarbonyl derivatives as inhibitors of Trypanosoma cruzi trypanothione reductase: binding mode revised by docking and GRIND2-based 3D-QSAR procedures

Trypanothione reductase has long been investigated as a promising target for chemotherapeutic intervention in Chagas disease, since it is an enzyme of a unique metabolic pathway that is exclusively present in the pathogen but not in the human host, which has the analog Glutathione reductase. In spite of the present data-set includes a small number of compounds, a combined use of flexible docking, pharmacophore perception, ligand binding site prediction, and Grid-Independent Descriptors GRIND2-based 3D-Quantitative Structure-Activity Relationships (QSAR) procedures allowed us to rationalize the different biological activities of a series of 11 aryl beta-aminocarbonyl derivatives, which are inhibitors of Trypanosoma cruzi trypanothione reductase (TcTR). Three QSAR models were built and validated using different alignments, which are based on docking with the TcTR crystal structure, pharmacophore, and molecular interaction fields. The high statistical significance of the models thus obtained assures the robustness of this second generation of GRIND descriptors here used, which were able to detect the most important residues of such enzyme for binding the aryl beta-aminocarbonyl derivatives, besides to rationalize distances among them. Finally, a revised binding mode has been proposed for our inhibitors and independently supported by the different methodologies here used, allowing further optimization of the lead compounds with such combined structure- and ligand-based approaches in the fight against the Chagas disease.

The 3D structure and function of digestive cathepsin L-like proteinases of Tenebrio molitor larval midgut

Cathepsin L-like proteinases (CAL) are major digestive proteinases in the beetle Tenebrio molitor. Procathepsin Ls 2 (pCAL2) and 3 (pCAL3) were expressed as recombinant proteins in Escherichia coil, purified and activated under acidic conditions. Immunoblot analyses of different T. molitor larval tissues demonstrated that a polyclonal antibody to pCAL3 recognized pCAL3 and cathepsin L 3 (CAD) only in the anterior two-thirds of midgut tissue and midgut luminal contents of T. molitor larvae. Furthermore, immunocytolocalization data indicated that pCAL3 occurs in secretory vesicles and microvilli in anterior midgut Therefore CAL3, like cathepsin L 2 (CAL2), is a digestive enzyme secreted by T. molitor anterior midgut CAD hydrolyses Z-FR-MCA and Z-RR-MCA (typical cathepsin substrates), whereas CAL2 hydrolyses only Z-FR-MCA. Active site mutants (pCAL2C25S and pCAL3C265) were constructed by replacing the catalytic cysteine with serine to prevent autocatalytic processing. Recombinant pCAL2 and pCAL3 mutants (pCAL2C25S and pCAL3C26S) were prepared, crystallized and their 3D structures determined at 1.85 and 2.1 angstrom, respectively. While the overall structure of these enzymes is similar to other members of the papain superfamily, structural differences in the S2 subsite explain their substrate specificities. The data also supported models for CAL trafficking to lysosomes and to secretory vesicles to be discharged into midgut contents. (C) 2012 Elsevier Ltd. All rights reserved.

Experimental Basis and New Insights for Cell Therapy in Chronic Obstructive Pulmonary Disease

Chronic Obstructive Pulmonary Disease (COPD) can be briefly described as air flow limitation and chronic dyspnea associated to an inflammatory response of the respiratory tract to noxious particles and gases. Its main feature is the obstruction of airflow and consequent chronic dyspnea. Despite recent advances, and the development of new therapeutic, medical and clinical approaches, a curative therapy is yet to be achieved. Therapies involving the use of tissue-specific or donor derived cells present a promising alternative in the treatment of degenerative diseases and injuries. Recent studies demonstrate that mesenchymal stem cells have the capacity to modulate immune responses in acute lung injury and pulmonary fibrosis in animal models, as well as in human patients. Due to these aspects, different groups raised the possibility that the stem cells from different sources, such as those found in bone marrow or adipose tissue, could act preventing the emphysematous lesion progression. In this paper, it is proposed a review of the current state of the art and future perspectives on the use of cell therapy in obstructive lung diseases.

Identifying regulational alterations in gene regulatory networks by state space representation of vector autoregressive models and variational annealing

Background: In the analysis of effects by cell treatment such as drug dosing, identifying changes on gene network structures between normal and treated cells is a key task. A possible way for identifying the changes is to compare structures of networks estimated from data on normal and treated cells separately. However, this approach usually fails to estimate accurate gene networks due to the limited length of time series data and measurement noise. Thus, approaches that identify changes on regulations by using time series data on both conditions in an efficient manner are demanded. Methods: We propose a new statistical approach that is based on the state space representation of the vector autoregressive model and estimates gene networks on two different conditions in order to identify changes on regulations between the conditions. In the mathematical model of our approach, hidden binary variables are newly introduced to indicate the presence of regulations on each condition. The use of the hidden binary variables enables an efficient data usage; data on both conditions are used for commonly existing regulations, while for condition specific regulations corresponding data are only applied. Also, the similarity of networks on two conditions is automatically considered from the design of the potential function for the hidden binary variables. For the estimation of the hidden binary variables, we derive a new variational annealing method that searches the configuration of the binary variables maximizing the marginal likelihood. Results: For the performance evaluation, we use time series data from two topologically similar synthetic networks, and confirm that our proposed approach estimates commonly existing regulations as well as changes on regulations with higher coverage and precision than other existing approaches in almost all the experimental settings. For a real data application, our proposed approach is applied to time series data from normal Human lung cells and Human lung cells treated by stimulating EGF-receptors and dosing an anticancer drug termed Gefitinib. In the treated lung cells, a cancer cell condition is simulated by the stimulation of EGF-receptors, but the effect would be counteracted due to the selective inhibition of EGF-receptors by Gefitinib. However, gene expression profiles are actually different between the conditions, and the genes related to the identified changes are considered as possible off-targets of Gefitinib. Conclusions: From the synthetically generated time series data, our proposed approach can identify changes on regulations more accurately than existing methods. By applying the proposed approach to the time series data on normal and treated Human lung cells, candidates of off-target genes of Gefitinib are found. According to the published clinical information, one of the genes can be related to a factor of interstitial pneumonia, which is known as a side effect of Gefitinib.

Proteasome inhibition and ROS generation by 4-nerolidylcatechol induces melanoma cell death

Induction of apoptotic cell death in response to chemotherapy and other external stimuli has proved extremely difficult in melanoma, leading to tumor progression, metastasis formation and resistance to therapy. A promising approach for cancer chemotherapy is the inhibition of proteasomal activity, as the half-life of the majority of cellular proteins is under proteasomal control and inhibitors have been shown to induce cell death programs in a wide variety of tumor cell types. 4-Nerolidylcatechol (4-NC) is a potent antioxidant whose cytotoxic potential has already been demonstrated in melanoma tumor cell lines. Furthermore, 4-NC was able to induce the accumulation of ubiquitinated proteins, including classic targets of this process such as Mcl-1. As shown for other proteasomal inhibitors in melanoma, the cytotoxic action of 4-NC is time-dependent upon the pro-apoptotic protein Noxa, which is able to bind and neutralize Mcl-1. We demonstrate the role of 4-NC as a potent inducer of ROS and p53. The use of an artificial skin model containing melanoma also provided evidence that 4-NC prevented melanoma proliferation in a 3D model that more closely resembles normal human skin.