The implications of alternative splicing in the ENCODE protein complement.

Alternative premessenger RNA splicing enables genes to generate more than one gene product. Splicing events that occur within protein coding regions have the potential to alter the biological function of the expressed protein and even to create new protein functions. Alternative splicing has been suggested as one explanation for the discrepancy between the number of human genes and functional complexity. Here, we carry out a detailed study of the alternatively spliced gene products annotated in the ENCODE pilot project. We find that alternative splicing in human genes is more frequent than has commonly been suggested, and we demonstrate that many of the potential alternative gene products will have markedly different structure and function from their constitutively spliced counterparts. For the vast majority of these alternative isoforms, little evidence exists to suggest they have a role as functional proteins, and it seems unlikely that the spectrum of conventional enzymatic or structural functions can be substantially extended through alternative splicing.

The activation process of the alpha1B-adrenergic receptor: potential role of protonation and hydrophobicity of a highly conserved aspartate.

In this study, a quantitative approach was used to investigate the role of D142, which belongs to the highly conserved E/DRY sequence, in the activation process of the alpha1B-adrenergic receptor (alpha1B-AR). Experimental and computer-simulated mutagenesis were performed by substituting all possible natural amino acids at the D142 site. The resulting congeneric set of proteins together with the finding that all the receptor mutants show various levels of constitutive (agonist-independent) activity enabled us to quantitatively analyze the relationships between structural/dynamic features and the extent of constitutive activity. Our results suggest that the hydrophobic/hydrophilic character of D142, which could be regulated by protonation/deprotonation of this residue, is an important modulator of the transition between the inactive (R) and active (R*) state of the alpha1B-AR. Our study represents an example of quantitative structure-activity relationship analysis of the activation process of a G protein-coupled receptor.

Identification of novel and cell type enriched cofactors of the transcription activation domain of RelA (p65 NF-kappaB).

RelA (NF-kappaB) is a transcription factor inducible by distinct stimuli in many different cell types. To find new cell type specific cofactors of NF-kappaB dependent transcription, we isolated RelA transcription activation domain binding proteins from the nuclear extracts of three different cell types. Analysis by electrophoresis and liquid chromatography tandem mass spectrometry identified several novel putative molecular partners. Some were strongly enriched in the complex formed from the nuclear extracts of specific cell types.

Unique spectrum of activity of prosimian TRIM5alpha against exogenous and endogenous retroviruses.

Lentiviruses, the genus of retrovirus that includes HIV-1, rarely endogenize. Some lemurs uniquely possess an endogenous lentivirus called PSIV ("prosimian immunodeficiency virus"). Thus, lemurs provide the opportunity to study the activity of host defense factors, such as TRIM5α, in the setting of germ line invasion. We characterized the activities of TRIM5α proteins from two distant lemurs against exogenous retroviruses and a chimeric PSIV. TRIM5α from gray mouse lemur, which carries PSIV in its genome, exhibited the narrowest restriction activity. One allelic variant of gray mouse lemur TRIM5α restricted only N-tropic murine leukemia virus (N-MLV), while a second variant restricted N-MLV and, uniquely, B-tropic MLV (B-MLV); both variants poorly blocked PSIV. In contrast, TRIM5α from ring-tailed lemur, which does not contain PSIV in its genome, revealed one of the broadest antiviral activities reported to date against lentiviruses, including PSIV. Investigation into the antiviral specificity of ring-tailed lemur TRIM5α demonstrated a major contribution of a 32-amino-acid expansion in variable region 2 (v2) of the B30.2/SPRY domain to the breadth of restriction. Data on lemur TRIM5α and the prediction of ancestral simian sequences hint at an evolutionary scenario where antiretroviral specificity is prominently defined by the lineage-specific expansion of the variable loops of B30.2/SPRY.

Genomic roles of the HCF transcriptional regulator in "Drosophila"

Abstract : Transcriptional regulation is the result of a combination of positive and negative effectors, such as transcription factors, cofactors and chromatin modifiers. During my thesis project I studied chromatin association, and transcriptional and cell cycle regulatory functions of dHCF, the Drosophila homologue of the human protein HCF-1 (host cell factor-1). The human and Drosophila HCF proteins are synthesized as large polypeptides that are cleaved into two subunits (HCFN and HCFC), which remain associated with one another by non covalent interactions. Studies in mammalian cells over the past 20 years have been devoted to understanding the cellular functions of HCF-1 and have revealed that it is a key regulator of transcription and cell cycle regulation. In human cells, HCF-1 interacts with the histone methyltransferase Set1/Ash2 and MLL/Ash2 complexes and the histone deacetylase Sin3 complex, which are involved in transcriptional activation and repression, respectively. HCF-1 is also recruited to promoters to regulate G1 -to-S phase progression during the cell cycle by the activator transcription factors E2F1 and E2F3, and by the repressor transcription factor E2F4. HCF-1 protein structure and these interactions between HCP-1 and E2F transcriptional regulator proteins are also conserved in Drosophila. In this doctoral thesis, I use proliferating Drosophila SL2 cells to study both the genomic-binding sites of dHCF, using a combination of chromatin immunoprecipitation and ultra high throughput sequencing (ChIP-seq) analysis, and dHCF regulated genes, employing RNAi and microarray expression analysis. I show that dHCF is bound to over 7500 chromosomal sites in proliferating SL2 cells, and is located at +-200 bp relative to the transcriptional start sites of about 30% of Drosophila genes. There is also a direct relationship between dHCF promoter association and promoter- associated transcriptional activity. Thus, dHCF binding levels at promoters correlated directly with transcriptional activity. In contrast, expression studies showed that dHCF appears to be involved in both transcriptional activation and repression. Analysis of dHCF-binding sites identified nine dHCF-associated motifs, four of them linked dHCF to (i) two insulator proteins, GAGA and BEAF, (ii) the E-box motif, and (iii) a degenerated TATA-box. The dHCF-associated motifs allowed the organization of the dHCF-bound genes into five biological processes: differentiation, cell cycle and gene expression, regulation of endocytosis, and cellular localization. I further show that different mechanisms regulate dHCF association with chromatin. Despite that after dHCF cleavage the dHCFN and dHCFC subunits remain associated, the two subunits showed different affinities for chromatin and differential binding to a set of tested promoters, suggesting that dHCF could target specific promoters through each of the two subunits. Moreover, in addition to the interaction between dHCF and E2F transcription factors, the dHCF binding pattern is correlated with dE2F2 genomic 4 distribution. I show that dE2F factors are necessary for recruitment of dHCF to the promoter of a set of dHCF regulated genes. Therefore dHCF, as in mammals, is involved in regulation of G1 to S phase progression in collaboration with the dE2Fs transcription factors. In addition, gene expression arrays reveal that dHCF could indirectly regulate cell cycle progression by promoting expression of genes involved in gene expression and protein synthesis, and inhibiting expression of genes involved in cell-cell adhesion. Therefore, dHCF is an evolutionary conserved protein, which binds to many specific sites of the Drosophila genome via interaction with DNA of chromatin-binding proteins to regulate the expression of genes involved in many different cellular functions. Résumé : La regulation de la transcription est le résultat des effets positifs et négatifs des facteurs de transcription, cofacteurs et protéines effectrices qui modifient la chromatine. Pendant mon projet de thèse, j'ai étudié l'association a la chromatine, ainsi que la régulation de la transcription et du cycle cellulaire par dHCF, l'homologue chez la drosophile de la protéine humaine HCF-1 (host cell factor-1). Chez 1'humain et la V drosophile, les deux protéines HCF sont synthétisées sous la forme d'un long polypeptide, qui est ensuite coupé en deux sous-unités au centre de la protéine. Les deux sous-unités restent associées ensemble grâce a des interactions non-covalentes. Des études réalisées pendant les 20 dernières années ont permit d'établir que HCF-l et un facteur clé dans la régulation de la transcription et du cycle cellulaire. Dans les cellules humaines, HCF-1 active et réprime la transcription en interagissant avec des complexes de protéines qui activent la transcription en méthylant les histones (HMT), comme par Set1/Ash2 et MLL/Ash2, et d'autres complexes qui répriment la transcription et sont responsables de la déacétylation des histones (HDAC) comme la protéine Sin3. HCF-l est aussi recruté aux promoteurs par les activateurs de la transcription E2F l et E2F3a, et par le répresseur de la transcription E2F4 pour réguler la transition entre les phases G1 et S du cycle cellulaire. La structure de HCF-1 et les interactions entre HCF-l et les régulateurs de la transcription sont conservées chez la drosophile. Pendant ma these j'ai utilisé les cellules de la drosophile, SL2 en culture, pour étudier les endroits de liaisons de HCF-l à la chromatine, grâce a immunoprecipitation de la chromatine et du séquençage de l'ADN massif ainsi que les gènes régulés par dHCF 3 grâce a la technique de RNAi et des microarrays. Mes résultats on montré que dHCF se lie à environ 7565 endroits, et estimé a 1200 paire de bases autour des sites d'initiation de la transcription de 30% des gènes de la drosophile. J 'ai observe une relation entre dHCF et le niveau de la transcription. En effet, le niveau de liaison dHCF au promoteur corrèle avec l'activité de la transcription. Cependant, mes études d'expression ont montré que dHCF est implique dans le processus d'activation et mais aussi de répression de la transcription. L'analyse des séquences d'ADN liées par dHCF a révèle neuf motifs, quatre de ces motifs ont permis d'associer dl-ICF a deux protéines isolatrices GAGA et BEAF, au motif pour les E-boxes et a une TATA-box dégénérée. Les neuf motifs associes à dHCF ont permis d'associer les gènes lies par dHCF au promoteur a cinq processus biologiques: différentiation, cycle cellulaire, expression de gènes, régulation de l'endocytosis et la localisation cellulaire, J 'ai aussi montré qu'il y a plusieurs mécanismes qui régulent l'association de dHCF a la chromatine, malgré qu'après clivage, les deux sous-unites dHCFN and dHCFC, restent associées, elles montrent différentes affinités pour la chromatine et lient différemment un group de promoteurs, les résultats suggèrent que dHCF peut se lier aux promoteurs en utilisant chacune de ses sous-unitées. En plus de l'association de dHCF avec les facteurs de transcription dE2F s, la distribution de dHCF sur le génome corrèle avec celle du facteur de transcription dE2F2. J'ai aussi montré que les dE2Fs sont nécessaires pour le recrutement de dHCF aux promoteurs d'un sous-groupe de gènes régules par dHCF. Mes résultats ont aussi montré que chez la drosophile comme chez les humains, dl-ICF est implique dans la régulation de la progression de la phase G1 a la phase S du cycle cellulaire en collaboration avec dE2Fs. D'ailleurs, les arrays d'expression ont suggéré que dHCF pourrait réguler le cycle cellulaire de façon indirecte en activant l'expression de gènes impliqués dans l'expression génique et la synthèse de protéines, et en inhibant l'expression de gènes impliqués dans l'adhésion cellulaire. En conclusion, dHCF est une protéine, conservée dans l'évolution, qui se lie spécifiquement a beaucoup d'endroits du génome de Drosophile, grâce à l'interaction avec d'autres protéines, pour réguler l'expression des gènes impliqués dans plusieurs fonctions cellulaires.

Constitutive activity and inverse agonism at the alpha1adrenoceptors.

Mutations of G protein-coupled receptors (GPCR) can increase their constitutive (agonist-independent) activity. Some of these mutations have been artificially introduced by site-directed mutagenesis, others occur spontaneously in human diseases. The alpha(1B)adrenoceptor was the first GPCR in which point mutations were shown to trigger receptor activation. This article briefly summarizes some of the findings reported in the last several years on constitutive activity of the alpha(1)adrenoceptor subtypes, the location where mutations have been found in the receptors, the spontaneous activity of native receptors in recombinant as well as physiological systems. In addition, it will highlight how the analysis of the pharmacological and molecular properties of the constitutively active adrenoceptor mutants provided an important contribution to our understanding of the molecular mechanisms underlying the mechanism of receptor activation and inverse agonism.

The malaria circumsporozoite protein has two functional domains, each with distinct roles as sporozoites journey from mosquito to mammalian host.

Plasmodium sporozoites make a remarkable journey from the mosquito midgut to the mammalian liver. The sporozoite's major surface protein, circumsporozoite protein (CSP), is a multifunctional protein required for sporozoite development and likely mediates several steps of this journey. In this study, we show that CSP has two conformational states, an adhesive conformation in which the C-terminal cell-adhesive domain is exposed and a nonadhesive conformation in which the N terminus masks this domain. We demonstrate that the cell-adhesive domain functions in sporozoite development and hepatocyte invasion. Between these two events, the sporozoite must travel from the mosquito midgut to the mammalian liver, and N-terminal masking of the cell-adhesive domain maintains the sporozoite in a migratory state. In the mammalian host, proteolytic cleavage of CSP regulates the switch to an adhesive conformation, and the highly conserved region I plays a critical role in this process. If the CSP domain architecture is altered such that the cell-adhesive domain is constitutively exposed, the majority of sporozoites do not reach their target organs, and in the mammalian host, they initiate a blood stage infection directly from the inoculation site. These data provide structure-function information relevant to malaria vaccine development.

Stable alpha-synuclein oligomers strongly inhibit chaperone activity of the Hsp70 system by weak interactions with J-domain co-chaperones.

α-Synuclein aggregation and accumulation in Lewy bodies are implicated in progressive loss of dopaminergic neurons in Parkinson disease and related disorders. In neurons, the Hsp70s and their Hsp40-like J-domain co-chaperones are the only known components of chaperone network that can use ATP to convert cytotoxic protein aggregates into harmless natively refolded polypeptides. Here we developed a protocol for preparing a homogeneous population of highly stable ^6;-sheet enriched toroid-shaped α-Syn oligomers with a diameter typical of toxic pore-forming oligomers. These oligomers were partially resistant to in vitro unfolding by the bacterial Hsp70 chaperone system (DnaK, DnaJ, GrpE). Moreover, both bacterial and human Hsp70/Hsp40 unfolding/refolding activities of model chaperone substrates were strongly inhibited by the oligomers but, remarkably, not by unstructured α-Syn monomers even in large excess. The oligomers acted as a specific competitive inhibitor of the J-domain co-chaperones, indicating that J-domain co-chaperones may preferably bind to exposed bulky misfolded structures in misfolded proteins and, thus, complement Hsp70s that bind to extended segments. Together, our findings suggest that inhibition of the Hsp70/Hsp40 chaperone system by α-Syn oligomers may contribute to the disruption of protein homeostasis in dopaminergic neurons, leading to apoptosis and tissue loss in Parkinson disease and related neurodegenerative diseases.

Microglial responsiveness as a sensitive marker for trimethyltin (TMT) neurotoxicity.

Activation of microglia is a well-documented phenomenon associated with diverse pathological conditions of the central nervous system. In order to investigate the involvement of microglial cells in the neurotoxic action of the heavy metal compound trimethyltin, three-dimensional brain cell cultures were treated during an early developmental period, using concentrations at or below the limit of cytotoxicity. Microglial cells were studied by cytochemical staining, using horseradish peroxidase-conjugated B4 isolectin (GSI-B4). In parallel, neurotoxic effects were assessed by determining the content of synaptophysin and synapsin I, both in the total homogenates and in the synaptosomal fraction of the cultures. Changes in the content of the specific growth cone protein, GAP-43, were also analyzed. It was found that low, non-cytotoxic concentrations of TMT (10(-9) to 10(-8) M) caused a significant increase in the number and/or the clustering of microglial cells. A decrease in the synaptic protein (synapsin I, synaptophysin) content was detected at 10(-8) M of TMT in synaptosomal fractions, whereas in the total homogenates, changes in synaptic proteins and GAP-43 were observed only at the cytotoxic TMT concentration (10(-6) M). Although it remains to be shown whether the microglial response is caused by direct or indirect action of TMT, the present findings show that microglial responsiveness can be detected prior to any sign of neuronal degeneration, and may serve as a sensitive indicator for heavy metal neurotoxicity in the brain.

EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization.

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.

Identification of SPLUNC1's ENaC-inhibitory domain yields novel strategies to treat sodium hyperabsorption in cystic fibrosis airways.

The epithelial sodium channel (ENaC) is responsible for Na+ and fluid absorption across colon, kidney, and airway epithelia. We have previously identified SPLUNC1 as an autocrine inhibitor of ENaC. We have now located the ENaC inhibitory domain of SPLUNC1 to SPLUNC1's N terminus, and a peptide corresponding to this domain, G22-A39, inhibited ENaC activity to a similar degree as full-length SPLUNC1 (W64;2.5 fold). However, G22-A39 had no effect on the structurally related acid-sensing ion channels, indicating specificity for ENaC. G22-A39 preferentially bound to the ^6;-ENaC subunit in a glycosylation-dependent manner. ENaC hyperactivity is contributory to cystic fibrosis (CF) lung disease. Addition of G22-A39 to CF human bronchial epithelial cultures (HBECs) resulted in an increase in airway surface liquid height from 4.2±0.6 to 7.9±0.6 _6;m, comparable to heights seen in normal HBECs, even in the presence of neutrophil elastase. Our data also indicate that the ENaC inhibitory domain of SPLUNC1 may be cleaved away from the main molecule by neutrophil elastase, which suggests that it may still be active during inflammation or neutrophilia. Furthermore, the robust inhibition of ENaC by the G22-A39 peptide suggests that this peptide may be suitable for treating CF lung disease.

The contact region between three domains of the extracellular loop of ASIC1a is critical for channel function.

Acid-sensing ion channels are members of the epithelial Na(+) channel/degenerin family. They are neuronal nonvoltage-gated Na(+) channels that are activated by extracellular acidification. In this study, we investigated the role of a highly conserved region of the extracellular part of ASIC1a that forms the contact between the finger domain, the adjacent beta-ball, and the upper palm domain in ASIC1a. The finger domain contributes to the pH-dependent gating and is linked via this contact zone to the rest of the protein. We found that mutation to Cys of residues in this region led to decreased channel expression and current amplitudes. Exposure of the engineered Cys residues to Cd(2+) or to charged methane thiosulfonate sulfhydryl reagents further reduced current amplitudes. This current inhibition was not due to changes in acid-sensing ion channel pH dependence or unitary conductance and was likely due to a decrease of the probability of channel opening. For some mutants, the effect of sulfhydryl reagents depended on the pH of exposure in the range 7.4 to 6.8, suggesting that this zone undergoes conformational changes during inactivation. Our study identifies a region in ASIC1a whose integrity is required for normal channel function.

Structural assessment of single amino acid mutations: application to TP53 function.

Single amino acid substitution is the type of protein alteration most related to human diseases. Current studies seek primarily to distinguish neutral mutations from harmful ones. Very few methods offer an explanation of the final prediction result in terms of the probable structural or functional effect on the protein. In this study, we describe the use of three novel parameters to identify experimentally-verified critical residues of the TP53 protein (p53). The first two parameters make use of a surface clustering method to calculate the protein surface area of highly conserved regions or regions with high nonlocal atomic interaction energy (ANOLEA) score. These parameters help identify important functional regions on the surface of a protein. The last parameter involves the use of a new method for pseudobinding free-energy estimation to specifically probe the importance of residue side-chains to the stability of protein fold. A decision tree was designed to optimally combine these three parameters. The result was compared to the functional data stored in the International Agency for Research on Cancer (IARC) TP53 mutation database. The final prediction achieved a prediction accuracy of 70% and a Matthews correlation coefficient of 0.45. It also showed a high specificity of 91.8%. Mutations in the 85 correctly identified important residues represented 81.7% of the total mutations recorded in the database. In addition, the method was able to correctly assign a probable functional or structural role to the residues. Such information could be critical for the interpretation and prediction of the effect of missense mutations, as it not only provided the fundamental explanation of the observed effect, but also helped design the most appropriate laboratory experiment to verify the prediction results.

A novel population of human melanoma-specific CD8 T cells recognizes Melan-AMART-1 immunodominant nonapeptide but not the corresponding decapeptide.

HLA-A2-restricted cytolytic T cells specific for the immunodominant human tumor Ag Melan-A(MART-1) can kill most HLA-matched melanoma cells, through recognition of two naturally occurring antigenic variants, i.e., Melan-A nonamer AAGIGILTV and decamer EAAGIGILTV peptides. Several previous studies have suggested a high degree of TCR cross-reactivity to the two peptides. In this study, we describe for the first time that some T cell clones are exclusively nonamer specific, because they are not labeled by A2/decamer-tetramers and do not recognize the decamer when presented endogenously. Functional assays with peptides gave misleading results, possibly because decamers were cleaved by exopeptidases. Interestingly, nonapeptide-specific T cell clones were rarely Valpha2.1 positive (only 1 of 19 clones), in contrast to the known strong bias for Valpha2.1-positive TCRs found in decamer-specific clones (59 of 69 clones). Molecular modeling revealed that nonapeptide-specific TCRs formed unfavorable interactions with the decapeptide, whereas decapeptide-specific TCRs productively created a hydrogen bond between CDR1alpha and glutamic acid (E) of the decapeptide. Ex vivo analysis of T cells from melanoma metastases demonstrated that both nonamer and decamer-specific T cells were enriched to substantial frequencies in vivo, and representative clones showed efficient tumor cell recognition and killing. We conclude that the two peptides should be regarded as distinct epitopes when analyzing tumor immunity and developing immunotherapy against melanoma.

EuroPhenome: a repository for high-throughput mouse phenotyping data.

The broad aim of biomedical science in the postgenomic era is to link genomic and phenotype information to allow deeper understanding of the processes leading from genomic changes to altered phenotype and disease. The EuroPhenome project (http://www.EuroPhenome.org) is a comprehensive resource for raw and annotated high-throughput phenotyping data arising from projects such as EUMODIC. EUMODIC is gathering data from the EMPReSSslim pipeline (http://www.empress.har.mrc.ac.uk/) which is performed on inbred mouse strains and knock-out lines arising from the EUCOMM project. The EuroPhenome interface allows the user to access the data via the phenotype or genotype. It also allows the user to access the data in a variety of ways, including graphical display, statistical analysis and access to the raw data via web services. The raw phenotyping data captured in EuroPhenome is annotated by an annotation pipeline which automatically identifies statistically different mutants from the appropriate baseline and assigns ontology terms for that specific test. Mutant phenotypes can be quickly identified using two EuroPhenome tools: PhenoMap, a graphical representation of statistically relevant phenotypes, and mining for a mutant using ontology terms. To assist with data definition and cross-database comparisons, phenotype data is annotated using combinations of terms from biological ontologies.