Transcriptional Profiling Supports the Notochordal Origin of Chordoma and Its Dependence on a TGFB1-TBXT Network.

Chordoma is a rare malignant tumor demonstrating notochordal differentiation. It is dependent on brachyury (TBXT), a hallmark notochordal gene and transcription factor, and shares histologic features and the same anatomic location as the notochord. This study involved a molecular comparison of chordoma and notochord to identify dysregulated cellular pathways. The lack of a molecular reference from appropriate control tissue limits our understanding of chordoma and its relationship to notochord. Therefore, an unbiased comparison of chordoma, human notochord, and an atlas of normal and cancerous tissue was conducted using gene expression profiling to clarify the chordoma/notochord relationship and potentially identify novel drug targets. The study found striking consistency in gene expression profiles between chordoma and notochord, supporting the hypothesis that chordoma develops from notochordal remnants. A 12-gene diagnostic chordoma signature was identified and the TBXT/transforming growth factor beta (TGF-ß)/SOX6/SOX9 pathway was hyperactivated in the tumor, suggesting that pathways associated with chondrogenesis were a central driver of chordoma development. Experimental validation in chordoma cells confirmed these findings and emphasized the dependence of chordoma proliferation and survival on TGF-ß. The computational and experimental evidence provided the first molecular connection between notochord and chordoma and identified core members of a chordoma regulatory pathway involving TBXT. This pathway provides new therapeutic targets for this unique malignant neoplasm and highlights TGF-ß as a prime druggable candidate.

The IFITM proteins mediate cellular resistance to influenza A H1N1 virus, West Nile virus, and dengue virus.

Influenza viruses exploit host cell machinery to replicate, resulting in epidemics of respiratory illness. In turn, the host expresses antiviral restriction factors to defend against infection. To find host cell modifiers of influenza A H1N1 viral infection, we used a functional genomic screen and identified over 120 influenza A virus-dependency factors with roles in endosomal acidification, vesicular trafficking, mitochondrial metabolism, and RNA splicing. We discovered that the interferon-inducible transmembrane proteins IFITM1, 2, and 3 restrict an early step in influenza A viral replication. The IFITM proteins confer basal resistance to influenza A virus but are also inducible by interferons type I and II and are critical for interferon's virustatic actions. Further characterization revealed that the IFITM proteins inhibit the early replication of flaviviruses, including dengue virus and West Nile virus. Collectively this work identifies a family of antiviral restriction factors that mediate cellular innate immunity to at least three major human pathogens.

Genetic Analysis of Giant Cell Lesions of the Maxillofacial and Axial/Appendicular Skeletons.

PURPOSE: To compare the genetic and protein expression of giant cell lesions (GCLs) of the maxillofacial (MF) and axial/appendicular (AA) skeletons. We hypothesized that when grouped according to biologic behavior and not simply by location, MF and AA GCLs would exhibit common genetic characteristics. MATERIALS AND METHODS: This was a prospective and retrospective study of patients with GCLs treated at Massachusetts General Hospital from 1993 to 2008. In a preliminary prospective study, fresh tissue from 6 aggressive tumors each from the MF and AA skeletons (n = 12 tumors) was obtained. RNA was extracted and amplified from giant cells (GCs) and stromal cells first separated by laser capture microdissection. Genes highly expressed by GCs and stroma at both locations were determined using an Affymetrix GeneChip analysis. As confirmation, a tissue microarray (TMA) was created retrospectively from representative tissue of preserved pathologic specimens to assess the protein expression of the commonly expressed genes found in the prospective study. Quantification of immunohistochemical staining of MF and AA lesions was performed using Aperio image analysis to determine whether immunoreactivity was predictive of aggressive or nonaggressive behavior. RESULTS: Five highly ranked genes were found commonly in GCs and stroma at each location: matrix metalloproteinase-9 (MMP-9), cathepsin K (CTSK), T-cell immune regulator-1 (TCIRG1), C-type lectin domain family-11, and zinc finger protein-836. MF (n = 40; 32 aggressive) and AA (n = 48; 28 aggressive) paraffin-embedded tumors were included in the TMA. The proteins CTSK, MMP-9, and TCIRG1 were confirmed to have abundant expression within both MF and AA lesions. Only the staining levels for TCIRG1 within the GCs predicted the clinical behavior of the MF lesions. CONCLUSIONS: MMP-9, CTSK, and TCIRG1 are commonly expressed by GCLs of the MF and AA skeletons. This supports the hypothesis that these lesions are similar but at different locations. TCIRG1 has not been previously associated with GCLs and could be a potential target for molecular diagnosis and/or therapy.

Genome-wide analysis of immune system genes by expressed sequence Tag profiling.

Profiling studies of mRNA and microRNA, particularly microarray-based studies, have been extensively used to create compendia of genes that are preferentially expressed in the immune system. In some instances, functional studies have been subsequently pursued. Recent efforts such as the Encyclopedia of DNA Elements have demonstrated the benefit of coupling RNA sequencing analysis with information from expressed sequence tags (ESTs) for transcriptomic analysis. However, the full characterization and identification of transcripts that function as modulators of human immune responses remains incomplete. In this study, we demonstrate that an integrated analysis of human ESTs provides a robust platform to identify the immune transcriptome. Beyond recovering a reference set of immune-enriched genes and providing large-scale cross-validation of previous microarray studies, we discovered hundreds of novel genes preferentially expressed in the immune system, including noncoding RNAs. As a result, we have established the Immunogene database, representing an integrated EST road map of gene expression in human immune cells, which can be used to further investigate the function of coding and noncoding genes in the immune system. Using this approach, we have uncovered a unique metabolic gene signature of human macrophages and identified PRDM15 as a novel overexpressed gene in human lymphomas. Thus, we demonstrate the utility of EST profiling as a basis for further deconstruction of physiologic and pathologic immune processes.

Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology.

Genome-wide association studies (GWAS) have defined over 150 genomic regions unequivocally containing variation predisposing to immune-mediated disease. Inferring disease biology from these observations, however, hinges on our ability to discover the molecular processes being perturbed by these risk variants. It has previously been observed that different genes harboring causal mutations for the same Mendelian disease often physically interact. We sought to evaluate the degree to which this is true of genes within strongly associated loci in complex disease. Using sets of loci defined in rheumatoid arthritis (RA) and Crohn's disease (CD) GWAS, we build protein-protein interaction (PPI) networks for genes within associated loci and find abundant physical interactions between protein products of associated genes. We apply multiple permutation approaches to show that these networks are more densely connected than chance expectation. To confirm biological relevance, we show that the components of the networks tend to be expressed in similar tissues relevant to the phenotypes in question, suggesting the network indicates common underlying processes perturbed by risk loci. Furthermore, we show that the RA and CD networks have predictive power by demonstrating that proteins in these networks, not encoded in the confirmed list of disease associated loci, are significantly enriched for association to the phenotypes in question in extended GWAS analysis. Finally, we test our method in 3 non-immune traits to assess its applicability to complex traits in general. We find that genes in loci associated to height and lipid levels assemble into significantly connected networks but did not detect excess connectivity among Type 2 Diabetes (T2D) loci beyond chance. Taken together, our results constitute evidence that, for many of the complex diseases studied here, common genetic associations implicate regions encoding proteins that physically interact in a preferential manner, in line with observations in Mendelian disease.

Peripheral blood cellular dynamics of rheumatoid arthritis treatment informs about efficacy of response to disease modifying drugs.

Rheumatoid arthritis (RA) is an autoimmune disease characterized by systemic inflammation and is mediated by multiple immune cell types. In this work, we aimed to determine the relevance of changes in cell proportions in peripheral blood mononuclear cells (PBMCs) during the development of disease and following treatment. Samples from healthy blood donors, newly diagnosed RA patients, and established RA patients that had an inadequate response to MTX and were about to start tumor necrosis factor inhibitors (TNFi) treatment were collected before and after 3 months of treatment. We used in parallel a computational deconvolution approach based on RNA expression and flow cytometry to determine the relative cell-type frequencies. Cell-type frequencies from deconvolution of gene expression indicate that monocytes (both classical and non-classical) and CD4+ cells (Th1 and Th2) were increased in RA patients compared to controls, while NK cells and B cells (naïve and mature) were significantly decreased in RA patients. Treatment with MTX caused a decrease in B cells (memory and plasma cell), and a decrease in CD4 Th cells (Th1 and Th17), while treatment with TNFi resulted in a significant increase in the population of B cells. Characterization of the RNA expression patterns found that most of the differentially expressed genes in RA subjects after treatment can be explained by changes in cell frequencies (98% and 74% respectively for MTX and TNFi).

Gene enrichment profiles reveal T-cell development, differentiation, and lineage-specific transcription factors including ZBTB25 as a novel NF-AT repressor.

The identification of transcriptional regulatory networks, which control tissue-specific development and function, is of central importance to the understanding of lymphocyte biology. To decipher transcriptional networks in T-cell development and differentiation we developed a browsable expression atlas and applied a novel quantitative method to define gene sets most specific to each of the represented cell subsets and tissues. Using this system, body atlas size datasets can be used to examine gene enrichment profiles from a cell/tissue perspective rather than gene perspective, thereby identifying highly enriched genes within a cell type, which are often key to cellular differentiation and function. A systems analysis of transcriptional regulators within T cells during different phases of development and differentiation resulted in the identification of known key regulators and uncharacterized coexpressed regulators. ZBTB25, a BTB-POZ family transcription factor, was identified as a highly T cell-enriched transcription factor. We provide evidence that ZBTB25 functions as a negative regulator of nuclear factor of activated T cells (NF-AT) activation, such that RNA interference mediated knockdown resulted in enhanced activation of target genes. Together, these findings suggest a novel mechanism for NF-AT mediated gene expression and the compendium of expression data provides a quantitative platform to drive exploration of gene expression across a wide range of cell/tissue types.

Unique functions of CHK1 and WEE1 underlie synergistic anti-tumor activity upon pharmacologic inhibition.

BACKGROUND: Inhibition of kinases involved in the DNA damage response sensitizes cells to genotoxic agents by abrogating checkpoint-induced cell cycle arrest. CHK1 and WEE1 act in a pathway upstream of CDK1 to inhibit cell cycle progression in response to damaged DNA. Therapeutic targeting of either CHK1 or WEE1, in combination with chemotherapy, is under clinical evaluation. These studies examine the overlap and potential for synergy when CHK1 and WEE1 are inhibited in cancer cell models. METHODS: Small molecules MK-8776 and MK-1775 were used to selectively and potently inhibit CHK1 and WEE1, respectively. RESULTS: In vitro, the combination of MK-8776 and MK-1775 induces up to 50-fold more DNA damage than either MK-8776 or MK-1775 alone at a fixed concentration. This requires aberrant cyclin-dependent kinase activity but does not appear to be dependent on p53 status alone. Furthermore, DNA damage takes place primarily in S-phase cells, implying disrupted DNA replication. When dosed together, the combination of MK-8776 and MK-1775 induced more intense and more durable DNA damage as well as anti-tumor efficacy than either MK-8776 or MK-1775 dosed alone. DNA damage induced by the combination was detected in up to 40% of cells in a treated xenograft tumor model. CONCLUSIONS: These results highlight the roles of WEE1 and CHK1 in maintaining genomic integrity. Importantly, the strong synergy observed upon inhibition of both kinases suggests unique yet complimentary anti-tumor effects of WEE1 and CHK1 inhibition. This demonstration of DNA double strand breaks in the absence of a DNA damaging chemotherapeutic provides preclinical rationale for combining WEE1 and CHK1 inhibitors as a cancer treatment regimen.

LRRK2 is involved in the IFN-gamma response and host response to pathogens.

LRRK2 was previously identified as a defective gene in Parkinson's disease, and it is also located in a risk region for Crohn's disease. In this study, we aim to determine whether LRRK2 could be involved in immune responses. We show that LRRK2 expression is enriched in human immune cells. LRRK2 is an IFN-γ target gene, and its expression increased in intestinal tissues upon Crohn's disease inflammation. In inflamed intestinal tissues, LRRK2 is detected in the lamina propria macrophages, B-lymphocytes, and CD103-positive dendritic cells. Furthermore, LRRK2 expression enhances NF-κB-dependent transcription, suggesting its role in immune response signaling. Endogenous LRRK2 rapidly translocates near bacterial membranes, and knockdown of LRRK2 interferes with reactive oxygen species production during phagocytosis and bacterial killing. These observations indicate that LRRK2 is an IFN-γ target gene, and it might be involved in signaling pathways relevant to Crohn's disease pathogenesis.

An integrative genomics approach identifies Hypoxia Inducible Factor-1 (HIF-1)-target genes that form the core response to hypoxia.

The transcription factor Hypoxia-inducible factor 1 (HIF-1) plays a central role in the transcriptional response to oxygen flux. To gain insight into the molecular pathways regulated by HIF-1, it is essential to identify the downstream-target genes. We report here a strategy to identify HIF-1-target genes based on an integrative genomic approach combining computational strategies and experimental validation. To identify HIF-1-target genes microarrays data sets were used to rank genes based on their differential response to hypoxia. The proximal promoters of these genes were then analyzed for the presence of conserved HIF-1-binding sites. Genes were scored and ranked based on their response to hypoxia and their HIF-binding site score. Using this strategy we recovered 41% of the previously confirmed HIF-1-target genes that responded to hypoxia in the microarrays and provide a catalogue of predicted HIF-1 targets. We present experimental validation for ANKRD37 as a novel HIF-1-target gene. Together these analyses demonstrate the potential to recover novel HIF-1-target genes and the discovery of mammalian-regulatory elements operative in the context of microarray data sets.

Transcriptome profiling and network analysis of genetically hypertensive mice identifies potential pharmacological targets of hypertension.

Hypertension is a condition with major cardiovascular and renal complications, affecting nearly a billion patients worldwide. Few validated gene targets are available for pharmacological intervention, so there is a need to identify new biological pathways regulating blood pressure and containing novel targets for treatment. The genetically hypertensive "blood pressure high" (BPH), normotensive "blood pressure normal" (BPN), and hypotensive "blood pressure low" (BPL) inbred mouse strains are an ideal system to study differences in gene expression patterns that may represent such biological pathways. We profiled gene expression in liver, heart, kidney, and aorta from BPH, BPN, and BPL mice and determined which biological processes are enriched in observed organ-specific signatures. As a result, we identified multiple biological pathways linked to blood pressure phenotype that could serve as a source of candidate genes causal for hypertension. To distinguish in the kidney signature genes whose differential expression pattern may cause changes in blood pressure from those genes whose differential expression pattern results from changes in blood pressure, we integrated phenotype-associated genes into Genetic Bayesian networks. The integration of data from gene expression profiling and genetics networks is a valuable approach to identify novel potential targets for the pharmacological treatment of hypertension.

HIF-1alpha and HIF-2alpha have divergent roles in colon cancer.

Hypoxia-inducible factor (HIF)-1 and HIF-2 are heterodimeric transcription factors that mediate the cellular response to hypoxia. Their key regulatory subunits, HIF-1alpha and HIF-2alpha, are induced similarly by hypoxia, but their functional roles in cancer may be distinct and isoform-specific. SW480 colon cancer cells with stable expression of siRNA to HIF-1alpha or HIF-2alpha or both were established. HIF-1alpha-deficient cells displayed lower rates of proliferation and migration, but HIF-2alpha-deficient cells exhibited enhanced anchorage independent growth in a soft agar assay. Xenograft studies revealed that HIF-1alpha deficiency inhibited overall tumor growth, whereas deficiency of HIF-2alpha stimulated tumor growth. In human colon cancer tissues, expression of HIF-1alpha and to a lesser extent, HIF-2alpha, was linked to upregulation of VEGF and tumor angiogenesis. However, loss of expression of HIF-2alpha but not HIF-1alpha was strongly correlated with advanced tumor stage. DNA microarray analysis identified distinct sets of HIF-1alpha and HIF-2alpha target genes that may explain these phenotypic differences. Collectively, these findings suggest that HIF isoforms may have differing cellular functions in colon cancer. In particular, HIF-1alpha promoted the growth of SW480 colon cancer cells but HIF-2alpha appeared to restrain growth. Consequently, therapeutic approaches that target HIF may need to consider these isoform-specific properties.

SPTLC1 binds ABCA1 to negatively regulate trafficking and cholesterol efflux activity of the transporter.

ABCA1 transport of cholesterol and phospholipids to nascent HDL particles plays a central role in lipoprotein metabolism and macrophage cholesterol homeostasis. ABCA1 activity is regulated both at the transcriptional level and at the post-translational level. To explore mechanisms involved in the post-translational regulation of the transporter, we have used affinity purification and mass spectrometry to identify proteins that bind ABCA1 and influence its activity. Previously, we demonstrated that an interaction between beta1-syntrophin stimulated ABCA1 activity, at least in part, be slowing the degradation of the transporter. This work demonstrates that one subunit of the serine palmitoyltransferase enzyme, SPTLC1, but not subunit 2 (SPTLC2), is copurified with ABCA1 and negatively regulates its function. In human THP-I macrophages and in mouse liver, the ABCA1-SPTLC1 complex was detected by co-immunoprecipitation, demonstrating that the interaction occurs in cellular settings where ABCA1 activity is critical for HDL genesis. Pharmacologic inhibition of SPTLC1 with myriocin, which resulted in the disruption of the SPTLC1-ABCA1 complex, and siRNA knockdown of SPTLC1 expression both stimulated ABCA1 efflux by nearly 60% ( p < 0.05). In contrast, dominant-negative mutants of SPTLC1 inhibited ABCA1 efflux, indicating that a reduced level of sphingomyelin synthesis could not explain the effect of myriocin on ABCA1 activity. In 293 cells, the SPTLC1 inhibition of ABCA1 activity led to the blockade of the exit of ABCA1 from the endoplasmic reticulum. In contrast, myriocin treatment of macrophages increased the level of cell surface ABCA1. In composite, these results indicate that the physical interaction of ABCA1 and SPTLC1 results in reduction of ABCA1 activity and that inhibition of this interaction produces enhanced cholesterol efflux.

Regionalized GC content of template DNA as a predictor of PCR success.

A set of 1438 human exons was subjected to nested PCR. The initial success rate using a standard PCR protocol required for ligation-independent cloning was 83.4%. Logistic regression analysis was conducted on 27 primer- and template-related characteristics, of which most could be ignored apart from those related to the GC content of the template. Overall GC content of the template was a good predictor for PCR success; however, specificity and sensitivity values for predicted outcome were improved to 84.3 and 94.8%, respectively, when regionalized GC content was employed. This represented a significant improvement in predictability with respect to GC content alone (P < 0.001; chi(2)) and is expected to increase in relative sensitivity as template size increases. Regionalized GC was calculated with respect to a threshold of 61% GC content and a sliding window of 21 bp across the target sequence. Fine-tuning of PCR conditions is not practicable for all target sequences whenever a large number of genes of different lengths and GC content are to be amplified in parallel, particularly if total open reading frame or domain coverage is essential for recombinant protein synthesis. Thus, the present method is proposed as a means of grouping subsets of genes possessing potentially difficult target sequences so that PCR conditions can be optimized separately in order to obtain improved outcomes.

Development of a High-Throughput Gene Expression Screen for Modulators of RAS-MAPK Signaling in a Mutant RAS Cellular Context.

The RAS-MAPK pathway controls many cellular programs, including cell proliferation, differentiation, and apoptosis. In colorectal cancers, recurrent mutations in this pathway often lead to increased cell signaling that may contribute to the development of neoplasms, thereby making this pathway attractive for therapeutic intervention. To this end, we developed a 26-member gene signature of RAS-MAPK pathway activity utilizing the Affymetrix QuantiGene Plex 2.0 reagent system and performed both primary and confirmatory gene expression-based high-throughput screens (GE-HTSs) using KRAS mutant colon cancer cells (SW837) and leveraging a highly annotated chemical library. The screen achieved a hit rate of 1.4% and was able to enrich for hit compounds that target RAS-MAPK pathway members such as MEK and EGFR. Sensitivity and selectivity performance measurements were 0.84 and 1.00, respectively, indicating high true-positive and true-negative rates. Active compounds from the primary screen were confirmed in a dose-response GE-HTS assay, a GE-HTS assay using 14 additional cancer cell lines, and an in vitro colony formation assay. Altogether, our data suggest that this GE-HTS assay will be useful for larger unbiased chemical screens to identify novel compounds and mechanisms that may modulate the RAS-MAPK pathway.

An Unbiased Oncology Compound Screen to Identify Novel Combination Strategies.

Combination drug therapy is a widely used paradigm for managing numerous human malignancies. In cancer treatment, additive and/or synergistic drug combinations can convert weakly efficacious monotherapies into regimens that produce robust antitumor activity. This can be explained in part through pathway interdependencies that are critical for cancer cell proliferation and survival. However, identification of the various interdependencies is difficult due to the complex molecular circuitry that underlies tumor development and progression. Here, we present a high-throughput platform that allows for an unbiased identification of synergistic and efficacious drug combinations. In a screen of 22,737 experiments of 583 doublet combinations in 39 diverse cancer cell lines using a 4 by 4 dosing regimen, both well-known and novel synergistic and efficacious combinations were identified. Here, we present an example of one such novel combination, a Wee1 inhibitor (AZD1775) and an mTOR inhibitor (ridaforolimus), and demonstrate that the combination potently and synergistically inhibits cancer cell growth in vitro and in vivo This approach has identified novel combinations that would be difficult to reliably predict based purely on our current understanding of cancer cell biology. Mol Cancer Ther; 15(6); 1155-62. ©2016 AACR.

Preclinical evaluation of the WEE1 inhibitor MK-1775 as single-agent anticancer therapy.

Inhibition of the DNA damage checkpoint kinase WEE1 potentiates genotoxic chemotherapies by abrogating cell-cycle arrest and proper DNA repair. However, WEE1 is also essential for unperturbed cell division in the absence of extrinsic insult. Here, we investigate the anticancer potential of a WEE1 inhibitor, independent of chemotherapy, and explore a possible cellular context underlying sensitivity to WEE1 inhibition. We show that MK-1775, a potent and selective ATP-competitive inhibitor of WEE1, is cytotoxic across a broad panel of tumor cell lines and induces DNA double-strand breaks. MK-1775-induced DNA damage occurs without added chemotherapy or radiation in S-phase cells and relies on active DNA replication. At tolerated doses, MK-1775 treatment leads to xenograft tumor growth inhibition or regression. To begin addressing potential response markers for MK-1775 monotherapy, we focused on PKMYT1, a kinase functionally related to WEE1. Knockdown of PKMYT1 lowers the EC(50) of MK-1775 by five-fold but has no effect on the cell-based response to other cytotoxic drugs. In addition, knockdown of PKMYT1 increases markers of DNA damage, γH2AX and pCHK1(S345), induced by MK-1775. In a post hoc analysis of 305 cell lines treated with MK-1775, we found that expression of PKMYT1 was below average in 73% of the 33 most sensitive cell lines. Our findings provide rationale for WEE1 inhibition as a potent anticancer therapy independent of a genotoxic partner and suggest that low PKMYT1 expression could serve as an enrichment biomarker for MK-1775 sensitivity.

A functional genomic screen identifies cellular cofactors of hepatitis C virus replication.

Hepatitis C virus (HCV) chronically infects 3% of the world's population, and complications from HCV are the leading indication for liver transplantation. Given the need for better anti-HCV therapies, one strategy is to identify and target cellular cofactors of the virus lifecycle. Using a genome-wide siRNA library, we identified 96 human genes that support HCV replication, with a significant number of them being involved in vesicle organization and biogenesis. Phosphatidylinositol 4-kinase PI4KA and multiple subunits of the COPI vesicle coat complex were among the genes identified. Consistent with this, pharmacologic inhibitors of COPI and PI4KA blocked HCV replication. Targeting hepcidin, a peptide critical for iron homeostasis, also affected HCV replication, which may explain the known dysregulation of iron homeostasis in HCV infection. The host cofactors for HCV replication identified in this study should serve as a useful resource in delineating new targets for anti-HCV therapies.

Identification of host proteins required for HIV infection through a functional genomic screen.

HIV-1 exploits multiple host proteins during infection. We performed a large-scale small interfering RNA screen to identify host factors required by HIV-1 and identified more than 250 HIV-dependency factors (HDFs). These proteins participate in a broad array of cellular functions and implicate new pathways in the viral life cycle. Further analysis revealed previously unknown roles for retrograde Golgi transport proteins (Rab6 and Vps53) in viral entry, a karyopherin (TNPO3) in viral integration, and the Mediator complex (Med28) in viral transcription. Transcriptional analysis revealed that HDF genes were enriched for high expression in immune cells, suggesting that viruses evolve in host cells that optimally perform the functions required for their life cycle. This effort illustrates the power with which RNA interference and forward genetics can be used to expose the dependencies of human pathogens such as HIV, and in so doing identify potential targets for therapy.

Analysis of high throughput protein expression in Escherichia coli.

The ability to efficiently produce hundreds of proteins in parallel is the most basic requirement of many aspects of proteomics. Overcoming the technical and financial barriers associated with high throughput protein production is essential for the development of an experimental platform to query and browse the protein content of a cell (e.g. protein and antibody arrays). Proteins are inherently different one from another in their physicochemical properties; therefore, no single protocol can be expected to successfully express most of the proteins. Instead of optimizing a protocol to express a specific protein, we used sequence analysis tools to estimate the probability of a specific protein to be expressed successfully using a given protocol, thereby avoiding a priori proteins with a low success probability. A set of 547 proteins, to be used for antibody production and selection, was expressed in Escherichia coli using a high throughput protein production pipeline. Protein properties derived from sequence alone were correlated to successful expression, and general guidelines are given to increase the efficiency of similar pipelines. A second set of 68 proteins was expressed to investigate the link between successful protein expression and inclusion body formation. More proteins were expressed in inclusion bodies; however, the formation of inclusion bodies was not a requirement for successful expression.