Publications

MicroRNAs (miRs) regulate a variety of cellular processes, and their impaired expression is involved in cancer. Silencing of tumor-suppressive miRs in cancer can occur through epigenetic modifications, including DNA methylation and histone deacetylation. We performed comparative miR profiling on cultured lung cancer cells before and after treatment with 5aza-deoxycytidine plus Trichostatin A to reverse DNA methylation and histone deacetylation, respectively. Several tens of miRs were strongly induced by such 'epigenetic therapy'. Two representatives, miR-512-5p (miR-512) and miR-373, were selected for further analysis. Both miRs were secreted in exosomes. Re-expression of both miRs augmented cisplatin-induced apoptosis and inhibited cell migration; miR-512 also reduced cell proliferation. TEAD4 mRNA was confirmed as a direct target of miR-512; likewise, miR-373 was found to target RelA and PIK3CA mRNA directly. Our results imply that miR-512 and miR-373 exert cell-autonomous and non-autonomous tumor-suppressive effects in lung cancer cells, where their re-expression may benefit epigenetic cancer therapy.

The number of Genome Wide Association Studies (GWAS) of schizophrenia is rapidly growing. However, the small effect of individual genes limits the number of reliably implicated genes, which are too few and too diverse to perform reliable pathway analysis; hence the biological roles of the genes implicated in schizophrenia are unclear. To overcome these limitations we combine GWAS with genome-wide expression data from human post-mortem brain samples of schizophrenia patients and controls, taking these steps: 1) Identify 36 GWAS-based genes which are expressed in our dataset. 2) Find a cluster of 19 genes with highly correlated expression. We show that this correlation pattern is robust and statistically significant. 3) GO-enrichment analysis of these 19 genes reveals significant enrichment of ion channels and calcium-related processes. This finding (based on analyzing a small number of coherently expressed genes) is validated and enhanced in two ways: First, the emergence of calcium channels and calcium signaling is corroborated by identifying proteins that interact with those encoded by the cluster of 19. Second, extend the 19 cluster genes into 1028 genes, whose expression is highly correlated with the cluster's average profile. When GO-enrichment analysis is performed on this extended set, many schizophrenia related pathways appear, with calcium-related processes enriched with high statistical significance. Our results give further, expression-based validation to GWAS results, support a central role of calcium-signaling in the pathogenesis of schizophrenia, and point to additional pathways potentially related to the disease.

Accurate prognosis and prediction of response to therapy are essential for personalized treatment of cancer. Even though many prognostic gene lists and predictors have been proposed, especially for breast cancer, highthroughput "omic" methods have so far not revolutionized clinical practice, and their clinical utility has not been satisfactorily established. Different prognostic gene lists have very few shared genes, the biological meaning of most signatures is unclear, and the published success rates are considered to be overoptimistic. This review examines critically the manner in which prognostic classifiers are derived using machine-learning methods and suggests reasons for the shortcomings and problems listed above. Two approaches that may hold hope for obtaining improved prognosis are presented. Both are based on using existing prior knowledge; one proposes combining molecular "omic" predictors with established clinical ones, and the second infers biologically relevant pathway deregulation scores for each tumor from expression data, and uses this representation to study and stratify individual tumors. Approaches such as the second one are referred to in the physics literature as "phenomenology"; they will, hopefully, play a significant role in future studies of cancer.

Systemic lupus erythematosus (SLE) is an autoimmune disease that can attack many different body organs; the triggering event is unknown. SLE has been associated with more than 100 different autoantibody reactivities - anti-dsDNA is prominent. Nevertheless, autoantibodies to dsDNA occur in only two-thirds of SLE patients. We previously reported the use of an antigen microarray to characterize SLE serology. We now report the results of an expanded study of serology in SLE patients and scleroderma (SSc) patients compared with healthy controls. The analysis validated and extended previous findings: two-thirds of SLE patients reacted to a large spectrum of self-molecules that overlapped with their reactivity to dsDNA; moreover, some SLE patients manifested a deficiency of natural IgM autoantibodies. Most significant was the finding that many SLE patients who were negative for autoantibodies to dsDNA manifested abnormal antibody responses to Epstein-Barr virus (EBV): these subjects made IgG antibodies to EBV antigens to which healthy subjects did not respond or they failed to make antibodies to EBV antigens to which healthy subjects did respond. This observation suggests that SLE may be associated with a defective immune response to EBV. The SSc patients shared many of these serological abnormalities with SLE patients, but differed from them in increased IgG autoantibodies to topoisomerase and centromere B; 84% of SLE patients and 58% of SSc patients could be detected by their abnormal antibodies to EBV. Hence an aberrant immune response to a ubiquitous viral infection such as EBV might set the stage for an autoimmune disease.

Pemphigus vulgaris (PV) is an autoimmune skin disease, which has been characterized by IgG autoantibodies to desmoglein 3. Here we studied the antibody signatures of PV patients compared with healthy subjects and with patients with two other autoimmune diseases with skin manifestations (systemic lupus erythematosus and scleroderma), using an antigen microarray and informatics analysis. We now report a previously unobserved phenomenon - patients with PV, compared with the healthy subjects and the two other diseases, show a significant decrease in IgG autoantibodies to a specific set of self-antigens. This novel finding demonstrates that an autoimmune disease may be associated with a loss of specific, healthy IgG autoantibodies and not only with a gain of specific, pathogenic IgG autoantibodies.

Signal-induced transcript isoform variation (TIV) includes alternative promoter usage as well as alternative splicing and alternative polyadenylation of mRNA. To assess the phenotypic relevance of signal-induced TIV, we employed exon arrays and breast epithelial cells, which migrate in response to the epidermal growth factor (EGF). We show that EGF rapidly - within one hour - induces widespread TIV in a significant fraction of the transcriptome. Importantly, TIV characterizes many genes that display no differential expression upon stimulus. In addition, similar EGF-dependent changes are shared by a panel of mammary cell lines. A functional screen, which utilized isoform-specific siRNA oligonucleotides, indicated that several isoforms play essential, non-redundant roles in EGF-induced mammary cell migration. Taken together, our findings highlight the importance of TIV in the rapid evolvement of a phenotypic response to extracellular signals. Copyright:

The transmembrane neural cell adhesion receptor L1 is a Wnt/b-catenin target gene expressed in many tumor types. In human colorectal cancer, L1 localizes preferentially to the invasive front of tumors and when overexpressed in colorectal cancer cells, it facilitates their metastasis to the liver. In this study, we investigated genes that are regulated in human colorectal cancer and by the L1-NF-kB pathway that has been implicated in liver metastasis. c-Kit was the most highly suppressed gene in both colorectal cancer tissue and the L1-NF-kB pathway. c-Kit suppression that resulted from L1-mediated signaling relied upon NF-kB, which directly inhibited the transcription of SP1, a major activator of the c-Kit gene promoter. Reconstituting c-Kit expression in L1- transfected cells blocked the biological effects conferred by L1 overexpression in driving motility and liver metastasis. We found that c-Kit expression in colorectal cancer cells is associated with a more pronounced epithelial morphology, along with increased expression of E-cadherin and decreased expression of Slug. Although c-Kit overexpression inhibited the motility and metastasis of L1-expressing colorectal cancer cells, it enhanced colorectal cancer cell proliferation and tumorigenesis, arguing that separate pathways mediate tumorigenicity and metastasis by c-Kit. Our findings provide insights into how colorectal cancer metastasizes to the liver, the most common site of dissemination in this cancer.

Motivation: Real time quantitative polymerase chain reaction (qPCR) is an important tool in quantitative studies of DNA and RNA molecules; especially in transcriptome studies, where different primer combinations allow identification of specific transcripts such as splice variants or precursor messenger RNA. Several softwares that implement various rules for optimal primer design are available. Nevertheless, as designing qPCR primers needs to be done manually, the repeated task is tedious, time consuming and prone to errors.Results: We used a set of rules to automatically design all possible exon-exon and intron-exon junctions in the human and mouse transcriptomes. The resulting database is included as a track in the UCSC genome browser, making it widely accessible and easy to use.

Deregulated proliferation is a hallmark of cancer cells. Here, we show that microRNA-10b* is a master regulator of breast cancer cell proliferation and is downregulated in tumoural samples versus matched peritumoural counterparts. Two canonical CpG islands (5kb) upstream from the precursor sequence are hypermethylated in the analysed breast cancer tissues. Ectopic delivery of synthetic microRNA-10b* in breast cancer cell lines or into xenograft mouse breast tumours inhibits cell proliferation and impairs tumour growth in vivo, respectively. We identified and validated in vitro and in vivo three novel target mRNAs of miR-10b* (BUB1, PLK1 and CCNA2), which play a remarkable role in cell cycle regulation and whose high expression in breast cancer patients is associated with reduced disease-free survival, relapse-free survival and metastasis-free survival when compared to patients with low expression. This also suggests that restoration of microRNA-10b* expression might have therapeutic promise.

Embryonic stem cells (ESCs) maintain high genomic plasticity, which is essential for their capacity to enter diverse differentiation pathways. Posttranscriptional modifications of chromatin histones play a pivotal role in maintaining this plasticity. We now report that one such modification, monoubiquitylation of histone H2B on lysine 120 (H2Bub1), catalyzed by the E3 ligase RNF20, increases during ESC differentiation and is required for efficient execution of this process. This increase is particularly important for the transcriptional induction of relatively long genes during ESC differentiation. Furthermore, we identify the deubiquitinase USP44 as a negative regulator of H2B ubiquitylation, whose downregulation during ESC differentiation contributes to the increase in H2Bub1. Our findings suggest that optimal ESC differentiation requires dynamic changes in H2B ubiquitylation patterns, which must occur in a timely and well-coordinated manner.

Aberrant activation of kinases has emerged to be a key event along with tumor progression, maintenance of tumor phenotype and response to anticancer treatments. This study documents the existence of an oncogenic autoregulatory feedback loop that includes the polo-like kinase-2 (Snk/Plk2) and mutant p53 proteins. Plk2 protein binds to and phosphorylates mutant p53, thereby potentiating its oncogenic activities. Phosphorylated mutant p53 binds more efficiently to p300, consequently strengthening its own transcriptional activity. Plk2 gene is regulated at a transcriptional level by both wt- and mutant p53 proteins. This leads to growth suppression or enhanced cell proliferation and chemoresistance, respectively. In turn, the siRNA-mediated knockdown of either mutant p53 or Plk2 proteins significantly curtails the growth properties of tumor cells and their chemo-resistance to anticancer treatments. Therefore, this paper identifies a novel tumor network including Plk2 and mutant p53 proteins whose triggering in response to DNA damage might disclose important implications for the treatment of human cancers.

Gene expression-based prediction of genomic copy number aberrations in the chromosomal region 12q13 to 12q15 that is flanked by MDM2 and CDK4 identified Wnt inhibitory factor 1 (WIF1) as a candidate tumor suppressor gene in glioblastoma. WIF1 encodes a secreted Wnt antagonist and was strongly downregulated in most glioblastomas as compared with normal brain, implying deregulation of Wnt signaling, which is associated with cancer. WIF1 silencing was mediated by deletion (7/69, 10%) or epigenetic silencing by promoter hypermethylation (29/110, 26%). Co-amplification of MDM2 and CDK4 that is present in 10% of glioblastomas was associated in most cases with deletion of the whole genomic region enclosed, including the WIF1 locus. This interesting pathogenetic constellation targets the RB and p53 tumor suppressor pathways in tandem, while simultaneously activating oncogenic Wnt signaling. Ectopic expression of WIF1 in glioblastoma cell lines revealed a dose-dependent decrease of Wnt pathway activity. Furthermore, WIF1 expression inhibited cell proliferation in vitro, reduced anchorage-independent growth in soft agar, and completely abolished tumorigenicity in vivo. Interestingly, WIF1 overexpression in glioblastoma cells induced a senescence-like phenotype that was dose dependent. These results provide evidence that WIF1 has tumor suppressing properties. Downregulation of WIF1 in 75% of glioblastomas indicates frequent involvement of aberrant Wnt signaling and, hence, may render glioblastomas sensitive to inhibitors of Wnt signaling, potentially by diverting the tumor cells into a senescence-like state.

Normal cells require continuous exposure to growth factors in order to cross a restriction point and commit to cell-cycle progression. This can be replaced by two short, appropriately spaced pulses of growth factors, where the first pulse primes a process, which is completed by the second pulse, and enables restriction point crossing. Through integration of comprehensive proteomic and transcriptomic analyses of each pulse, we identified three processes that regulate restriction point crossing: (1) The first pulse induces essential metabolic enzymes and activates p53-dependent restraining processes. (2) The second pulse eliminates, via the PI3K/AKT pathway, the suppressive action of p53, as well as (3) sets an ERK-EGR1 threshold mechanism, which digitizes graded external signals into an all-or-none decision obligatory for S phase entry. Together, our findings uncover two gating mechanisms, which ensure that cells ignore fortuitous growth factors and undergo proliferation only in response to consistent mitogenic signals.

The fact that there is very little if any overlap between the genes of different prognostic signatures for early-discovery breast cancer is well documented. The reasons for this apparent discrepancy have been explained by the limits of simple machine-learning identification and ranking techniques, and the biological relevance and meaning of the prognostic gene lists was questioned. Subsequently, proponents of the prognostic gene lists claimed that different lists do capture similar underlying biological processes and pathways. The present study places under scrutiny the validity of this claim, for two important gene lists that are at the focus of current large-scale validation efforts. We performed careful enrichment analysis, controlling the effects of multiple testing in a manner which takes into account the nested dependent structure of gene ontologies. In contradiction to several previous publications, we find that the only biological process or pathway for which statistically significant concordance can be claimed is cell proliferation, a process whose relevance and prognostic value was well known long before gene expression profiling. We found that the claims reported by others, of wider concordance between the biological processes captured by the two prognostic signatures studied, were found either to be lacking statistical rigor or were in fact based on addressing some other question.

A mutation within one allele of the p53 tumor suppressor gene can inactivate the remaining wild-type allele in a dominant-negative manner and in some cases can exert an additional oncogenic activity, known as mutant p53 gain of function (GOF). To study the role of p53 mutations in prostate cancer and to discriminate between the dominant-negative effect and the GOF activity of mutant p53, we measured, using microarrays, the expression profiles of three immortalized prostate epithelial cultures expressing wild-type, inactivated p53 or mutated p53. Analysis of these gene expression profiles showed that both inactivated p53 and p53 R175H mutant expression resulted in the upregulation of cell cycle progression genes. A second group, which was upregulated exclusively by mutant p53 R175H, was predominantly enriched in developmental genes. This group of genes included the Twist1, a regulator of metastasis and epithelial-mesenchymal transition (EMT). Twist1 levels were also elevated in metastatic prostate cancer-derived cell line DU145, in immortalized lung fibroblasts and in a subset of lung cancer samples, all in a mutant p53-dependent manner. p53 R175H mutant bearing immortalized epithelial cells showed typical features of EMT, such as higher expression of mesenchymal markers, lower expression of epithelial markers and enhanced invasive properties in vitro. The mechanism by which p53 R175H mutant induces Twist1 expression involves alleviation of the epigenetic repression. Our data suggest that Twist1 expression might be upregulated following p53 mutation in cancer cells.

The recent decade has witnessed a surge of physicists to biology. Some of the activities of the participating groups focus on bona fide physics questions, posed on biological systems (such as the physics of molecular motors, for example). Another kind of research in which physicists take part alongside computer scientists and applied mathematicians, deals with questions that are of direct interest to biologists; they come under the umbrella of computational and systems biology. The topic of these lectures lies at the most biological end of this spectrum, addressing problems of clinical relevance which were posed and initiated by biologists. The objective of these lectures is to help the curious physicist to learn and to understand more about this emerging, highly interdisciplinary field of research, by providing brief introductions to molecular biology and cancer research. This is followed by a cursory review of some recent research done by the "Domany group" and its collaborations with biological and clinical labs. Furthermore, we mention (mainly in footnotes) a small subset of studies in which physicists have contributed to this field during the past years. A more detailed review of recent contributions by physicists is beyond the scope of this introductory text. The introductory nature of these lecture notes naturally induces a strong bias regarding publications cited; consequently, these lecture notes do not provide a fair, historically correct and updated review of relevant literature.

Highly regenerative tissues such as blood must possess effective DNA damage responses (DDR) that balance long-term regeneration with protection from leukemogenesis. Hematopoietic stem cells (HSCs) sustain life-long blood production, yet their response to DNA damage remains largely unexplored. We report that human HSCs exhibit delayed DNA double-strand break rejoining, persistent γH2AX foci, and enhanced p53- and ASPP1-dependent apoptosis after γ-radiation compared to progenitors. p53 inactivation or Bcl-2 overexpression reduced radiation-induced apoptosis and preserved in vivo repopulating HSC function. Despite similar protection from irradiation-induced apoptosis, only Bcl-2-overexpressing HSCs showed higher self-renewal capacity, establishing that intact p53 positively regulates self-renewal independently from apoptosis. The reduced self-renewal of HSCs with inactivated p53 was associated with increased spontaneous γH2AX foci in secondary transplants of HSCs. Our data reveal distinct physiological roles of p53 that together ensure optimal HSC function: apoptosis regulation and prevention of γH2AX foci accumulation upon HSC self-renewal.

Patients with systemic lupus erythematosus (SLE) produce antibodies to many different self-antigens. Here, we investigated antibodies in SLE sera using an antigen microarray containing many hundreds of antigens, mostly self-antigens. The aim was to detect sets of antibody reactivities characteristic of SLE patients in each of various clinical states - SLE patients with acute lupus nephritis, SLE patients in renal remission, and SLE patients who had never had renal involvement. The analysis produced two novel findings: (i) an SLE antibody profile persists independently of disease activity and despite long-term clinical remission, and (ii) this SLE antibody profile includes increases in four specific immunoglobulin G (IgG) reactivities to double-stranded DNA (dsDNA), single-stranded DNA (ssDNA), Epstein-Barr virus (EBV) and hyaluronic acid; the profile also includes decreases in specific IgM reactivities to myeloperoxidase (MPO), CD99, collagen III, insulin-like growth factor binding protein 1 (IGFBP1) and cardiolipin. The reactivities together showed high sensitivity (> 93%) and high specificity for SLE (> 88%). A healthy control subject who had the SLE antibody profile was later found to develop clinical SLE. The present study did not detect antibody reactivities that differentiated among the various subgroups of SLE subjects with statistical significance. Thus, SLE is characterized by an enduring antibody profile irrespective of clinical state. The association of SLE with decreased IgM natural autoantibodies suggests that these autoantibodies might enhance resistance to SLE.

ID4 (inhibitor of DNA binding 4) is a member of a family of proteins that function as dominant-negative regulators of basic helix-loop-helix transcription factors. Growing evidence links ID proteins to cell proliferation, differentiation and tumorigenesis. Here we identify ID4 as a transcriptional target of gain-of-function p53 mutants R175H, R273H and R280K. Depletion of mutant p53 protein severely impairs ID4 expression in proliferating tumor cells. The protein complex mutant p53-E2F1 assembles on specific regions of the ID4 promoter and positively controls ID4 expression. The ID4 protein binds to and stabilizes mRNAs encoding pro-angiogenic factors IL8 and GRO-α. This results in the increase of the angiogenic potential of cancer cells expressing mutant p53. These findings highlight the transcriptional axis mutant p53, E2F1 and ID4 as a still undefined molecular mechanism contributing to tumor neo-angiogenesis.

Prostate cancer (PC) is a heterogeneous disease whose aggressive phenotype is the second leading cause of cancer-related death in men. The identification of key molecules and pathways that play a pivotal role in PC progression towards an aggressive form is crucial. A major effort towards this end has been taken by global analyses of gene expression profiles. However, the large body of data did not provide a definitive idea about the genes which are associated with the aggressive growth of PC. In order to identify such genes, we performed an interspecies comparison between several human data sets and high quality microarray data that we generated from the transgenic adenocarcinoma of mouse prostate (TRAMP) strain. TheTRAMP PC mimics the histological and pathological appearance as well as the aggressive phenotype of human PC (huPC). Analysis of the microarray data, derived from microdissected TRAMP specimens removed at different stages of the disease yielded genetic signatures delineating the TRAMP PC development and progression. Comparison of the TRAMP data with a set of genes representing the core expression signature of huPC yielded a limited set genes. Some of these genes are known predictors of poor prognosis in huPC. Interestingly, the modulation of genes responsible for the invasive phenotype of huPC occurs in TRAMP already during the transition to prostate intraepithelial neoplasia (PIN) and onwards to localized tumors. We therefore suggest that critical oncogenic events leading to an aggressive phenotype of huPC can be studied in the PIN stage of TRAMP.

Aging is often associated with a decline in hippocampus-dependent spatial memory. Here, we show that functional cell-mediated immunity is required for the maintenance of hippocampus-dependent spatial memory. Sudden imposition of immune compromise in young mice caused spatial memory impairment, whereas immune reconstitution reversed memory deficit in immune-deficient mice. Analysis of hippocampal gene expression suggested that immune-dependent spatial memory performance was associated with the expression of insulin-like growth factor (Igf1) and of genes encoding proteins related to presynaptic activity (Syt10, Cplx2). We further showed that memory loss in aged mice could be attributed to age-related attenuation of the immune response and could be reversed by immune system activation. Homeostatic-driven proliferation of lymphocytes, which expands the existing T cell repertoire, restored spatial memory deficits in aged mice. Thus, our results identify a novel function of the immune system in the maintenance of spatial memory and suggest an original approach for arresting or reversing age-associated memory loss.

We developed a method for estimating the positional distribution of transcription factor (TF) binding sites using ChIP-chip data, and applied it to recently published experiments on binding sites of nine TFs: OCT4, SOX2, NANOG, HNF1A, HNF4A, HNF6, FOXA2, USF1 and CREB1. The data were obtained from a genome-wide coverage of promoter regions from 8-kb upstream of the transcription start site (TSS) to 2-kb downstream. The number of target genes of each TF ranges from few hundred to several thousand. We found that for each of the nine TFs the estimated binding site distribution is closely approximated by a mixture of two components: a narrow peak, localized within 300-bp upstream of the TSS, and a distribution of almost uniform density within the tested region. Using Gene Ontology (GO) and Enrichment analysis, we were able to associate (for each of the TFs studied) the target genes of both types of binding with known biological processes. Most GO terms were enriched either among the proximal targets or among those with a uniform distribution of binding sites. For example, the three stemness-related TFs have several hundred target genes that belong to 'development' and 'morphogenesis' whose binding sites belong to the uniform distribution.

We investigate the distributions of the link overlap, P(Q), in three-dimensional Ising spin glasses. We use clustering methodology to identify a set of pairs of states from different Gibbs states and calculate its contribution to P(Q). We find that the distribution over this set does not become trivial as the system size increases.

L1-CAM, a neuronal cell adhesion receptor, is also expressed in a variety of cancer cells. Recent studies identified L1-CAM as a target gene of β-catenin-T-cell factor (TCF) signaling expressed at the invasive front of human colon cancer tissue. We found that L1-CAM expression in colon cancer cells lacking L1-CAM confers metastatic capacity, and mice injected in their spleen with such cells form liver metastases. We identified ADAM10, a metalloproteinase that cleaves the L1-CAM extracellular domain, as a novel target gene of β-catenin-TCF signaling. ADAM10 overexpression in colon cancer cells displaying endogenous L1-CAM enhanced L1-CAM cleavage and induced liver metastasis, and ADAM10 also enhanced metastasis in colon cancer cells stably transfected with L1-CAM. DNA microarray analysis of genes induced by L1-CAM in colon cancer cells identified a cluster of genes also elevated in a large set of human colon carcinoma tissue samples. Expression of these genes in normal colon epithelium was low. These results indicate that there is a gene program induced by L1-CAM in colon cancer cells that is also present in colorectal cancer tissue and suggest that L1-CAM can serve as target for colon cancer therapy.

Background. Transcription factors (TF) regulate expression by binding to specific DNA sequences. A binding event is functional when it affects gene expression. Functionality of a binding site is reflected in conservation of the binding sequence during evolution and in over represented binding in gene groups with coherent biological functions. Functionality is governed by several parameters such as the TF-DNA binding strength, distance of the binding site from the transcription start site (TSS), DNA packing, and more. Understanding how these parameters control functionality of different TFs in different biological contexts is a must for identifying functional TF binding sites and for understanding regulation of transcription. Methodology/Principal Findings. We introduce a novel method to screen the promoters of a set of genes with shared biological function (obtained from the functional Gene Ontology (GO) classification) against a precompiled library of motifs, and find those motifs which are statistically over-represented in the gene set. More than 8000 human (and 23,000 mouse) genes, were assigned to one of 134 GO sets. Their promoters were searched (from 200 bp downstream to 1000 bp upstream the TSS) for 414 known DNA motifs. We optimized the sequence similarity score threshold, independently for every location window, taking into account nucleotide heterogeneity along the promoters of the target genes. The method, combined with binding sequence and location conservation between human and mouse, identifies with high probability functional binding sites for groups of functionally-related genes. We found many location-sensitive functional binding events,and showed that they clustered close to the TSS. Our method and findings were tested experimentally. Conclusions/Significance. We identified reliably functional TF binding sites. This is an essential step towards constructing regulatory networks. The promoter region proximal to the TSS is of central importance for regulation of transcription in human and mouse, just as it is in bacteria and yeast.

Dietary antioxidants and selenium compounds were shown to have a therapeutic effect against hepatocellular carcinoma in several mouse models. We tested the effects of tannic acid and selenomethionine on hepatocellular carcinoma development in Mdr2 knockout (Mdr2-KO) mice. Mdr2-KO and age-matched Mdr2 heterozygous control mice were fed with tannic acid or selenomethionine during the first 3 months of life. Then, several mice from each group were sacrificed, and liver tissue samples were removed for analysis. The remaining mice were fed a regular diet until the age of 16 months, at which time the number and size of liver tumors were determined. Liver tissue samples of 3-month-old mice were subjected to gene expression profiling analysis using cDNA macroarrays containing probes for 240 genes that regulate responses to oxidative stress and inflammation or lipid metabolism. Both tannic acid and selenomethionine had partial chemopreventive effect on development of hepatocellular carcinoma in Mdr2-KO mice: they reduced the incidence of large tumor nodules (diameter > 1 cm) at age 16 months. Both agents inhibited gene expression and reversed up-regulation of many genes that control inflammation or response to oxidative stress in Mdr2-KO livers at age 3 months. This inhibitory effect on gene expression correlated with the ability of agents to reduce incidence of large tumors: selenomethionine was more active than tannic acid in both aspects. Understanding the molecular mechanism of chemoprevention effect could improve our therapeutic modalities while using these agents.

Signaling pathways invoke interplays between forward signaling and feedback to drive robust cellular response. In this study, we address the dynamics of growth factor signaling through profiling of protein phosphorylation and gene expression, demonstrating the presence of a kinetically defined cluster of delayed early genes that function to attenuate the early events of growth factor signaling. Using epidermal growth factor receptor signaling as the major model system and concentrating on regulation of transcription and mRNA stability, we demonstrate that a number of genes within the delayed early gene cluster function as feedback regulators of immediate early genes. Consistent with their role in negative regulation of cell signaling, genes within this cluster are downregulated in diverse tumor types, in correlation with clinical outcome. More generally, our study proposes a mechanistic description of the cellular response to growth factors by defining architectural motifs that underlie the function of signaling networks.

A bioinformatics-based analysis of endochondral bone formation model detected several genes upregulated in this process. Among these genes the dickkopf homolog 3 (Dkk3) was upregulated and further studies showed that its expression affects in vitro and in vivo osteogenesis. This study indicates a possible role of Dkk3 in regulating bone formation. Introduction: Endochondral bone formation is a complex biological process involving numerous chondrogenic, osteogenic, and angiogenic proteins, only some of which have been well studied. Additional key genes may have important roles as well. We hypothesized that to identify key genes and signaling pathways crucial for bone formation, a comprehensive gene discovery strategy should be applied to an established in vivo model of osteogenesis. Materials and Methods: We used in vivo implanted C3H10T1/2 cells that had been genetically engineered to express human bone morphogenetic protein-2 (BMP2) in a tetracycline-regulated system that controls osteogenic differentiation. Oligonucleotide microarray data from the implants (n = 4 repeats) was analyzed using coupled two-way clustering (CTWC) and statistical methods. For studying the effects of dickkopf homolog 3 (Dkk3) in chondrogenesis and osteogenesis. C3H10T1/2 mesenchymal progenitors were used. Results: The CTWC revealed temporal expression of Dkk3 with other chondrogenesis-, osteogenesis-, and Wnt-related genes. Quantitative RT-PCR confirmed the expression of Dkk3 in the implants. C3H10T1/2 cells that expressed Dkk3 in the presence of BMP2 displayed lower levels of alkaline phosphatase and collagen 1 mRNA expression than control C3H10T1/2 cells that did not express Dkk3. Interestingly, the levels of collagen II mRNA expression, Alcian blue staining, and glucose aminoglycans (GAGs) production were not influenced by Dkk3 expression. In vivo μCT and bioluminescence imaging revealed that co-expression of Dkk3 and BMP2 by implanted C3H10T1/2 cells induced the formation of significantly lower quantities of bone than cells expressing only BMP2. Conclusions: A bioinformatics analysis enabled the identification of Dkk3 as a pivotal gene with a novel function in endochondral bone formation. Our results showed that Dkk3 might have inhibitory effects on osteogenesis, but no effect on chondrogenesis, indicating that Dkk3 plays a regulatory role in endochondral bone formation. Further mechanistic studies are required to reveal the mechanism of action of Dkk3 in endochondral bone formation.

Background: The human genome contains over one million Alu repeat elements whose distribution is not uniform. While metabolism-related genes were shown to be enriched with Alu, in structural genes Alu elements are under-represented. Such observations led researchers to suggest that Alu elements were involved in gene regulation and were selected to be present in some genes and absent from others. This hypothesis is gaining strength due to findings that indicate involvement of Alu elements in a variety of functions; for example, Alu sequences were found to contain several functional transcription factor (TF) binding sites (BSs). We performed a search for new putative BSs on Alu elements, using a database of Position Specific Score Matrices (PSSMs). We searched consensus Alu sequences as well as specific Alu elements that appear on the 5 Kbp regions upstream to the transcription start site (TSS) of about 14000 genes. Results: We found that the upstream regions of the TSS are enriched with Alu elements, and the Alu consensus sequences contain dozens of putative BSs for TFs. Hence several TFs have Alu-associated BSs upstream of the TSS of many genes. For several TFs most of the putative BSs reside on Alu; a few of these were previously found and their association with Alu was also reported. In four cases the fact that the identified BSs resided on Alu went unnoticed, and we report this association for the first time. We found dozens of new putative BSs. Interestingly, many of the corresponding TFs are associated with early markers of development, even though the upstream regions of development-related genes are Alu-poor, compared with translational and protein biosynthesis related genes, which are Alu-rich. Finally, we found a correlation between the mouse B1 and human Alu densities within the corresponding upstream regions of orthologous genes. Conclusion: We propose that evolution used transposable elements to insert TF binding motifs into promoter regions. We observed enrichment of biosynthesis genes with Alu-associated BSs of developmental TFs. Since development and cell proliferation (of which biosynthesis is an essential component) were proposed to be opposing processes, these TFs possibly play inhibitory roles, suppressing proliferation during differentiation.

Predicting at the time of discovery the prognosis and metastatic potential of cancer is a major challenge in current clinical research. Numerous recent studies searched for gene expression signatures that outperform traditionally used clinical parameters in outcome prediction. Finding such a signature will free many patients of the suffering and toxicity associated with adjuvant chemotherapy given to them under current protocols, even though they do not need such treatment. A reliable set of predictive genes also will contribute to a better understanding of the biological mechanism of metastasis. Several groups have published lists of predictive genes and reported good predictive performance based on them. However, the gene lists obtained for the same clinical types of patients by different groups differed widely and had only very few genes in common. This lack of agreement raised doubts about the reliability and robustness of the reported predictive gene lists, and the main source of the problem was shown to be the small number of samples that were used to generate the gene lists. Here, we introduce a previously undescribed mathematical method, probably approximately correct (PAC) sorting, for evaluating the robustness of such lists. We calculate for several published data sets the number of samples that are needed to achieve any desired level of reproducibility. For example, to achieve a typical overlap of 50% between two predictive lists of genes, breast cancer studies would need the expression profiles of several thousand early discovery patients.

Keywords: Oncology; Radiology, Nuclear Medicine & Medical Imaging

Traditionally, immunologic diagnosis has been based on an attempt to correlate each disease with a specific immune reactivity, such as an antibody or a T-cell response to a single antigen specific for the disease entity. The state of the body, however, appears to be encoded by the immune system in collectives of reactivities and not by single reactivities. Here we describe our use of microarray technology and informatics to develop an antigen chip capable of detecting global patterns of antibodies binding to hundreds of antigens simultaneously. The patterns fashion diagnostic signatures.

Specific HPV DNA sequences are associated with more than 90% of invasive carcinomas of the uterine cervix. Viral E6 and E7 oncogenes are key mediators in cell transformation by disrupting TP53 and RB pathways. To investigate molecular mechanisms involved in the progression of invasive cervical carcinoma, we performed a gene expression study on cases selected according to viral and clinical parameters. Using Coupled Two-Way Clustering and Sorting Points Into Neighbourhoods methods, we identified a 'cervical cancer proliferation cluster' composed of 163 highly correlated transcripts. Most of these transcripts corresponded to E2F pathway genes controlling cell division or proliferation, whereas none was known as TP53 primary target. The average expression level of the genes of this cluster was higher in tumours with an early relapse than in tumours with a favourable course (P=0.026). Moreover, we found that E6/E7 mRNA expression level was positively correlated with the expression level of the cluster genes and with viral DNA load. These findings suggest that HPV E6/E7 expression level plays a key role in the progression of invasive carcinoma of the uterine cervix via the deregulation of cellular genes controlling tumour cell proliferation. HPV expression level may thus provide a biological marker useful for prognosis assessment and specific therapy of the disease.

On the basis of epidemiological studies, infection was suggested to play a role in the etiology of human cancer. While for some cancers such a role was indeed demonstrated, there is no direct biological support for the role of viral pathogens in the pathogenesis of childhood leukemia. Using a novel bioinformatic tool that alternates between clustering and standard statistical methods of analysis, we performed a 'double-blind' search of published gene expression data of subjects with different childhood acute lymphoblastic leukemia (ALL) subtypes, looking for unanticipated partitions of patients, induced by unexpected groups of genes with correlated expression. We discovered a group of about 30 genes, related to the interferon response pathway, whose expression levels divide the ALL samples into two subgroups; high in 50, low in 285 patients. Leukemic subclasses prevalent in early childhood (the age most susceptible to infection) are over-represented in the high-expression subgroup. Similar partitions, induced by the same genes, were found also in breast and ovarian cancer but not in lung cancer, prostate cancer and lymphoma. About 40% of breast cancer samples expressed the 'interferon-related' signature. It is of interest that several studies demonstrated mouse mammary tumor virus-like sequences in about 40% of breast cancer samples. Our discovery of an unanticipated strong signature of an interferon-induced pathway provides molecular support for a role for either inflammation or viral infection in the pathogenesis of childhood leukemia as well as breast and ovarian cancer.

The difficulty to dissect a complex phenotype of established malignant cells to several critical transcriptional programs greatly impends our understanding of the malignant transformation. The genetic elements required to transform some primary human cells to a tumorigenic state were described in several recent studies. We took the advantage of the global genomic profiling approach and tried to go one step further in the dissection of the transformation network. We sought to identify the genetic signatures and key target genes, which underlie the genetic alterations in p53, Ras, INK4A locus, and telotnerase, introduced in a stepwise manner into primary human fibroblasts. Here, we show that these are the minimally required genetic alterations for sarcomagenesis in vivo. A genome-wide expression profiling identified distinct genetic signatures corresponding to the genetic alterations listed above. Most importantly, unique transformation hallmarks, such as differentiation block, aberrant mitotic progression, increased angiogenesis, and invasiveness, were identified and coupled with genetic signatures assigned for the genetic alterations in the p53, INK4A locus, and H-Ras, respectively. Furthermore, a transcriptional program that defines the cellular response to p53 inactivation was an excellent predictor of metastasis development and bad prognosis in breast cancer patients. Deciphering these transformation fingerprints, which are affected by the most common oncogenic mutations, provides considerable insight into regulatory circuits controlling malignant transformation and will hopefully open new avenues for rational therapeutic decisions.

Motivation: Genes are often characterized dichotomously as either housekeeping or single-tissue specific. We conjectured that crucial functional information resides in genes with midrange profiles of expression. Results: To obtain such novel information genome-wide, we have determined the mRNA expression levels for one of the largest hitherto analyzed set of 62839 probesets in 12 representative normal human tissues. Indeed, when using a newly defined graded tissue specificity index τ, valued between 0 for housekeeping genes and 1 for tissue-specific genes, genes with midrange profiles having 0.15 50% of all expression patterns. We developed a binary classification, indicating for every gene the l_B tissues in which it is overly expressed, and the 12 - l_B tissues in which it shows low expression. The 85 dominant midrange patterns with l_B = 2-11 were found to be bimodally distributed, and to contribute most significantly to the definition of tissue specification dendrograms. Our analyses provide a novel route to infer expression profiles for presumed ancestral nodes in the tissue dendrogram. Such definition has uncovered an unsuspected correlation, whereby de novo enhancement and diminution of gene expression go hand in hand. These findings highlight the importance of gene suppression events, with implications to the course of tissue specification in ontogeny and phylogeny.

We calculate the Shannon entropy rate of a binary Hidden Markov Process (HMP), of given transition rate and noise ε (emission), as a series expansion in ε. The first two orders are calculated exactly. We then evaluate, for finite histories, simple upper-bounds of Cover and Thomas. Surprisingly, we find that for a fixed order k and history of n steps, the bounds become independent of n for large enough n. This observation is the basis of a conjecture, that the upper-bound obtained for n ≥ (k + 3)/2 gives the exact entropy rate for any desired order k of ε.

Motivation: Predicting the metastatic potential of primary malignant tissues has direct bearing on the choice of therapy. Several microarray studies yielded gene sets whose expression profiles successfully predicted survival. Nevertheless, the overlap between these gene sets is almost zero. Such small overlaps were observed also in other complex diseases, and the variables that could account for the differences had evoked a wide interest. One of the main open questions in this context is whether the disparity can be attributed only to trivial reasons such as different technologies, different patients and different types of analyses. Results: To answer this question, we concentrated on a single breast cancer dataset, and analyzed it by a single method, the one which was used by van't Veer et al. to produce a set of outcome-predictive genes. We showed that, in fact, the resulting set of genes is not unique; it is strongly influenced by the subset of patients used for gene selection. Many equally predictive lists could have been produced from the same analysis. Three main properties of the data explain this sensitivity: (1) many genes are correlated with survival; (2) the differences between these correlations are small; (3) the correlations fluctuate strongly when measured over different subsets of patients. A possible biological explanation for these properties is discussed.

Human embryonic stem cells (ESC) are undifferentiated and are endowed with the capacities of self-renewal and pluripotential differentiation. Adult stem cells renew their own tissue, but whether they can transdifferentiate to other tissues is still controversial. To understand the genetic program that underlies the pluripotency of stem cells, we compared the transcription profile of ESC with that of progenitor/stem cells of human hematopoietic and keratinocytic origins, along with their mature cells to be viewed as snapshots along tissue differentiation. ESC gene profiles show higher complexity with significantly more highly expressed genes than adult cells. We hypothesize that ESC use a strategy of expressing genes that represent various differentiation pathways and selection of only a few for continuous expression upon differentiation to a particular target. Such a strategy may be necessary for the pluripotency of ESC. The progenitors of either hematopoietic or keratinocytic cells also follow the same design principle. Using advanced clustering, we show that many of the ESC expressed genes are turned off in the progenitors/stem cells followed by a further down-regulation in adult tissues. Concomitantly, genes specific to the target tissue are up-regulated toward mature cells of skin or blood.

One's present repertoire of antibodies encodes the history of one's past immunological experience. Can the present autoantibody repertoire be consulted to predict resistance or susceptibility to the future development of an autoimmune disease? Here, we developed an antigen microarray chip and used bioinformatic analysis to study a model of type 1 diabetes developing in nonobese diabetic male mice in which the disease was accelerated and synchronized by exposing the mice to cyclophosphamide at 4 weeks of age. We obtained sera from 19 individual mice, treated the mice to induce cyclophosphamide-accelerated diabetes (CAD), and found, as expected, that 9 mice became severely diabetic, whereas 10 mice permanently resisted diabetes. We again obtained serum from each mouse after CAD induction. We then analyzed, by using rank-order and superparamagnetic clustering, the patterns of antibodies in individual mice to 266 different antigens spotted on the chip. A selected panel of 27 different antigens (10% of the array) revealed a pattern of IgG antibody reactivity in the pre-CAD sera that discriminated between the mice resistant or susceptible to CAD with 100% sensitivity and 82% specificity (P = 0.017). Surprisingly, the set of IgG antibodies that was informative before CAD induction did not separate the resistant and susceptible groups after the onset of CAD; new antigens became critical for post-CAD repertoire discrimination. Thus, at least for a model disease, present antibody repertoires can predict future disease, predictive and diagnostic repertoires can differ, and decisive information about immune system behavior can be mined by bioinformatic technology. Repertoires matter.

Informatic methodologies are being applied successfully to analyze the complexity of the genome. But beyond the genome, the immune system reflects the state of the body in health and disease. Traditionally, immunologists have reduced the immune system, where possible, to one-to-one relationships between particular antigens and particular antibodies or T-cell clones. Autoimmune diseases, caused by an immune attack against a body component, are usually investigated by following the response to single self-antigens. In this study, we apply informatics to analyze patterns of autoantibodies rather than single species of autoantibodies. This study was designed not to replace traditional approaches to immune diagnosis, but to test whether meaningful patterns of autoantibodies might exist. Using an unbiased solid-phase ELISA antibody test, we detected serum IgG and IgM antibodies in the sera of 20 healthy persons and 20 persons with type 1 diabetes mellitus binding to an array of 87 different antigens, mostly self-antigens. The healthy subjects manifested autoantibodies to a variety of self-antigens, many known to be associated with autoimmune diseases. We investigated the patterns of these autoantibodies using a coupled two-way clustering algorithm developed for analyzing data from gene arrays. We now report that the reactivity patterns of autoantibodies to particular subsets of self-antigens exhibited non-trivial structure, which significantly discriminated between healthy persons and persons with type 1 diabetes. The results show that despite the wide prevalence of autoantibodies, the patterns of reactivity to defined subsets of self-antigens can provide information about the state of the body.

Summary: The CTWC server provides access to the software, CTWC1.00, that implements Coupled Two Way Clustering (Getz et al., 2000), a method designed to mine gene expression data.

Using DNA microarray and clustering of expressed genes we have analyzed the mechanism of inhibition of wild-type p53-induced apoptosis by the cytokine interleukin 6 (IL-6) and the calcium mobilizer thapsigargin (TG). Clustering analysis of 1,786 genes, the expression level of which changed after activation of wild-type p53 in the absence or presence of IL-6 or TG, showed that these compounds did not cause a general inhibition of the ability of p53 to up-regulate or down-regulate gene expression. Expression of various p53 targets implicated as mediators of p53-induced apoptosis was also not affected by IL-6 or TG. These compounds thus can bypass the effect of wild-type p53 on gene expression and inhibit apoptosis. IL-6 and TG activated different p53-independent pathways of gene expression that include up-regulation of antiapoptotic genes. IL-6 and TG also activated different differentiationassociated genes. The ability of compounds such as cytokines and calcium mobilizers to inhibit p53-mediated apoptosis without generally inhibiting gene expression regulated by p53 can facilitate tumor development and tumor resistance to radiation and chemotherapy in cells that retain wild-type p53.

The expression levels of many thousands of genes can be measured simultaneously by DNA microarrays (chips). This novel experimental tool has revolutionized research in molecular biology and generated considerable excitement. A typical experiment uses a few tens of such chips, each dedicated to a single sample-such as tissue extracted from a particular tumor. The results of such an experiment contain several hundred thousand numbers, that come in the form of a table, of several thousand rows (one for each gene) and 50-100 columns (one for each sample). We developed a clustering methodology to mine such data. In this review I provide a very basic introduction to the subject, aimed at a physics audience with no prior knowledge of either gene expression or clustering methods. I explain what genes are, what is gene expression and how it is measured by DNA chips. Next I explain what is meant by "clustering" and how we analyze the massive amounts of data from such experiments, and present results obtained from analysis of data from colon cancer, brain tumors and breast cancer.

We present an automated procedure to assign CATH and SCOP classifications to proteins whose FSSP score is available. CATH classification is assigned down to the topology level, and SCOP classification is assigned to the fold level. Because the FSSP database is updated weekly, this method makes it possible to update also CATH and SCOP with the same frequency. Our predictions have a nearly perfect success rate when ambiguous cases are discarded. These ambiguous cases are intrinsic in any protein structure classification that relies on structural information alone. Hence, we introduce the "twilight zone for structure classification." We further suggest that to resolve these ambiguous cases, other criteria of classification, based also on information about sequence and function, must be used.

We introduce a method for validation of results obtained by clustering analysis of data. The method is based on resampling the available data. A figure of merit that measures the stability of clustering solutions against resampling is introduced. Clusters that are stable against resampling give rise to local maxima of this figure of merit. This is presented first for a one-dimensional data set, for which an analytic approximation for the figure of merit is derived and compared with numerical measurements. Next, the applicability of the method is demonstrated for higher-dimensional data, including gene microarray expression data.

The unbiased samples of ground states (GS) of the ±J short-range Ising spin glass possess a characteristics of the mean-field (MF) description, such as nontrivial and a hierarchical but nonultrametric, treelike structure. Clustering the GS showed that this structure was induced by correlated spin domains, which served as the cores of macroscopic zero energy excitations, taking the system from one state-cluster to another.

Two methods were proposed recently to derive energy parameters from known native protein conformations and corresponding sets of decoys. One is based on finding, by means of a perceptron learning scheme, energy parameters such that the native conformations have lower energies than the decoys. The second method maximizes the difference between the native energy and the average energy of the decoys, measured in terms of the width of the decoys' energy distribution (Z-score). Whereas the perceptron method is sensitive mainly to 'outlier' (i.e., extremal) decoys, the Z-score optimization is governed by the high density regions in decoy-space. We compare the two methods by deriving contact energies for two very different sets of decoys: the first obtained for model lattice proteins and the second by threading. We find that the potentials derived by the two methods are of similar quality and fairly closely related. This finding indicates that standard, naturally occurring sets of decoys are distributed in a way that yields robust energy parameters (that are quite insensitive to the particular method used to derive them). The main practical implication of this finding is that it is not necessary to fine-tune the potential search method to the particular set of decoys used. (C) 2000 Wiley-Liss, Inc.

We analyzed several energy functions for predicting the native state of proteins from an energy minimization procedure. We derived the parameters of a given energy function by imposing the basic requirement that the energy of the native conformation of a protein is lower than that of any conformation chosen from a set of decoys. Our work is motivated by a recent result which proved that the simple pairwise contact approximation of the energy is insufficient to satisfy simultaneously such a basic requirement for all the proteins in a database. Here, we investigate the reasons of such negative results and show how to improve the predictive power of methods based on energy minimization. We generated decoys by gapless threading, and we derive energy parameters by perceptron learning. We first considered hydrophobic contributions to the energy, defined in several ways, and showed that the additional hydrophobic terms enlarge slightly the number of proteins that can be stabilized together. Next, we performed various modifications of the pairwise energy term. We introduced (1) a distinction between inter-residue contacts on the surface and in the core of a protein and (2) a simple distance-dependent pairwise interaction in which a two-tier definition of contact replaces the original (single-tier) one. Our results suggest that a detailed treatment of the pairwise potential is likely to be more relevant than the consideration of other forces. (C) 2000 Wiley-Liss, Inc.

A simple model for flowing sand on an inclined plane is introduced. The model is related to recent experiments by Douady and Daerr and reproduces some of the experimentally observed features. Avalanches of intermediate size appear to be compact, placing the critical behavior of the model into the universality class of compact directed percolation. On very large scales, however, the avalanches break up into several branches, leading to a crossover from compact to ordinary directed percolation. Thus, systems of flowing granular matter on an inclined plane could serve as a first physical realization of directed percolation.

The aim of clustering is to partition data according to natural classes present in it. We proposed recently a method that makes no explicit assumption about the structure of the data and under very general and natural assumptions solves the clustering problem by evaluating thermal properties of a disordered (granular) magnet. The method was tested successfully on a variety of artificial and real-life problems; here we emphasize its application to analyze results obtained by a novel method of computer vision. The combination of these two techniques provides a powerful tool that succeeded to cluster properly 90 images of 6 objects on the basis of their pairwise dissimilarities. These dissimilarities, which constitute a highly non-metric set of pairwise distances between the images, form the input for clustering. A hierarchical organization of the images that agrees with human intuition, was obtained without assigning to the images coordinates in some abstract space.

We introduce and investigate a simple model to describe recent experiments by Douady and Daerr on flowing sand. The model reproduces experimentally observed compact avalanches, whose opening angle decreases linearly as a threshold is approached. On large scales the model exhibits a crossover from compact directed percolation to directed percolation; we predict similar behavior for the experimental system. We estimate the regime where \u201ctrue\u201d directed percolation morphology and exponents will be observed, providing the first experimental realization for this class of models. 1999

Damage spreading for Ising cluster dynamics is investigated numerically by using random numbers in a way that conforms with the notion of submitting the two evolving replicas to the same thermal noise. Two damage spreading transitions are found; damage does not spread either at low or high temperatures. We determine some critical exponents at the high-temperature transition point, which seem consistent with directed percolation.

We discuss the results of our attempt to predict the folded state of proteins by using the contact map representation. In order to use contact maps, one has to solve a few methodological problems. First, we propose an efficient way to explore the space of contact maps to generate creditable candidates for the native state. Second, we introduce a procedure to ensure that the generated maps are physical. Third, and most important, we present a method to derive contact energy parameters based on perceptron learning. The energy function must be able to discriminate between the native state and a set of decoys, simultaneously for all the proteins in a given database. We show that such an energy function exists if the candidates are produced by gapless threading, for a database of 153 proteins. If, however, we use as decoys maps of low energy that are generated by our dynamical scheme, not even one single protein can be stabilized. We conclude by discussing the perspectives of our approach.

Background: Two problems are of major importance in protein fold prediction: how to generate plausible conformations, and how to choose an energy function to identify the native state. Contact maps are a simple representation of protein structure and offer a promising framework to address these two issues. Results: In this work we develop Monte Carlo dynamics in contact map space. The procedure is divided into four steps: non- local dynamics, in which large-scale 'cluster' moves are performed (clusters are in approximate correspondence with secondary structure elements); local dynamics, in which secondary structure location is optimized; reconstruction, in which the physicality of the contact map is restored; and refinement, which consists of a further Monte Carlo energy minimization in real space. We demonstrate that such a dynamical procedure is effective in producing uncorrelated low-energy states. Conclusions: The procedure introduced in this paper very effectively generates a representative ensemble of conformations. We are able to show that existing sets of pairwise contact energy parameters are not suitable to single out the native state within this ensemble. The remaining outstanding issue in protein folding is to find an energy function that can discriminate the native state from decoys.

We simulated site dilute Ising models in d = 3 dimensions for several lattice sizes L. For each L singular thermodynamic quantities X were measured at criticality and their distributions P(X) were determined for ensembles of several thousand random samples. For L → ∞ the relative width of P(X) tends to a universal constant: there is no self-averaging. The width of the distribution of the sample (i) dependent pseudocritical temperatures T_c(i, L) scales as δT_c(L)∼L^-1/v and not as ∼L^-d/2. The sample dependence of X_i(T, L) enters dominantly, but not exclusively, via T_c(i, L).

We present a general definition of damage spreading in a pair of models. Using this general framework, one can define damage spreading in an objective manner that does not depend on the particular dynamic procedure that is being used. The formalism can be used for any spin-model or cellular automaton, with sequential or parallel update rules. At this point we present its application to the Domany-Kinzel cellular automaton in one dimension, this being the simplest model in which damage spreading has been found and studied extensively. We show that the active phase of this model consists of three subphases characterized by different damage-spreading properties.

We present two interesting results regarding damage spreading in ferromagnetic Ising models. First, we show that a damage spreading transition can occur in an Ising chain that evolves in contact with a thermal reservoir. Damage heals at low temperature and spreads at high T. The dynamic rules for the system's evolution for which such a transition is observed are as legitimate as the conventional rules (Glauber, Metropolis, heat bath). Our second result is that such transitions are not always in the directed percolation universality class.

We present a new approach for clustering, based on the physical properties of an inhomogeneous ferromagnetic model. We do not assume any structure of the underlying distribution of the data. A Potts spin is assigned to each data point and short range interactions between neighboring points are introduced. Spin-spin correlations, measured (by Monte Carlo procedure) in a superparamagnetic regime in which aligned domains appear, serve to partition the data points into clusters. Our method outperforms other algorithms for toy problems as well as for real data.

Keywords: Computer Science, Information Systems; Neurosciences

We investigate an antiferromagnetic triangular Ising model with anisotropic ferromagnetic interactions between further neighbors, originally proposed by Kitatani and Oguchi [J. Phys. Soc. Jpn. 57, 1344 (1988)]. The phase diagram as a function of temperature and the ratio between first- and second-neighbor interaction strengths is thoroughly examined. We search for a Kosterlitz-Thouless transition to a state with algebraic decay of correlations, calculating the correlation lengths on strips of width up to 15 sites by transfer-matrix methods. Phenomenological renormalization, conformal invariance arguments, the Roomany-Wyld approximation, and a direct analysis of the scaled mass gaps are used. Our results provide limited evidence that a Kosterlitz-Thouless phase is present. Alternative scenarios are discussed.

The critical behavior of a bond-disordered Ashkin-Teller model on a square lattice is investigated by intensive Monte Carlo simulations. A duality transformation is used to locate a critical plane of the disordered model. This critical plane corresponds to the line of critical points of the pure model, along which critical exponents vary continuously. Along this line the scaling exponent corresponding to randomness φ=(α/ν) varies continuously and is positive so that the randomness is relevant, and different critical behavior is expected for the disordered model. We use a cluster algorithm for the Monte Carlo simulations based on the Wolff embedding idea, and perform a finite size scaling study of several critical models, extrapolating between the critical bond-disordered Ising and bond-disordered four-state Potts models. The critical behavior of the disordered model is compared with the critical behavior of an anisotropic Ashkin-Teller model, which is used as a reference pure model. We find no essential change in the order parameters' critical exponents with respect to those of the pure model. The divergence of the specific heat C is changed dramatically. Our results favor a logarithmic type divergence at Tc, C∼lnL for the random-bond Ashkin-Teller and four-state Potts models and C∼ln lnL for the random-bond Ising model.

A simple model describing nucleation and growth following quench into a liquid-solid coexistence region is presented. The model is compared with small-angle X-ray scattering results on di-block copolymers.

The authors report calculations of the critical capacity of perceptrons that are subject to pattern-dependent noise. They also calculate the capacity of a perceptron whose weights take values on a shifted sphere, and show how this system resembles 'noisy' behaviour in some limits.

We present a novel cluster algorithm for Monte Carlo simulations of the fully frustrated Ising model on the square lattice. The new method does not suffer from problems of metastability, and is extremely efficient even at T=0. Our algorithm is a special case of a more general Monte Carlo simulation scheme. The general scheme unifies many cluster algorithms that were developed recently in order to accelerate Monte Carlo simulations.

A model-layered feed-forward neural network with an iterated learning stage is described. The system is capable of self-organizing its internal representations so that a large number of patterns can be perfectly stored and retrieved. The performance of the network is found to be robust against various degrading operations. Comparison of performance with an exactly solvable model is made.

The phase diagram of an Ising model on the hcp lattice with antiferromagnetic nearest-neighbor interactions J in the basal plane and coupling D (ferromagnetic or antiferromagnetic) between planes is analyzed. Previously obtained exact results predict the existence of a Lifshitz tricritical point in its phase diagram. In mean-field theory the model is shown to exhibit a twofold Lifshitz point and a phase diagram with paramagnetic, ferromagnetic, and modulated phases. Monte Carlo simulations that clearly indicate the tricritical nature of the Lifshitz point are reported, and series-expansion techniques are used to identify a commensurate phase in a low-temperature portion of the phase diagram.

The problem of recognizing noisy patterns is trivially solved by identifying the key pattern which has the largest overlap with the input. We show that a procedure that does this can be implemented in a simple layered neural network. Recognition is free of spurious states, the network never reaches a state of total confusion. Processing time scales as log M, where M is the number of stored patterns.

We have recently found a phase transition in the dynamics of a particle that diffuses in the presence of a hierarchical array of barriers. The transition occurs as R, a parameter that controls the relative size of the barriers, is varied. For R>Rc, the system is diffusive; as R Rc+ the diffusion constant vanishes. For R

A renormalization-group analysis of critical behavior near a Lifshitz tricritical point (LTP) is presented, with emphasis on the role played by new, momentum-dependent quartic terms. These result in new stable fixed points which determine the critical behavior. For some values of n (the number of order-parameter components) and m (the dimensionality of the softsubspace, characterized by quartic momentum-dependent inverse correlation functions), the renormalization-group recursion relations have two stable and accessible fixed points. However, one of these can never be reached in practice, due to a thermodynamic instability which results in a first-order phase transition. For m=d-1, one of the fixed points describes the critical dynamics of the usual n-vector spin model in (d-1) dimensions. This dynamic fixed point also characterizes LTP behavior for large n and n=1. In all other cases, the LTP has new exponents, which are not related to the dynamic model. Our results may be relevant to Lifshitz tricritical behavior in RbCaF3 and in some liquid-crystal systems.

Critical dynamics in d dimensions can be mapped onto a supersymmetric equilibrium problem in d+1 with a Lifshitz tricritical point. The existence of such a point is demonstrated for the Ising model on the hcp lattice by Monte Carlo simulation.

An x-y model with nearest-neighbor interaction V()=2J{1-[cos2(2)]p2} is studied by Monte Carlo simulation. The model has a transition from a paramagnetic to an algebraic (massless) phase. An increase of the parameter p, which changes the shape of the potential, brings about a change in the nature of the transition from continuous to first order.

We present a multiple scattering formalism that was developed in order to calculate elastic wave scattering by complex defects. The formalism is based on an expansion of the T matrix in terms of single scatterer T matrices; we include all single and double scattering processes. The method has been numerically implemented for the case of two spherical cavities; the interaction terms were found to be small, which indicates that other problems of interest, such as a subsurface spherical cavity or crack can also be treated.

Using the definition due to Landau, we present a closed-form expression for the optical-activity tensor for a cholesteric liquid crystal. The rotatory power is calculated in the pretransitional regime (above the isotropic-cholesteric transition), as well as in the ordered blue phase. The expression obtained for the first case reduces, in the low-chirality limit, to that of Dolganov et al. The result for the blue phase is new, and indicates the possibility of distinguishing between some of the ordered cubic structures that were recently proposed by various workers.

The Schrödinger equation is solved on a variety of fractal lattices using a recursive technique. In this method, the energy levels and wave functions on a lattice with Nn sites is calculated in terms of the corresponding quantities on a smaller lattice with Nn-1 sites, via a kind of decimation process. As n the resulting energy levels are discrete, very closely spaced, and highly degenerate. Smoothed densities of states have a wide variety of singularities.

For isotropic, centrosymmetric, ferroelectric systems described by an n>1 order parameter p, a new term in the Hamiltonian density, w[p2p-p(p2)], is allowed. It is relevant, in the renormalization-group sense, in d=4 dimensions, when n

The critical properties of Potts models with random bonds are considered in two dimensions. A position-space renormalization-group procedure, based on the Migdal-Kadanoff method, is developed. While all previous position-space calculations satisfied the Harris criterion and the resulting scaling relation only approximately, we found conditions under which these relations are exactly satisfied, and constructed our renormalization-group procedure accordingly. Numerical results for phase diagrams and thermodynamic functions for various random-bond Potts models are presented. In addition, some exact results obtained using a duality transformation, as well as an heuristic derivation of scaling properties that correspond to the percolation problem are given.

A position-space renormalization-group method, suitable for studying the localization properties of electrons in a disordered system, was developed. Two different approximations to a well-defined exact procedure were used. The first method is a perturbative treatment to lowest order in the intercell couplings. This yields a localization edge in three dimensions, with a fixed point at the band center (E=0) at a critical disorder σc7.0. In the neighborhood of the fixed point the localization length L is predicted to diverge as L∼(σ-σc+βE2)-ν. In two dimensions no fixed point is found, indicating localization even for small randomness, in agreement with Abrahams, Anderson, Licciardello, and Ramakrishnan. The second method is an application of the finite-lattice approximation, in which the intercell hopping between two (or more) cells is treated to infinite order in perturbation theory. To our knowledge, this method has not been previously used for quantum systems. Calculations based on this approximation were carried out in two dimensions only, yielding results that are in agreement with those of the lowest-order approximation.

The q=3-state Potts model undergoes a first-order transition in three dimensions, and a continuous one in two dimensions. For a system of finite thickness L, the transition will be either first order (for L>L^x) or continuous (for L

We study directed percolation in two dimensions. Critical properties, such as the threshold pc and the correlation length exponent , depend on the direction · at which percolation is observed. We observe crossover from (·=0) to (·0); similar crossover is found when the probabilities assigned to vertical and horizontal bonds, pV and pH, become unequal. For the case pH=1 we solved the problem exactly, obtaining pc (·)=[1+cot(14- ·)]-1,=2.