SURFMAP: A Software for Mapping in Two Dimensions Protein Surface Features.

Molecular cartography using two-dimensional (2D) representation of protein surfaces has been shown to be very promising for protein surface analysis. Here, we present SURFMAP, a free standalone and easy-to-use software that enables the fast and automated 2D projection of either predefined features of protein surface (i.e., electrostatic potential, hydrophobicity, stickiness, and surface relief) or any descriptor encoded in the temperature factor column of a PDB file. SURFMAP proposes three different "equal-area" projections that have the advantage of preserving the area measures. It provides the user with (i) 2D maps that enable the easy and visual analysis of protein surface features of interest and (ii) maps in a text file format allowing the fast and straightforward quantitative comparison of 2D maps of homologous proteins.

Protein Interaction Energy Landscapes are Shaped by Functional and also Non-functional Partners.

In the crowded cell, a strong selective pressure operates on the proteome to limit the competition between functional and non-functional protein-protein interactions. We developed an original theoretical framework in order to interrogate how this competition constrains the behavior of proteins with respect to their partners or random encounters. Our theoretical framework relies on a two-dimensional (2D) representation of interaction energy landscapes, with 2D energy maps, which reflect in a synthetic way the spatial distribution of the interaction propensity of a protein surface for another protein. We realized the interaction propensity mapping of proteins' surfaces in interaction with functional and arbitrary partners and asked whether the distribution of their interaction propensity is conserved during evolution. Therefore, we performed several thousands of cross-docking simulations to systematically characterize the energy landscapes of 103 proteins interacting with different sets of homologs, corresponding to their functional partner's family or arbitrary protein families. Then, we systematically compared the energy maps resulting from the docking of each protein with the different protein families of the dataset. Strikingly, we show that the interaction propensity not only of the binding sites but also of the rest of the surface is conserved for docking partners belonging to the same protein family. Interestingly, this observation holds for docked proteins corresponding to true but also arbitrary partners. Our theoretical framework enables the characterization of the energy behavior of a protein in interaction with hundreds of proteins and opens the way for the characterization of the behavior of proteins in a specific environment.

Functional study of the Hap4-like genes suggests that the key regulators of carbon metabolism HAP4 and oxidative stress response YAP1 in yeast diverged from a common ancestor.

The transcriptional regulator HAP4, induced by respiratory substrates, is involved in the balance between fermentation and respiration in S. cerevisiae. We identified putative orthologues of the Hap4 protein in all ascomycetes, based only on a conserved sixteen amino acid-long motif. In addition to this motif, some of these proteins contain a DNA-binding motif of the bZIP type, while being nonetheless globally highly divergent. The genome of the yeast Hansenula polymorpha contains two HAP4-like genes encoding the protein HpHap4-A which, like ScHap4, is devoid of a bZIP motif, and HpHap4-B which contains it. This species has been chosen for a detailed examination of their respective properties. Based mostly on global gene expression studies performed in the S. cerevisiae HAP4 disruption mutant (ScΔhap4), we show here that HpHap4-A is functionally equivalent to ScHap4, whereas HpHap4-B is not. Moreover HpHAP4-B is able to complement the H2O2 hypersensitivity of the ScYap1 deletant, YAP1 being, in S. cerevisiae, the main regulator of oxidative stress. Finally, a transcriptomic analysis performed in the ScΔyap1 strain overexpressing HpHAP4-B shows that HpHap4-B acts both on oxidative stress response and carbohydrate metabolism in a manner different from both ScYap1 and ScHap4. Deletion of these two genes in their natural host, H. polymorpha, confirms that HpHAP4-A participates in the control of the fermentation/respiration balance, while HpHAP4-B is involved in oxidative stress since its deletion leads to hypersensitivity to H2O2. These data, placed in an evolutionary context, raise new questions concerning the evolution of the HAP4 transcriptional regulation function and suggest that Yap1 and Hap4 have diverged from a unique regulatory protein in the fungal ancestor.

A head and neck cancer tumor response-specific gene signature for cisplatin, 5-fluorouracil induction chemotherapy fails with added taxanes.

BACKGROUND: It is a major clinical challenge to predict which patients, with advanced stage head and neck squamous cell carcinoma, will not exhibit a reduction in tumor size following induction chemotherapy in order to avoid toxic effects of ineffective chemotherapy and delays for instituting other therapeutic options. Further, it is of interest to know to what extent a gene signature, which identifies patients with tumors that will not respond to a particular induction chemotherapy, is applicable when additional chemotherapeutic agents are added to the regimen. METHODOLOGY/PRINCIPAL FINDINGS: To identify genes that predict tumor resistance to induction with cisplatin/5-fluorouracil (PF) or PF and a taxane, we analyzed patient tumor biopsies with whole genome microarrays and quantitative reverse transcriptase-PCR (TLDA) cards. A leave one out cross-validation procedure allowed evaluation of the prediction tool. A ten-gene microarray signature correctly classified 12/13 responders and 7/10 non-responders to PF (92% specificity, 82.6% accuracy). TLDA analysis (using the same classifier) of the patients correctly classified 12/12 responders and 8/10 non-responders (100% specificity, 90.9% accuracy). Further, TLDA analysis correctly predicted the response of 5 new patients and, overall, 12/12 responders and 13/15 non-responders (100% specificity, 92.6% accuracy). The protein products of the genes constituting the signature physically associate with 27 other proteins, involved in regulating gene expression, constituting an interaction network. In contrast, TLDA-based prediction (with the same gene signature) of responses to induction with PF and either of two taxanes was poor (0% specificity, 25% accuracy and 33.3% specificity, 25% accuracy). CONCLUSIONS/SIGNIFICANCE: Successful transfer of the microarray-based gene signature to an independent, PCR-based technology suggests that TLDA-based signatures could be a useful hospital-based technology for determining therapeutic options. Although highly specific for tumor responses to PF induction, the gene signature is unsuccessful when taxanes are added. The results illustrate the subtlety in developing "personalized medicine".

Functional study of genes essential for autogamy and nuclear reorganization in Paramecium.

Like all ciliates, Paramecium tetraurelia is a unicellular eukaryote that harbors two kinds of nuclei within its cytoplasm. At each sexual cycle, a new somatic macronucleus (MAC) develops from the germ line micronucleus (MIC) through a sequence of complex events, which includes meiosis, karyogamy, and assembly of the MAC genome from MIC sequences. The latter process involves developmentally programmed genome rearrangements controlled by noncoding RNAs and a specialized RNA interference machinery. We describe our first attempts to identify genes and biological processes that contribute to the progression of the sexual cycle. Given the high percentage of unknown genes annotated in the P. tetraurelia genome, we applied a global strategy to monitor gene expression profiles during autogamy, a self-fertilization process. We focused this pilot study on the genes carried by the largest somatic chromosome and designed dedicated DNA arrays covering 484 genes from this chromosome (1.2% of all genes annotated in the genome). Transcriptome analysis revealed four major patterns of gene expression, including two successive waves of gene induction. Functional analysis of 15 upregulated genes revealed four that are essential for vegetative growth, one of which is involved in the maintenance of MAC integrity and another in cell division or membrane trafficking. Two additional genes, encoding a MIC-specific protein and a putative RNA helicase localizing to the old and then to the new MAC, are specifically required during sexual processes. Our work provides a proof of principle that genes essential for meiosis and nuclear reorganization can be uncovered following genome-wide transcriptome analysis.

Mutations in the Saccharomyces cerevisiae kinase Cbk1p lead to a fertility defect that can be suppressed by the absence of Brr1p or Mpt5p (Puf5p), proteins involved in RNA metabolism.

In Saccharomyces cerevisiae the protein kinase Cbk1p is a member of the regulation of Ace2p and cellular morphogenesis (RAM) network that is involved in cell separation after cytokinesis, cell integrity, and cell polarity. In cell separation, the RAM network promotes the daughter cell-specific localization of the transcription factor Ace2p, resulting in the asymmetric transcription of genes whose products are necessary to digest the septum joining the mother and the daughter cell. RAM and SSD1 play a role in the maintenance of cell integrity. In the presence of a wild-type SSD1 gene, deletion of any RAM component causes cell lysis. We show here that some mutations of CBK1 also lead to a reduced fertility and a reduced expression of some of the mating type-specific genes. As polarized growth is an integral part of the mating process, we have isolated suppressors of the fertility defect. Among these, mutations in BRR1 or MPT5 lead to a restoration of fertility and a more-or-less pronounced restoration of polarity; they also show genetic interactions with SSD1. Our experiments reveal a multilayered system controlling aspects of cell separation, cell integrity, mating, and polarized growth.

Multiple defects in the respiratory chain lead to the repression of genes encoding components of the respiratory chain and TCA cycle enzymes.

Respiratory complexes III, IV and V are formed by components of both nuclear and mitochondrial origin and are embedded in the inner mitochondrial membrane. Their assembly requires the auxiliary factor Oxa1, and the absence of this protein has severe consequences on these three major respiratory chain enzymes. We have studied, in the yeast Saccharomyces cerevisiae, the effect of the loss of Oxa1 function and of other respiratory defects on the expression of nuclear genes encoding components of the respiratory complexes and tricarboxylic acid cycle enzymes. We observed that the concomitant decrease in the level of two respiratory enzymes, complexes III and IV, led to their repression. These genes are known targets of the transcriptional activator complex Hap2/3/4/5 that plays a central role in the reprogramming of yeast metabolism when cells switch from a fermenting, glucose-repressed state to a respiring, derepressed state. We found that the Hap4 protein, the regulatory subunit of the transcriptional complex, was present at a lower level in the oxa1 mutants whereas no change in HAP4 transcript level was observed, suggesting a posttranscriptional modulation. In addition, an altered mitochondrial morphology was observed in mutants with decreased expression of Hap2/3/4/5 target genes. We suggest that the aberrant mitochondrial morphology, presumably caused by the severely decreased level of at least two respiratory enzymes, might be part of the signalling pathway linking the mitochondrial defect and Hap2/3/4/5.

Response of human renal tubular cells to cyclosporine and sirolimus: a toxicogenomic study.

The molecular mechanisms involved in the potentially nephrotoxic response of tubular cells to immunosuppressive drugs remain poorly understood. Transcriptional profiles of human proximal tubular cells exposed to cyclosporine A (CsA), sirolimus (SRL) or their combination, were established using oligonucleotide microarrays. Hierarchical clustering of genes implicated in fibrotic processes showed a clear distinction between expression profiles with CsA and CsA+SRL treatments on the one hand and SRL treatment on the other. Functional analysis found that CsA and CsA+SRL treatments preferentially alter biological processes located at the cell membrane, such as ion transport or signal transduction, whereas SRL modifies biological processes within the nucleus and related to transcriptional activity. Genome wide expression analysis suggested that CsA may induce an endoplasmic reticulum (ER) stress in tubular cells in vitro. Moreover we found that CsA exposure in vivo is associated with the upregulation of the ER stress marker BIP in kidney transplant biopsies. In conclusion, this toxicogenomic study highlights the molecular interaction networks that may contribute to the tubular response to CsA and SRL. These results may also offer a new working hypothesis for future research in the field of CsA nephrotoxicity. Further studies are needed to evaluate if ER stress detection in tubular cells in human biopsies can predict CsA nephrotoxicity.

Expression ratio evaluation in two-colour microarray experiments is significantly improved by correcting image misalignment.

MOTIVATION: Two-colour microarrays are widely used to perform transcriptome analysis. In most cases, it appears that the 'red' and 'green' images resulting from the scan of a microarray slide are slightly shifted one with respect to the other. To increase the robustness of the measurement of the fluorescent emission intensities, multiple acquisitions with the same or different PMT gains can be used. In these cases, a systematic correction of image shift is required. RESULTS: To accurately detect this shift, we first developed an approach using cross-correlation. Second, we evaluated the most appropriate interpolation method to be used to derive the registered image. Then, we quantified the effects of image shifts on spot quality, using two different quality estimators. Finally, we measured the benefits associated with a systematic image registration. In this study, we demonstrate that registering the two images prior to data extraction provides a more reliable estimate of the two colours' ratio and thus increases the accuracy of measurements of variations in gene expression. AVAILABILITY: http://bioinfome.cgm.cnrs-gif.fr/.