Functional and topological characterization of transcriptional cooperativity in yeast

This work was supported by grants from Spanish Ministry of Science andInnovation (MICINN) BIO2011-22568 & BIO2008-205.

Improving the prediction of the functional impact of cancer mutations by baseline tolerance transformation

High-throughput prioritization of cancer-causing mutations (drivers) is a key challenge of cancer genome projects, due to the number of somatic variants detected in tumors. One important step in this task is to assess the functional impact of tumor somatic mutations. A number of computational methods have been employed for that purpose, although most were originally developed to distinguish disease-related nonsynonymous single nucleotide variants (nsSNVs) from polymorphisms. Our new method, transformed Functional Impact score for Cancer (transFIC), improves the assessment of the functional impact of tumor nsSNVs by taking into account the baseline tolerance of genes to functional variants.

Prediction of cross-recognition of peptide-HLA A2 by Melan-A-specific cytotoxic T lymphocytes using three-dimensional quantitative structure-activity relationships.

The cross-recognition of peptides by cytotoxic T lymphocytes is a key element in immunology and in particular in peptide based immunotherapy. Here we develop three-dimensional (3D) quantitative structure-activity relationships (QSARs) to predict cross-recognition by Melan-A-specific cytotoxic T lymphocytes of peptides bound to HLA A*0201 (hereafter referred to as HLA A2). First, we predict the structure of a set of self- and pathogen-derived peptides bound to HLA A2 using a previously developed ab initio structure prediction approach [Fagerberg et al., J. Mol. Biol., 521-46 (2006)]. Second, shape and electrostatic energy calculations are performed on a 3D grid to produce similarity matrices which are combined with a genetic neural network method [So et al., J. Med. Chem., 4347-59 (1997)] to generate 3D-QSAR models. The models are extensively validated using several different approaches. During the model generation, the leave-one-out cross-validated correlation coefficient (q (2)) is used as the fitness criterion and all obtained models are evaluated based on their q (2) values. Moreover, the best model obtained for a partitioned data set is evaluated by its correlation coefficient (r = 0.92 for the external test set). The physical relevance of all models is tested using a functional dependence analysis and the robustness of the models obtained for the entire data set is confirmed using y-randomization. Finally, the validated models are tested for their utility in the setting of rational peptide design: their ability to discriminate between peptides that only contain side chain substitutions in a single secondary anchor position is evaluated. In addition, the predicted cross-recognition of the mono-substituted peptides is confirmed experimentally in chromium-release assays. These results underline the utility of 3D-QSARs in peptide mimetic design and suggest that the properties of the unbound epitope are sufficient to capture most of the information to determine the cross-recognition.

Functional analysis of Plasmodium vivax VIR proteins reveals different subcellular localizations and cytoadherence to the ICAM-1 endothelial receptor.

The subcellular localization and function of variant subtelomeric multigene families in Plasmodium vivax remain vastly unknown. Among them, the vir superfamily is putatively involved in antigenic variation and in mediating adherence to endothelial receptors. In the absence of a continuous in vitro culture system for P. vivax, we have generated P. falciparum transgenic lines expressing VIR proteins to infer location and function. We chose three proteins pertaining to subfamilies A (VIR17), C (VIR14) and D (VIR10), with domains and secondary structures that predictably traffic these proteins to different subcellular compartments. Here, we showed that VIR17 remained inside the parasite and around merozoites, whereas VIR14 and VIR10 were exported to the membrane of infected red blood cells (iRBCs) in an apparent independent pathway of Maurer's clefts. Remarkably, VIR14 was exposed at the surface of iRBCs and mediated adherence to different endothelial receptors expressed in CHO cells under static conditions. Under physiological flow conditions, however, cytoadherence was only observed to ICAM-1, which was the only receptor whose adherence was specifically and significantly inhibited by antibodies against conserved motifs of VIR proteins. Immunofluorescence studies using these antibodies also showed different subcellular localizations of VIR proteins in P. vivax-infected reticulocytes from natural infections. These data suggest that VIR proteins are trafficked to different cellular compartments and functionally demonstrates that VIR proteins can specifically mediate cytoadherence to the ICAM-1 endothelial receptor.

Predicting current and future spatial community patterns of plant functional traits

Community-level patterns of functional traits relate to community assembly and ecosystem functioning. By modelling the changes of different indices describing such patterns - trait means, extremes and diversity in communities - as a function of abiotic gradients, we could understand their drivers and build projections of the impact of global change on the functional components of biodiversity. We used five plant functional traits (vegetative height, specific leaf area, leaf dry matter content, leaf nitrogen content and seed mass) and non-woody vegetation plots to model several indices depicting community-level patterns of functional traits from a set of abiotic environmental variables (topographic, climatic and edaphic) over contrasting environmental conditions in a mountainous landscape. We performed a variation partitioning analysis to assess the relative importance of these variables for predicting patterns of functional traits in communities, and projected the best models under several climate change scenarios to examine future potential changes in vegetation functional properties. Not all indices of trait patterns within communities could be modelled with the same level of accuracy: the models for mean and extreme values of functional traits provided substantially better predictive accuracy than the models calibrated for diversity indices. Topographic and climatic factors were more important predictors of functional trait patterns within communities than edaphic predictors. Overall, model projections forecast an increase in mean vegetation height and in mean specific leaf area following climate warming. This trend was important at mid elevation particularly between 1000 and 2000 m asl. With this study we showed that topographic, climatic and edaphic variables can successfully model descriptors of community-level patterns of plant functional traits such as mean and extreme trait values. However, which factors determine the diversity of functional traits in plant communities remains unclear and requires more investigations.

Gene-disease network analysis reveals functional modules in mendelian, complex and environmental diseases

AbstractBACKGROUND: Scientists have been trying to understand the molecular mechanisms of diseases to design preventive and therapeutic strategies for a long time. For some diseases, it has become evident that it is not enough to obtain a catalogue of the disease-related genes but to uncover how disruptions of molecular networks in the cell give rise to disease phenotypes. Moreover, with the unprecedented wealth of information available, even obtaining such catalogue is extremely difficult.PRINCIPAL FINDINGS: We developed a comprehensive gene-disease association database by integrating associations from several sources that cover different biomedical aspects of diseases. In particular, we focus on the current knowledge of human genetic diseases including mendelian, complex and environmental diseases. To assess the concept of modularity of human diseases, we performed a systematic study of the emergent properties of human gene-disease networks by means of network topology and functional annotation analysis. The results indicate a highly shared genetic origin of human diseases and show that for most diseases, including mendelian, complex and environmental diseases, functional modules exist. Moreover, a core set of biological pathways is found to be associated with most human diseases. We obtained similar results when studying clusters of diseases, suggesting that related diseases might arise due to dysfunction of common biological processes in the cell.CONCLUSIONS: For the first time, we include mendelian, complex and environmental diseases in an integrated gene-disease association database and show that the concept of modularity applies for all of them. We furthermore provide a functional analysis of disease-related modules providing important new biological insights, which might not be discovered when considering each of the gene-disease association repositories independently. Hence, we present a suitable framework for the study of how genetic and environmental factors, such as drugs, contribute to diseases.AVAILABILITY: The gene-disease networks used in this study and part of the analysis are available at http://ibi.imim.es/DisGeNET/DisGeNETweb.html#Download

Rapid Perturbation in Viremia Levels Drives Increases in Functional Avidity of HIV-specific CD8 T Cells.

The factors determining the functional avidity and its relationship with the broad heterogeneity of antiviral T cell responses remain partially understood. We investigated HIV-specific CD8 T cell responses in 85 patients with primary HIV infection (PHI) or chronic (progressive and non-progressive) infection. The functional avidity of HIV-specific CD8 T cells was not different between patients with progressive and non-progressive chronic infection. However, it was significantly lower in PHI patients at the time of diagnosis of acute infection and after control of virus replication following one year of successful antiretroviral therapy. High-avidity HIV-specific CD8 T cells expressed lower levels of CD27 and CD28 and were enriched in cells with an exhausted phenotype, i.e. co-expressing PD-1/2B4/CD160. Of note, a significant increase in the functional avidity of HIV-specific CD8 T cells occurred in early-treated PHI patients experiencing a virus rebound after spontaneous treatment interruption. This increase in functional avidity was associated with the accumulation of PD-1/2B4/CD160 positive cells, loss of polyfunctionality and increased TCR renewal. The increased TCR renewal may provide the mechanistic basis for the generation of high-avidity HIV-specific CD8 T cells. These results provide insights on the relationships between functional avidity, viremia, T-cell exhaustion and TCR renewal of antiviral CD8 T cell responses.

A combined computational and functional approach identifies new residues involved in pH-dependent gating of ASIC1a.

Acid-sensing ion channels (ASICs) are key receptors for extracellular protons. These neuronal nonvoltage-gated Na(+) channels are involved in learning, the expression of fear, neurodegeneration after ischemia, and pain sensation. We have applied a systematic approach to identify potential pH sensors in ASIC1a and to elucidate the mechanisms by which pH variations govern ASIC gating. We first calculated the pK(a) value of all extracellular His, Glu, and Asp residues using a Poisson-Boltzmann continuum approach, based on the ASIC three-dimensional structure, to identify candidate pH-sensing residues. The role of these residues was then assessed by site-directed mutagenesis and chemical modification, combined with functional analysis. The localization of putative pH-sensing residues suggests that pH changes control ASIC gating by protonation/deprotonation of many residues per subunit in different channel domains. Analysis of the function of residues in the palm domain close to the central vertical axis of the channel allowed for prediction of conformational changes of this region during gating. Our study provides a basis for the intrinsic ASIC pH dependence and describes an approach that can also be applied to the investigation of the mechanisms of the pH dependence of other proteins.

Investigation of high-resolution functional magnetic resonance imaging by means of surface and array radiofrequency coils at 7 T.

In this investigation, high-resolution, 1x1x1-mm(3) functional magnetic resonance imaging (fMRI) at 7 T is performed using a multichannel array head coil and a surface coil approach. Scan geometry was optimized for each coil separately to exploit the strengths of both coils. Acquisitions with the surface coil focused on partial brain coverage, while whole-brain coverage fMRI experiments were performed with the array head coil. BOLD sensitivity in the occipital lobe was found to be higher with the surface coil than with the head array, suggesting that restriction of signal detection to the area of interest may be beneficial for localized activation studies. Performing independent component analysis (ICA) decomposition of the fMRI data, we consistently detected BOLD signal changes and resting state networks. In the surface coil data, a small negative BOLD response could be detected in these resting state network areas. Also in the data acquired with the surface coil, two distinct components of the positive BOLD signal were consistently observed. These two components were tentatively assigned to tissue and venous signal changes.

N-Glycosylation and conserved cysteine residues in RAMP3 play a critical role for the functional expression of CRLR/RAMP3 adrenomedullin receptor.

The calcitonin receptor-like receptor (CRLR) and receptor activity modifying protein-3 (RAMP3) can assemble into a CRLR/RAMP3 heterodimeric receptor that exhibits the characteristics of a high affinity adrenomedullin receptor. RAMP3 participates in adrenomedullin (AM) binding via its extracellular N-terminus characterized by the presence of six highly conserved cysteine residues and four N-glycosylation consensus sites. Here, we assessed the usage of these conserved residues in cotranslational modifications of RAMP3 and addressed their role in functional expression of the CRLR/RAMP3 receptor. Using a Xenopus oocyte expression system, we show that (i) RAMP3 is assembled with CRLR as a multiple N-glycosylated species in which two, three, or four consensus sites are used; (ii) elimination of all N-glycans in RAMP3 results in a significant inhibition of receptor [(125)I]AM binding and an increase in the EC(50) value for AM; (iii) several lines of indirect evidence indicate that each of the six cysteines is involved in disulfide bond formation; (iv) when all cysteines are mutated to serines, RAMP3 is N-glycosylated at all four consensus sites, suggesting that disulfide bond formation inhibits N-gylcosylation; and (v) elimination of all cysteines abolishes adrenomedullin binding and leads to a complete loss of receptor function. Our data demonstrate that cotranslational modifications of RAMP3 play a critical role in the function of the CRLR/RAMP3 adrenomedullin receptor.

Numerous potentially functional but non-genic conserved sequences on human chromosome 21.

The use of comparative genomics to infer genome function relies on the understanding of how different components of the genome change over evolutionary time. The aim of such comparative analysis is to identify conserved, functionally transcribed sequences such as protein-coding genes and non-coding RNA genes, and other functional sequences such as regulatory regions, as well as other genomic features. Here, we have compared the entire human chromosome 21 with syntenic regions of the mouse genome, and have identified a large number of conserved blocks of unknown function. Although previous studies have made similar observations, it is unknown whether these conserved sequences are genes or not. Here we present an extensive experimental and computational analysis of human chromosome 21 in an effort to assign function to sequences conserved between human chromosome 21 (ref. 8) and the syntenic mouse regions. Our data support the presence of a large number of potentially functional non-genic sequences, probably regulatory and structural. The integration of the properties of the conserved components of human chromosome 21 to the rapidly accumulating functional data for this chromosome will improve considerably our understanding of the role of sequence conservation in mammalian genomes.

Absence of p53 gene mutations in a tumor panel representative of pilocytic astrocytoma diversity using a p53 functional assay.

Although p53-gene mutations occur with significant frequency in diffuse low-grade and high-grade astrocytomas, and are postulated to play an important role in tumorigenesis in these cases, the role of the p53 gene in pilocytic astrocytomas remains unclear. Published data using DNA-based assays for p53-gene analysis in these tumors have shown contradictory results in mutation frequency (0-14%). It is not known whether these heterogeneous results stem from the biological diversity of this tumor group or from technical problems. To re-evaluate p53-gene status in pilocytic tumors, we analyzed 18 tumors chosen to represent the clinical and biological heterogeneity of this tumor type with respect to anatomical location, patient age, gender, ethnic origin (Caucasian or Japanese) and the concomitant occurrence of neurofibromatosis type 1 (NF1). All primary tumors were histologically diagnosed as pilocytic astrocytoma (WHO grade I), except for one anaplastic pilocytic astrocytoma (WHO grade III) which developed in an NF1 patient and recurred as glioblastoma multiforme (WHO grade IV). p53 mutations were detected using an assay in yeast which tests the transcriptional activity of p53 proteins synthesized from tumor mRNA-derived p53-cDNA templates. None of 18 tumors, including 3 NF1-related tumors, showed p53-gene mutations between and including exons 4 and 11. We conclude that p53-gene mutations are extremely rare findings in pilocytic astrocytomas, and are absent even in those exceptional cases in which malignant progression of such tumors has occurred.

Cloning, functional expression, and chromosomal localization of the human pancreatic islet glucose-dependent insulinotropic polypeptide receptor.

Glucose-dependent insulinotropic polypeptide (GIP) is a hormone secreted by the endocrine K-cells from the duodenum that stimulates glucose-induced insulin secretion. Here, we present the molecular characterization of the human pancreatic islet GIP receptor. cDNA clones for the GIP receptor were isolated from a human pancreatic islet cDNA library. They encoded two different forms of the receptor, which differed by a 27-amino acid insertion in the COOH-terminal cytoplasmic tail. The receptor protein sequence was 81% identical to that of the rat GIP receptor. When expressed in Chinese hamster lung fibroblasts, both forms of the receptor displayed high-affinity binding for GIP (180 and 600 pmol/l). GIP binding was displaced by < 20% by 1 mumol/l glucagon, glucagon-like peptide (GLP-I)(7-36) amide, vasoactive intestinal peptide, and secretin. However exendin-4 and exendin-(9-39) at 1 mumol/l displaced binding by approximately 70 and approximately 100% at 10 mumol/l. GIP binding to both forms of the receptor induced a dose-dependent increase in intracellular cAMP levels (EC50 values of 0.6-0.8 nmol/l) but no elevation of cytoplasmic calcium concentrations. Interestingly, both exendin-4 and exendin-(9-39) were antagonists of the receptor, inhibiting GIP-induced cAMP formation by up to 60% when present at a concentration of 10 mumol/l. Finally, the physical and genetic chromosomal localization of the receptor gene was determined to be on 19q13.3, close to the ApoC2 gene. These data will help study the physiology and pathophysiology of the human GIP receptor.