Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Drosophila Erect wing (Ewg) controls mitochondrial fusion during muscle growth and maintenance by regulation of the Opa1-like gene

Mitochondrial biogenesis and morphological changes are associated with tissue-specific functional demand, but the factors and pathways that regulate these processes have not been completely identified. A lack of mitochondrial fusion has been implicated in various developmental and pathological defects. The spatiotemporal regulation of mitochondrial fusion in a tissue such as muscle is not well understood. Here, we show in Drosophila indirect flight muscles (IFMs) that the nuclear-encoded mitochondrial inner membrane fusion gene, Opa1-like, is regulated in a spatiotemporal fashion by the transcription factor/co-activator Erect wing (Ewg). In IFMs null for Ewg, mitochondria undergo mitophagy and/or autophagy accompanied by reduced mitochondrial functioning and muscle degeneration. By following the dynamics of mitochondrial growth and shape in IFMs, we found that mitochondria grow extensively and fuse during late pupal development to form the large tubular mitochondria. Our evidence shows that Ewg expression during early IFM development is sufficient to upregulate Opa1-like, which itself is a requisite for both late pupal mitochondrial fusion and muscle maintenance. Concomitantly, by knocking down Opa1-like during early muscle development, we show that it is important for mitochondrial fusion, muscle differentiation and muscle organization. However, knocking down Opa1-like, after the expression window of Ewg did not cause mitochondrial or muscle defects. This study identifies a mechanism by which mitochondrial fusion is regulated spatiotemporally by Ewg through Opa1-like during IFM differentiation and growth.

DNA STRUCTURAL FEATURES AND ARCHITECTURE OF PROMOTER REGIONS PLAY A ROLE IN GENE RESPONSIVENESS OF S. CEREVISIAE

Gene expression is the most fundamental biological process, which is essential for phenotypic variation. It is regulated by various external (environment and evolution) and internal (genetic) factors. The level of gene expression depends on promoter architecture, along with other external factors. Presence of sequence motifs, such as transcription factor binding sites (TFBSs) and TATA-box, or DNA methylation in vertebrates has been implicated in the regulation of expression of some genes in eukaryotes, but a large number of genes lack these sequences. On the other hand, several experimental and computational studies have shown that promoter sequences possess some special structural properties, such as low stability, less bendability, low nucleosome occupancy, and more curvature, which are prevalent across all organisms. These structural features may play role in transcription initiation and regulation of gene expression. We have studied the relationship between the structural features of promoter DNA, promoter directionality and gene expression variability in S. cerevisiae. This relationship has been analyzed for seven different measures of gene expression variability, along with two different regulatory effect measures. We find that a few of the variability measures of gene expression are linked to DNA structural properties, nucleosome occupancy, TATA-box presence, and bidirectionality of promoter regions. Interestingly, gene responsiveness is most intimately correlated with DNA structural features and promoter architecture.

Efficient Dendrimer-DNA Complexation and Gene Delivery Vector Properties of Nitrogen-Core Poly(propyl ether imine) Dendrimer in Mammalian Cells

Dendrimers as vectors for gene delivery were established, primarily by utilizing few prominent dendrimer types so far. We report herein studies of DNA complexation efficacies and gene delivery vector properties of a nitrogen-core poly(propyl ether imine) (PETIM) dendrimer, constituted with 22 tertiary amine internal branches and 24 primary amines at the periphery. The interaction of the dendrimer with pEGFPDNA was evaluated through UV-vis, circular dichroism (CD) spectral studies, ethidium bromide fluorescence emission quenching, thermal melting, and gel retardation assays, from which most changes to DNA structure during complexation was found to occur at a weight ratio of dendrimer:DNA similar to 2:1. The zeta potential measurements further confirmed this stoichiometry at electroneutrality. The structure of a DNA oligomer upon dendrimer complexation was simulated through molecular modeling and the simulation showed that the dendrimer enfolded DNA oligomer along both major and minor grooves, without causing DNA deformation, in 1:1 and 2:1 dendrimer-to-DNA complexes. Atomic force microscopy (AFM) studies on dendrimer-pEGFP DNA complex showed an increase in the average z-height as a result of dendrimers decorating the DNA, without causing a distortion of the DNA structure. Cytotoxicity studies involving five different mammalian cell lines, using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide] (MTT) assay, reveal the dendrimer toxicity profile (IC50) values of similar to 400-1000 mu g mL(-1), depending on the cell line tested. Quantitative estimation, using luciferase assay, showed that the gene transfection was at least 100 times higher when compared to poly(ethylene imine) branched polymer, having similar number of cationic sites as the dendrimer. The present study establishes the physicochemical behavior of new nitrogen-core PETIM dendrimer-DNA complexes, their lower toxicities, and efficient gene delivery vector properties.

A 16-Gene Signature Distinguishes Anaplastic Astrocytoma from Glioblastoma

Anaplastic astrocytoma (AA; Grade III) and glioblastoma (GBM; Grade IV) are diffusely infiltrating tumors and are called malignant astrocytomas. The treatment regimen and prognosis are distinctly different between anaplastic astrocytoma and glioblastoma patients. Although histopathology based current grading system is well accepted and largely reproducible, intratumoral histologic variations often lead to difficulties in classification of malignant astrocytoma samples. In order to obtain a more robust molecular classifier, we analysed RT-qPCR expression data of 175 differentially regulated genes across astrocytoma using Prediction Analysis of Microarrays (PAM) and found the most discriminatory 16-gene expression signature for the classification of anaplastic astrocytoma and glioblastoma. The 16-gene signature obtained in the training set was validated in the test set with diagnostic accuracy of 89%. Additionally, validation of the 16-gene signature in multiple independent cohorts revealed that the signature predicted anaplastic astrocytoma and glioblastoma samples with accuracy rates of 99%, 88%, and 92% in TCGA, GSE1993 and GSE4422 datasets, respectively. The protein-protein interaction network and pathway analysis suggested that the 16-genes of the signature identified epithelial-mesenchymal transition (EMT) pathway as the most differentially regulated pathway in glioblastoma compared to anaplastic astrocytoma. In addition to identifying 16 gene classification signature, we also demonstrated that genes involved in epithelial-mesenchymal transition may play an important role in distinguishing glioblastoma from anaplastic astrocytoma.

Conditional silencing of topoisomerase I gene of Mycobacterium tuberculosis validates its essentiality for cell survival

Topoisomerases are an important class of enzymes for regulating the DNA transaction processes. Mycobacterium tuberculosis (Mtb) is one of the most formidable pathogens also posing serious challenges for therapeutic interventions. The organism contains only one type IA topoisomerase (Rv3646c), offering an opportunity to test its potential as a candidate drug target. To validate the essentiality of M.tuberculosis topoisomerase I (TopoI(Mt)) for bacterial growth and survival, we have generated a conditionally regulated strain of topoI in Mtb. The conditional knockdown mutant exhibited delayed growth on agar plate. In liquid culture, the growth was drastically impaired when TopoI expression was suppressed. Additionally, novobiocin and isoniazid showed enhanced inhibitory potential against the conditional mutant. Analysis of the nucleoid revealed its altered architecture upon TopoI depletion. These studies establish the essentiality of TopoI for the M.tuberculosis growth and open up new avenues for targeting the enzyme.

Predicting gene ontology annotations of orphan GWAS genes using protein-protein interactions

Background: The number of genome-wide association studies (GWAS) has increased rapidly in the past couple of years, resulting in the identification of genes associated with different diseases. The next step in translating these findings into biomedically useful information is to find out the mechanism of the action of these genes. However, GWAS studies often implicate genes whose functions are currently unknown; for example, MYEOV, ANKLE1, TMEM45B and ORAOV1 are found to be associated with breast cancer, but their molecular function is unknown. Results: We carried out Bayesian inference of Gene Ontology (GO) term annotations of genes by employing the directed acyclic graph structure of GO and the network of protein-protein interactions (PPIs). The approach is designed based on the fact that two proteins that interact biophysically would be in physical proximity of each other, would possess complementary molecular function, and play role in related biological processes. Predicted GO terms were ranked according to their relative association scores and the approach was evaluated quantitatively by plotting the precision versus recall values and F-scores (the harmonic mean of precision and recall) versus varying thresholds. Precisions of similar to 58% and similar to 40% for localization and functions respectively of proteins were determined at a threshold of similar to 30 (top 30 GO terms in the ranked list). Comparison with function prediction based on semantic similarity among nodes in an ontology and incorporation of those similarities in a k nearest neighbor classifier confirmed that our results compared favorably. Conclusions: This approach was applied to predict the cellular component and molecular function GO terms of all human proteins that have interacting partners possessing at least one known GO annotation. The list of predictions is available at http://severus.dbmi.pitt.edu/engo/GOPRED.html. We present the algorithm, evaluations and the results of the computational predictions, especially for genes identified in GWAS studies to be associated with diseases, which are of translational interest.

Role of DNA sequence based structural features of promoters in transcription initiation and gene expression

D Regulatory information for transcription initiation is present in a stretch of genomic DNA, called the promoter region that is located upstream of the transcription start site (TSS) of the gene. The promoter region interacts with different transcription factors and RNA polymerase to initiate transcription and contains short stretches of transcription factor binding sites (TFBSs), as well as structurally unique elements. Recent experimental and computational analyses of promoter sequences show that they often have non-B-DNA structural motifs, as well as some conserved structural properties, such as stability, bendability, nucleosome positioning preference and curvature, across a class of organisms. Here, we briefly describe these structural features, the differences observed in various organisms and their possible role in regulation of gene expression.

Multi-Temporal Satellite Image Analysis Using Gene Expression Programming

This paper discusses an approach for river mapping and flood evaluation to aid multi-temporal time series analysis of satellite images utilizing pixel spectral information for image classification and region-based segmentation to extract water covered region. Analysis of Moderate Resolution Imaging Spectroradiometer (MODIS) satellite images is applied in two stages: before flood and during flood. For these images the extraction of water region utilizes spectral information for image classification and spatial information for image segmentation. Multi-temporal MODIS images from ``normal'' (non-flood) and flood time-periods are processed in two steps. In the first step, image classifiers such as artificial neural networks and gene expression programming to separate the image pixels into water and non-water groups based on their spectral features. The classified image is then segmented using spatial features of the water pixels to remove the misclassified water region. From the results obtained, we evaluate the performance of the method and conclude that the use of image classification and region-based segmentation is an accurate and reliable for the extraction of water-covered region.

Ontology analysis of global gene expression differences of human bone marrow stromal cells cultured on 3D scaffolds or 2D films

Differences in gene expression of human bone marrow stromal cells (hBMSCs) during culture in three-dimensional (3D) nanofiber scaffolds or on two-dimensional (2D) films were investigated via pathway analysis of microarray mRNA expression profiles. Previous work has shown that hBMSC culture in nanofiber scaffolds can induce osteogenic differentiation in the absence of osteogenic supplements (OS). Analysis using ontology databases revealed that nanofibers and OS regulated similar pathways and that both were enriched for TGF-beta and cell-adhesion/ECM-receptor pathways. The most notable difference between the two was that nanofibers had stronger enrichment for cell-adhesion/ECM-receptor pathways. Comparison of nanofibers scaffolds with flat films yielded stronger differences in gene expression than comparison of nanofibers made from different polymers, suggesting that substrate structure had stronger effects on cell function than substrate polymer composition. These results demonstrate that physical (nanofibers) and biochemical (OS) signals regulate similar ontological pathways, suggesting that these cues use similar molecular mechanisms to control hBMSC differentiation. Published by Elsevier Ltd.

Methylation Silencing of ULK2, an Autophagy Gene, Is Essential for Astrocyte Transformation and Tumor Growth

Glioblastoma (GBM) is the most aggressive type of brain tumor and shows very poor prognosis. Here, using genome-wide methylation analysis, we show that G-CIMP+ and G-CIMP-subtypes enrich distinct classes of biological processes. One of the hypermethylated genes in GBM, ULK2, an upstream autophagy inducer, was found to be down-regulated in GBM. Promoter hypermethylation of ULK2 was confirmed by bisulfite sequencing. GBM and glioma cell lines had low levels of ULK2 transcripts, which could be reversed upon methylation inhibitor treatment. ULK2 promoter methylation and transcript levels showed significant negative correlation. Ectopic overexpression of ULK2-induced autophagy, which further enhanced upon nutrient starvation or temozolomide chemotherapy. ULK2 also inhibited the growth of glioma cells, which required autophagy induction as kinase mutant of ULK2 failed to induce autophagy and inhibit growth. Furthermore, ULK2 induced autophagy and inhibited growth in Ras-transformed immortalized Baby Mouse Kidney (iBMK) ATG5(+/+) but not in autophagy-deficient ATG5(-/-) cells. Growth inhibition due to ULK2 induced high levels of autophagy under starvation or chemotherapy utilized apoptotic cell death but not at low levels of autophagy. Growth inhibition by ULK2 also appears to involve catalase degradation and reactive oxygen species generation. ULK2 overexpression inhibited anchorage independent growth, inhibited astrocyte transformation in vitro and tumor growth in vivo. Of all autophagy genes, we found ULK2 and its homologue ULK1 were only down-regulated in all grades of glioma. Thus these results altogether suggest that inhibition of autophagy by ULK1/2 down-regulation is essential for glioma development.

A Novel Mutation in the Transglutaminase-1 Gene in an Autosomal Recessive Congenital Ichthyosis Patient

Structure-function implication on a novel homozygous Trp250/Gly mutation of transglutaminase-1 (TGM1) observed in a patient of autosomal recessive congenital ichthyosis is invoked from a bioinformatics analysis. Structural consequences of this mutation are hypothesized in comparison to homologous enzyme human factor XIIIA accepted as valid in similar structural analysis and are projected as guidelines for future studies at an experimental level on TGM1 thus mutated.

Co-liposomes comprising a lipidated multivalent RGD-peptide and a cationic gemini cholesterol induce selective gene transfection in alpha v beta 3 and alpha v beta 5 integrin receptor-rich cancer cells

The alpha v beta 3 and alpha v beta 5 integrins, transmembrane glycoprotein receptors, are over-expressed in numerous tumors and in endothelial cells that constitute tumor blood vessels. As this protein selectively binds to the Arg-Gly-Asp (RGD) sequence containing peptides, it is an attractive way to target tumors. Herein we have developed novel formulations for integrin mediated selective gene delivery. These formulations are composed of a novel palmitoylated tetrameric RGD containing scaffold (named RAFT-RGD), cationic gemini cholesterol (GL5) and a natural helper lipid 1,2-dioleoyl-L-alpha-glycero-3-phosphatidylethanolamine (DOPE). We have optimized a co-liposomal formulation to introduce the multivalent RGD-containing macromolecule in GL5: DOPE (GL5D) mixture to produce GL5D-RGD. We have unambiguously shown the selectivity of these formulations towards cancer cells that over express alpha v beta 3 and alpha v beta 5 integrins. Two reporter plasmids, pEGFP-C3 and PGL-3, were employed for the transfection experiments and it was shown that GL5D-RGD Liposomes increased exclusively the transfection in alpha v beta 3 and alpha v beta 5 overexpressing HeLa cells.