Expanding molecular modeling and design tools to non-natural sidechains.

Protein-protein interactions encode the wiring diagram of cellular signaling pathways and their deregulations underlie a variety of diseases, such as cancer. Inhibiting protein-protein interactions with peptide derivatives is a promising way to develop new biological and therapeutic tools. Here, we develop a general framework to computationally handle hundreds of non-natural amino acid sidechains and predict the effect of inserting them into peptides or proteins. We first generate all structural files (pdb and mol2), as well as parameters and topologies for standard molecular mechanics software (CHARMM and Gromacs). Accurate predictions of rotamer probabilities are provided using a novel combined knowledge and physics based strategy. Non-natural sidechains are useful to increase peptide ligand binding affinity. Our results obtained on non-natural mutants of a BCL9 peptide targeting beta-catenin show very good correlation between predicted and experimental binding free-energies, indicating that such predictions can be used to design new inhibitors. Data generated in this work, as well as PyMOL and UCSF Chimera plug-ins for user-friendly visualization of non-natural sidechains, are all available at http://www.swisssidechain.ch. Our results enable researchers to rapidly and efficiently work with hundreds of non-natural sidechains.

ARE PPARS NOVEL THERAPEUTIC TARGETS TO IMPROVE MELANOMA TREATMENT ?

Cutaneous melanoma is an aggressive malignant tumor of melanocytes, the pigment- producing cells of the epidermis, with a high incidence in developed countries. Despite some major clinical breakthroughs in the last few years, efficient therapies for metastatic melanoma, which portends a very bad prognosis, are still lacking. Among the potential therapeutic targets that have been attracting at-tention in melanoma are the peroxisome proliferator-activated receptors (PPARs). These members - a, ß and 7 - of the nuclear hormone receptor family, which are ligand-gated transcription factors endowed with a multitude of functions besides metabolism homeostasis, have displayed promising antitumor properties in a wide range of cancer cells, including melanoma. However, our knowledge of PPARs' functions in this skin cancer is far from complete, making the usefulness of any of the a, ß or 7 isotype as a therapeutic target uncertain. In this work, we showed that all three PPAR isotypes are expressed in normal melanocytes, in most melanoma cell lines and in primary and metastatic melanomas, and that PPAR/3 and 7 display transcriptional activity in normal melanocytes and melanoma cells. We also showed that the PPAR7 agonist rosiglitazone had anti-melanoma properties largely independent of PPAR7 expression, which was widely varying across the different cell lines and melanoma biopsies we evaluated and was not correlated with cell line stage. Consistent with the general view of PPAR7 as a tumor suppressor gene, we found that, in human samples, PPAR7 was less expressed in melanoma than in normal skin. Transcriptornic profiling of metastatic melanoma cells in which PPAR7 was pharmacologically modulated revealed an association with epithelial-to-mesenchymal transition, though the functional relevance of this finding remains to be determined. Collectively, our results suggests that PPAR7 activity in melanoma is highly complex and that a straightforward picture of PPAR7's role in this skin cancer is difficult to draw. In this study, we also provided compelling evidence that thioredoxin interacting protein (TXNIP) is, in melanoma, a bona fide PPAR7 target gene, the expression of which is repressed by PPAR7 activation. Although TXNIP is mostly known as an inhibitor of the major antioxidant thioredoxin, it has demonstrated a range of biological functions and is generally considered as a tumor suppressor gene. Consistently, we found that TXNIP expression is associated with growth arrest of melanoma cells in vitro and that forced expression of TXNIP strongly impairs cell proliferation. Interestingly, we also discovered that TXNIP favors melanoma cell migration while it diminishes their adhesion. Finally, we provided several lines of evidence that TXNIP may regulate these processes at the transcriptional level as well as by direct protein-protein interactions in the plasma membrane. Altogether, our findings suggest that the PPAR7 target TXNIP may be a double-edged sword in melanoma, hindering tumor growth but promoting invasion and dissemination. Experiments to evaluate the net biological outcome of TXNIP modulation in vivo are ongoing. -- Le mélanome cutané est une tumeur maligne agressive des mélanocytes, cellules de l'épiderme qui produisent la mélanine. Ce cancer présente un taux d'incidence élevé dans les pays développés et est grevé d'un pronostic très sombre une fois qu'il a disséminé. Malgré les importants progrès réalisés ces dernières années, aucune thérapie lie s'est encore montrée véritablement efficace contre le mélanome métastatique. Parmi les cibles thérapeutiques potentielles, nombre de groupes de recherche se sont penchés sur les peroxisome proliferator-activated receptors (PPARs). Ces récepteurs - a, ß et 7 - font partie de la famille des récepteurs nucléaires aux hormones, des facteurs de transcription activés par des ligands et dotés d'une multitude de fonctions en sus de la régulation du métabolisme. Ces protéines ont démontré des propriétés anti-tumorales prometteuses dans une large gamme de cellules cancéreuses, y compris le mélanome. Cependant, nous connaissons encore très mal les fonctions des PPARs dans ce cancer de la peau, rendant l'utilité thérapeutique de l'un des isotypes a, ß ou 7 incertaine. Dans ce travail, nous avons montré que les trois isotypes sont exprimés dans les mélanocytes normaux, dans la plupart des lignées de mélanome ainsi que dans des mélanomes primaires et métastatiques; nous avons aussi montré que PPAR/3 et 7 sont actifs sur le plan transcriptionnel dans les mélanocytes normaux et les cellules de mélanome. La rosiglitazone, un agoniste de PPAR7, a démontré des propriétés anti-mélanome essentiellement indépendantes de l'expression de PPAR7, qui semble très variable dans les lignées et les biopsies que nous avons évaluées; de plus, l'expression de PPAR7 n'est pas corrélée avec le stade de la lignée. En accord avec la vision communément admise de PPAR7 comme étant un gène suppresseur de tumeur, nous avons observé dans des échantillons humains que PPAR7 est moins exprimé dans les mélanomes que dans la peau normale. Une étude transcrip- tomique de cellules de mélanome métastatique a révélé que la modulation phar-macologique de PPAR7 est associée avec la transition épithélio-mésenchymateuse, même si la pertinence fonctionnelle de cette trouvaille reste à déterminer. Collec-tivement, ces résultats suggèrent que l'activité de PPAR/y dans le mélanome est hautement complexe et qu'une image claire du rôle de PPAR7 dans ce cancer est difficile à dessiner. Dans cette étude, nous avons également fourni de solides preuves que la thiore-doxin interacting protein (TXNIP) est, dans le mélanome, un gène cible bona fide de PPAR7 dont l'expression est réprimée par l'activation de PPAR7. Bien que TXNIP soit surtout connu comme un inhibiteur de la thiorédoxine -un anti-oxydant majeur - cette protéine a démontré une large gamme de fonctions biologiques et est généralement considérée comme un gène suppresseur de tumeur. En accord avec cette conception, nous avons trouvé que l'expression de TXNIP est associée avec l'arrêt de croissance des cellules de mélanome in vitro et que l'expression forcée de TXNIP freine considérablement la prolifération cellulaire. Nous avons aussi découvert que TXNIP favorise la migration des cellules de mélanome alors qu'elle diminue leur adhésion. Enfin, nous avons obtenu plusieurs preuves que TXNIP pourrait réguler ces processus tant au niveau transcriptionnel que par des interactions protéine-protéine au sein de la membrane plasmique. En conclusion, nos résultats suggèrent que la cible de PPAR7 TXNIP pourrait être une épée à double tranchant dans le mélanome, freinant la croissance tumorale mais favorisant l'invasion et la dissémination. Des expériences permettant d'évaluer l'effet biologique net de la modulation de TXNIP in vivo sont en cours.

Phosphorylation modulates FOXC2 transcriptional program by controlling the high-order interaction with DNA.

Phosphorylation of transcription factors is a rapid and reversible process linking cell signaling and control of gene expression, therefore understanding how it controls the transcription factor functions is one of the challenges of functional genomics. We performed such analysis for the forkhead transcription factor FOXC2 mutated in human hereditary disease lymphedemadistichiasis and important for the development of venous and lymphatic valves and lymphatic collecting vessels. We found that FOXC2 is phosphorylated in a cell-cycle dependent manner on eight evolutionary conserved serine/threonine residues, seven of which are clustered within a 70 amino acid domain. Surprisingly, the mutation of phosphorylation sites or a complete deletion of the domain did not affect the transcriptional activity of FOXC2 in a synthetic reporter assay. However, overexpression of the wild type or phosphorylation-deficient mutant resulted in overlapping but distinct gene expression profiles suggesting that binding of FOXC2 to individual sites under physiological conditions is affected by phosphorylation. To gain a direct insight into the role of FOXC2 phosphorylation, we performed comparative genome-wide location analysis (ChIP-chip) of wild type and phosphorylation-deficient FOXC2 in primary lymphatic endothelial cells. The effect of loss of phosphorylation on FOXC2 binding to genomic sites ranged from no effect to nearly complete inhibition of binding, suggesting a mechanism for how FOXC2 transcriptional program can be differentially regulated depending on FOXC2 phosphorylation status. Based on these results, we propose an extension to the enhanceosome model, where a network of genomic context-dependent DNA-protein and protein-protein interactions not only distinguishes a functional site from a nonphysiological site, but also determines whether binding to the functional site can be regulated by phosphorylation. Moreover, our results indicate that FOXC2 may have different roles in quiescent versus proliferating lymphatic endothelial cells in vivo.

Time-scale and other invariants of integrative mechanical behavior in living cells.

In dealing with systems as complex as the cytoskeleton, we need organizing principles or, short of that, an empirical framework into which these systems fit. We report here unexpected invariants of cytoskeletal behavior that comprise such an empirical framework. We measured elastic and frictional moduli of a variety of cell types over a wide range of time scales and using a variety of biological interventions. In all instances elastic stresses dominated at frequencies below 300 Hz, increased only weakly with frequency, and followed a power law; no characteristic time scale was evident. Frictional stresses paralleled the elastic behavior at frequencies below 10 Hz but approached a Newtonian viscous behavior at higher frequencies. Surprisingly, all data could be collapsed onto master curves, the existence of which implies that elastic and frictional stresses share a common underlying mechanism. Taken together, these findings define an unanticipated integrative framework for studying protein interactions within the complex microenvironment of the cell body, and appear to set limits on what can be predicted about integrated mechanical behavior of the matrix based solely on cytoskeletal constituents considered in isolation. Moreover, these observations are consistent with the hypothesis that the cytoskeleton of the living cell behaves as a soft glassy material, wherein cytoskeletal proteins modulate cell mechanical properties mainly by changing an effective temperature of the cytoskeletal matrix. If so, then the effective temperature becomes an easily quantified determinant of the ability of the cytoskeleton to deform, flow, and reorganize.

One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.

Rational design of 4-aryl-1,2,3-triazoles for indoleamine 2,3-dioxygenase 1 inhibition.

Indoleamine 2,3-dioxygenase 1 (IDO1) is an important therapeutic target for the treatment of diseases such as cancer that involve pathological immune escape. Starting from the scaffold of our previously discovered IDO1 inhibitor 4-phenyl-1,2,3-triazole, we used computational structure-based methods to design more potent ligands. This approach yielded highly efficient low molecular weight inhibitors, the most active being of nanomolar potency both in an enzymatic and in a cellular assay, while showing no cellular toxicity and a high selectivity for IDO1 over tryptophan 2,3-dioxygenase (TDO). A quantitative structure-activity relationship based on the electrostatic ligand-protein interactions in the docked binding modes and on the quantum chemically derived charges of the triazole ring demonstrated a good explanatory power for the observed activities.

The Bacillus subtilis bacteriophage SPP1 G39P delivers and activates the G40P DNA helicase upon interacting with the G38P-bound replication origin.

Initiation of Bacillus subtilis bacteriophage SPP1 replication requires the phage-encoded genes 38, 39 and 40 products (G38P, G39P and G40P). G39P, which does not bind DNA, interacts with the replisome organiser, G38P, in the absence of ATP and with the ATP-activated hexameric replication fork helicase, G40P. G38P, which specifically interacts with the phage replication origin (oriL) DNA, does not seem to form a stable complex with G40P in solution. G39P when complexed with G40P-ATP inactivates the single-stranded DNA binding, ATPase and unwinding activities of G40P, and such effects are reversed by increasing amounts of G38P. Unwinding of a forked substrate by G40P-ATP is increased about tenfold by the addition of G38P and G39P to the reaction mixture. The specific protein-protein interactions between oriL-bound G38P and the G39P-G40P-ATPgammaS complex are necessary for helicase delivery to the SPP1 replication origin. Formation of G38P-G39P heterodimers releases G40P-ATPgammaS from the unstable oriL-G38P-G39P-G40P-ATPgammaS intermediate. G40P-ATPgammaS binds to the origin region, the uncomplexed G38P fraction remains bound to oriL, and the G38P-G39P heterodimer is lost from the complex. We demonstrate that G39P is a component of an oligomeric nucleoprotein complex which plays an important role in the initiation of SPP1 replication.

Cognate peptide-receptor ligand mapping by directed phage display

Background: A rapid phage display method for the elucidation of cognate peptide specific ligand for receptors is described. The approach may be readily integrated into the interface of genomic and proteomic studies to identify biologically relevant ligands.Methods: A gene fragment library from influenza coat protein haemagglutinin (HA) gene was constructed by treating HA cDNA with DNAse I to create 50 ¿ 100 bp fragments. These fragments were cloned into plasmid pORFES IV and in-frame inserts were selected. These in-frame fragment inserts were subsequently cloned into a filamentous phage display vector JC-M13-88 for surface display as fusions to a synthetic copy of gene VIII. Two well characterized antibodies, mAb 12CA5 and pAb 07431, directed against distinct known regions of HA were used to pan the library. Results: Two linear epitopes, HA peptide 112 ¿ 126 and 162¿173, recognized by mAb 12CA5 and pAb 07431, respectively, were identified as the cognate epitopes.Conclusion: This approach is a useful alternative to conventional methods such as screening of overlapping synthetic peptide libraries or gene fragment expression libraries when searching for precise peptide protein interactions, and may be applied to functional proteomics.

Remorins form a novel family of coiled coil-forming oligomeric and filamentous proteins associated with apical, vascular and embryonic tissues in plants.

Remorins form a superfamily of plant-specific plasma membrane/lipid-raft-associated proteins of unknown structure and function. Using specific antibodies, we localized tomato remorin 1 to apical tissues, leaf primordia and vascular traces. The deduced remorin protein sequence contains a predicted coiled coil-domain, suggesting its participation in protein-protein interactions. Circular dichroism revealed that recombinant potato remorin contains an alpha-helical region that forms a functional coiled-coil domain. Electron microscopy of purified preparations of four different recombinant remorins, one from potato, two divergent isologs from tomato, and one from Arabidopsis thaliana , demonstrated that the proteins form highly similar filamentous structures. The diameters of the negatively-stained filaments ranged from 4.6-7.4 nm for potato remorin 1, 4.3-6.2 nm for tomato remorin 1, 5.7-7.5 nm for tomato remorin 2, and 5.7-8.0 nm for Arabidopsis Dbp. Highly polymerized remorin 1 was detected in glutaraldehyde-crosslinked tomato plasma membrane preparations and a population of the protein was immunolocalized in tomato root tips to structures associated with discrete regions of the plasma membrane.

Implicit methods for qualitative modeling of gene regulatory networks.

Advancements in high-throughput technologies to measure increasingly complex biological phenomena at the genomic level are rapidly changing the face of biological research from the single-gene single-protein experimental approach to studying the behavior of a gene in the context of the entire genome (and proteome). This shift in research methodologies has resulted in a new field of network biology that deals with modeling cellular behavior in terms of network structures such as signaling pathways and gene regulatory networks. In these networks, different biological entities such as genes, proteins, and metabolites interact with each other, giving rise to a dynamical system. Even though there exists a mature field of dynamical systems theory to model such network structures, some technical challenges are unique to biology such as the inability to measure precise kinetic information on gene-gene or gene-protein interactions and the need to model increasingly large networks comprising thousands of nodes. These challenges have renewed interest in developing new computational techniques for modeling complex biological systems. This chapter presents a modeling framework based on Boolean algebra and finite-state machines that are reminiscent of the approach used for digital circuit synthesis and simulation in the field of very-large-scale integration (VLSI). The proposed formalism enables a common mathematical framework to develop computational techniques for modeling different aspects of the regulatory networks such as steady-state behavior, stochasticity, and gene perturbation experiments.

Specific targeting of pro-death NMDA receptor signals with differing reliance on the NR2B PDZ ligand.

NMDA receptors (NMDARs) mediate ischemic brain damage, for which interactions between the C termini of NR2 subunits and PDZ domain proteins within the NMDAR signaling complex (NSC) are emerging therapeutic targets. However, expression of NMDARs in a non-neuronal context, lacking many NSC components, can still induce cell death. Moreover, it is unclear whether targeting the NSC will impair NMDAR-dependent prosurvival and plasticity signaling. We show that the NMDAR can promote death signaling independently of the NR2 PDZ ligand, when expressed in non-neuronal cells lacking PSD-95 and neuronal nitric oxide synthase (nNOS), key PDZ proteins that mediate neuronal NMDAR excitotoxicity. However, in a non-neuronal context, the NMDAR promotes cell death solely via c-Jun N-terminal protein kinase (JNK), whereas NMDAR-dependent cortical neuronal death is promoted by both JNK and p38. NMDAR-dependent pro-death signaling via p38 relies on neuronal context, although death signaling by JNK, triggered by mitochondrial reactive oxygen species production, does not. NMDAR-dependent p38 activation in neurons is triggered by submembranous Ca(2+), and is disrupted by NOS inhibitors and also a peptide mimicking the NR2B PDZ ligand (TAT-NR2B9c). TAT-NR2B9c reduced excitotoxic neuronal death and p38-mediated ischemic damage, without impairing an NMDAR-dependent plasticity model or prosurvival signaling to CREB or Akt. TAT-NR2B9c did not inhibit JNK activation, and synergized with JNK inhibitors to ameliorate severe excitotoxic neuronal loss in vitro and ischemic cortical damage in vivo. Thus, NMDAR-activated signals comprise pro-death pathways with differing requirements for PDZ protein interactions. These signals are amenable to selective inhibition, while sparing synaptic plasticity and prosurvival signaling.

The role of ubiquitin specific peptidase 15 (USP 15) in human glioblastoma

Glioblastoma (GBM) is the most common and most aggressive malignant primary brain tumour. Despite the aggressiveness of the applied therapy, the prognosis remains poor with a median survival to of about 15 months. It is important to identify new candidate genes that could have clinical application in this disease. Previous gene expression studies from human GBM samples in our laboratory, revealed Ubiquitin Specific Peptidase 15 (USP15) as a gene with low expression, significantly associated with genomic deletions of the chromosomal region encompassing the USP15 locus. USP15 belongs to the ubiquitin-specific protease (USPs) family of which the main role is the reversion of ubiquitination and thereby stabilization of substrates. Previously, USP15 has been suggested to have a tumour suppressor function via its substrates APC and Caspase 3. We established GBM cell lines that stably express USP15 wt or its catalytic mutant. USP15 expression impairs cell growth by inhibiting cell cycle progression. On the other hand USP15 depletion in GBM cell lines induces cell cycle progression and proliferation. In order to identify the molecular pathways in which USP15 is implicated we aimed to identify protein-binding partners in the GBM cell line LN-229 by Mass spectrometry. As a result we identified eight new proteins that interact with USP15. These proteins are involved in important cellular processes like cytokinesis, cell cycle, cellular migration, and apoptosis. Three of these protein interactions were confirmed by co-immunoprecipitation in four GBM cell lines LN-229, LN428, LN18, LN-Z308. One of the binding proteins is HECTD1 E3 ligase of which the murine homologue promotes the APC-Axin interaction to negatively regulate the Wnt pathway. USP15 can de-ubiquitinate HECTD1 in the LN229 cell line while its depletion led to decrease of HECTD1 in GBM cell lines suggesting stabilizing role for USP15. Moreover, HECTD1 stable expression in LN229 inhibits cell cycle, while its depletion induces cell cycle progression. These results suggest that the USP15-HECTD1 interaction might enhance the antiproliferative effect of HECTD1 in GBM cell lines. Using the TOPflash/FOPflash luciferase system we showed that HECTD1 and USP15 overexpression can attenuate WNT pathway activity, and decrease the Axin2 expression. These data indicate that this new protein interaction of USP15 with HECTD1 results in negative regulation of the WNT pathway in GBM cell lines. Further investigation of the regulation of this interaction or of the protein binding network of HECTD1 in GBM may allow the discovery of new therapeutic targets. Finally PTPIP51 and KIF15 are the other two identified protein partners of USP15. These two proteins are involved in cell proliferation and their depletion in LN-229 cell line led to induction of cell cycle progression. USP15 displays a stabilizing role for them. Hence, these results show that the tumour suppressive role of USP15 in GBM cell line via different molecular mechanisms indicating the multidimensional function of USP15. Résumé Le glioblastome (GBM) est la tumeur primaire la plus fréquente et la plus agressive du cervau caractérisée par une survie médiane d'environ à 15 mois. De précédant travaux effectués au sein de notre laboratoire portant sur l'étude de l'expression de gènes pour des échantillons humains de GBM ont montré que le gène Ubiquitin Specific Peptidase 15 (USP1S) était significativement associée à une délétion locales à 25% des cas. Initialement, les substrats protéiques APC et CaspaseS de USP15 ont conduit à considérer cette protéine comme un suppresseur de tumeur. USP15 appartient à la famille protèsse spécifique de l'ubiquitine (USPs) dont le rôle principal est la réversion de l'ubiquitination et la stabilisation de substrats. Par conséquent, nous avons établi des lignées de cellules de glioblastome qui expriment de manière stable USP15 ou bien son mutant catalytique. Ainsi, nous avons ainsi démontré que l'expression de l'USP15 empêche la croissance cellulaire en inhibant la progression du cycle cellulaire. Inversement, la suppression de l'expression du gène USP15 dans les lignées cellulaires de glioblastome induit la progression du cycle cellulaire et la prolifération. Afin d'identifier les voies moléculaires dans lesquelles sont impliquées USP15, nous avons cherché à identifier les partenaires de liaisons protéiques par spectrométrie de masse dans la lignée cellulaire LN-229. Ainsi, huit nouvelles protéines interagissant avec USP15 ont été identifiées dont la ligase E3 HECTD1. L'homologue murin de Hectdl favorise l'interaction APC-Axin en régulant négativement la voie de signalisation de Wnt. USP15 interagit en désubiquitinant HECTD1 dans la lignée cellulaire LN-229 et provoque ainsi l'atténuation de l'activité de cette voie de signalisation. En conclusion, HECTD1, en interagissant avec USP15, joue un rôle de suppresseur de tumeur dans les lignées cellulaire de GBM.

Distinct OGT-Binding Sites Promote HCF-1 Cleavage.

Human HCF-1 (also referred to as HCFC-1) is a transcriptional co-regulator that undergoes a complex maturation process involving extensive O-GlcNAcylation and site-specific proteolysis. HCF-1 proteolysis results in two active, noncovalently associated HCF-1N and HCF-1C subunits that regulate distinct phases of the cell-division cycle. HCF-1 O-GlcNAcylation and site-specific proteolysis are both catalyzed by O-GlcNAc transferase (OGT), which thus displays an unusual dual enzymatic activity. OGT cleaves HCF-1 at six highly conserved 26 amino acid repeat sequences called HCF-1PRO repeats. Here we characterize the substrate requirements for OGT cleavage of HCF-1. We show that the HCF-1PRO-repeat cleavage signal possesses particular OGT-binding properties. The glutamate residue at the cleavage site that is intimately involved in the cleavage reaction specifically inhibits association with OGT and its bound cofactor UDP-GlcNAc. Further, we identify a novel OGT-binding sequence nearby the first HCF-1PRO-repeat cleavage signal that enhances cleavage. These results demonstrate that distinct OGT-binding sites in HCF-1 promote proteolysis, and provide novel insights into the mechanism of this unusual protease activity.

Divergent evolution and purifying selection of the flaA gene sequences in Aeromonas

BACKGROUND: The bacterial flagellum is the most important organelle of motility in bacteria and plays a key role in many bacterial lifestyles, including virulence. The flagellum also provides a paradigm of how hierarchical gene regulation, intricate protein-protein interactions and controlled protein secretion can result in the assembly of a complex multi-protein structure tightly orchestrated in time and space. As if to stress its importance, plants and animals produce receptors specifically dedicated to the recognition of flagella. Aside from motility, the flagellum also moonlights as an adhesion and has been adapted by humans as a tool for peptide display. Flagellar sequence variation constitutes a marker with widespread potential uses for studies of population genetics and phylogeny of bacterial species. RESULTS: We sequenced the complete flagellin gene (flaA) in 18 different species and subspecies of Aeromonas. Sequences ranged in size from 870 (A. allosaccharophila) to 921 nucleotides (A. popoffii). The multiple alignment displayed 924 sites, 66 of which presented alignment gaps. The phylogenetic tree revealed the existence of two groups of species exhibiting different FlaA flagellins (FlaA1 and FlaA2). Maximum likelihood models of codon substitution were used to analyze flaA sequences. Likelihood ratio tests suggested a low variation in selective pressure among lineages, with an omega ratio of less than 1 indicating the presence of purifying selection in almost all cases. Only one site under potential diversifying selection was identified (isoleucine in position 179). However, 17 amino acid positions were inferred as sites that are likely to be under positive selection using the branch-site model. Ancestral reconstruction revealed that these 17 amino acids were among the amino acid changes detected in the ancestral sequence. CONCLUSION: The models applied to our set of sequences allowed us to determine the possible evolutionary pathway followed by the flaA gene in Aeromonas, suggesting that this gene have probably been evolving independently in the two groups of Aeromonas species since the divergence of a distant common ancestor after one or several episodes of positive selection. REVIEWERS: This article was reviewed by Alexey Kondrashov, John Logsdon and Olivier Tenaillon (nominated by Laurence D Hurst).

A Dependency Parsing Approach to Biomedical Text Mining

Biomedical research is currently facing a new type of challenge: an excess of information, both in terms of raw data from experiments and in the number of scientific publications describing their results. Mirroring the focus on data mining techniques to address the issues of structured data, there has recently been great interest in the development and application of text mining techniques to make more effective use of the knowledge contained in biomedical scientific publications, accessible only in the form of natural human language. This thesis describes research done in the broader scope of projects aiming to develop methods, tools and techniques for text mining tasks in general and for the biomedical domain in particular. The work described here involves more specifically the goal of extracting information from statements concerning relations of biomedical entities, such as protein-protein interactions. The approach taken is one using full parsing—syntactic analysis of the entire structure of sentences—and machine learning, aiming to develop reliable methods that can further be generalized to apply also to other domains. The five papers at the core of this thesis describe research on a number of distinct but related topics in text mining. In the first of these studies, we assessed the applicability of two popular general English parsers to biomedical text mining and, finding their performance limited, identified several specific challenges to accurate parsing of domain text. In a follow-up study focusing on parsing issues related to specialized domain terminology, we evaluated three lexical adaptation methods. We found that the accurate resolution of unknown words can considerably improve parsing performance and introduced a domain-adapted parser that reduced the error rate of theoriginal by 10% while also roughly halving parsing time. To establish the relative merits of parsers that differ in the applied formalisms and the representation given to their syntactic analyses, we have also developed evaluation methodology, considering different approaches to establishing comparable dependency-based evaluation results. We introduced a methodology for creating highly accurate conversions between different parse representations, demonstrating the feasibility of unification of idiverse syntactic schemes under a shared, application-oriented representation. In addition to allowing formalism-neutral evaluation, we argue that such unification can also increase the value of parsers for domain text mining. As a further step in this direction, we analysed the characteristics of publicly available biomedical corpora annotated for protein-protein interactions and created tools for converting them into a shared form, thus contributing also to the unification of text mining resources. The introduced unified corpora allowed us to perform a task-oriented comparative evaluation of biomedical text mining corpora. This evaluation established clear limits on the comparability of results for text mining methods evaluated on different resources, prompting further efforts toward standardization. To support this and other research, we have also designed and annotated BioInfer, the first domain corpus of its size combining annotation of syntax and biomedical entities with a detailed annotation of their relationships. The corpus represents a major design and development effort of the research group, with manual annotation that identifies over 6000 entities, 2500 relationships and 28,000 syntactic dependencies in 1100 sentences. In addition to combining these key annotations for a single set of sentences, BioInfer was also the first domain resource to introduce a representation of entity relations that is supported by ontologies and able to capture complex, structured relationships. Part I of this thesis presents a summary of this research in the broader context of a text mining system, and Part II contains reprints of the five included publications.