HTPSELEX--a database of high-throughput SELEX libraries for transcription factor binding sites.

HTPSELEX is a public database providing access to primary and derived data from high-throughput SELEX experiments aimed at characterizing the binding specificity of transcription factors. The resource is primarily intended to serve computational biologists interested in building models of transcription factor binding sites from large sets of binding sequences. The guiding principle is to make available all information that is relevant for this purpose. For each experiment, we try to provide accurate information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, assembled clone sequences (concatemers) and complete sets of in vitro selected protein-binding tags. In addition, we offer in-house derived binding sites models. HTPSELEX also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols. The FTP site contains the trace archives and database flatfiles. The web server offers user-friendly interfaces for viewing individual entries and quality-controlled download of SELEX sequence libraries according to a user-defined sequencing quality threshold. HTPSELEX is available from ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex/ and http://www.isrec.isb-sib.ch/htpselex.

High-throughput SELEX SAGE method for quantitative modeling of transcription-factor binding sites.

The ability to determine the location and relative strength of all transcription-factor binding sites in a genome is important both for a comprehensive understanding of gene regulation and for effective promoter engineering in biotechnological applications. Here we present a bioinformatically driven experimental method to accurately define the DNA-binding sequence specificity of transcription factors. A generalized profile was used as a predictive quantitative model for binding sites, and its parameters were estimated from in vitro-selected ligands using standard hidden Markov model training algorithms. Computer simulations showed that several thousand low- to medium-affinity sequences are required to generate a profile of desired accuracy. To produce data on this scale, we applied high-throughput genomics methods to the biochemical problem addressed here. A method combining systematic evolution of ligands by exponential enrichment (SELEX) and serial analysis of gene expression (SAGE) protocols was coupled to an automated quality-controlled sequence extraction procedure based on Phred quality scores. This allowed the sequencing of a database of more than 10,000 potential DNA ligands for the CTF/NFI transcription factor. The resulting binding-site model defines the sequence specificity of this protein with a high degree of accuracy not achieved earlier and thereby makes it possible to identify previously unknown regulatory sequences in genomic DNA. A covariance analysis of the selected sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism.

Aptamers binding to Bradykinin

Bradykinin is a peptide of the kinin group, involved in a number of receptor-mediated physiological actions, including inflammation and vasodilation, as well as neuromodulation, neuroprotection and promotion of neurogenesis. Bradykinin is the main ligand of the B2 receptor- the main kinin receptor- which is involved in the cardiac and renal protective effects of kinins in diseases. Antibodies have been considered for a long time as promising therapeutic agents in various fields, especially cancer-related ones. Aptamers, on the other hand, have proven to be an excellent alterative, since they have similar properties to those of monoclonal antibodies, such a high-specificity of recognition and high-affinity binding. Plus, they are developed using in vitro selection procedures and can be reproduced by enzymatic reactions. SELEX is a powerful tool for the development of both DNA and RNA aptamers. The main goal of this project was to design a method to select aptamers against bradykinin using capillary electrophoresis alongside the SELEX technique. The selection was done by comparing the aptamers’ (ssDNA-target complex) electrophoretic mobility with that of the ssDNA and the target, which allowed us to define an appropriate collection window that took into consideration the analytes’ detection time, thus enabling the collection of the desired oligonucleotides. After two selection rounds, the collected pool was sequenced, the affinity was measured and the aptamers’ secondary structure was predicted. We concluded that with only two selection cycles, the original DNA library’s bulk affinity grew around 0.4%. The structural characterization of the aptamers, performed with the aid of the Mfold software, revealed that there are many repetitive motifs amongst them, indicating that the selection process was successful. We have obtained 16 sequences of candidate aptamers as bradykinin ligands of similar sequences and secondary structures whose biological activity should be analyzed after synthesis; mainly in regard to their role as bradykinin inhibitors.

Development of a novel aptamer-based multi-sensor device for the detection of osteopontin

Doctoral Dissertation for PhD degree in Chemical and Biological Engineering

Quantitative modelling of transcription factor binding sites

Abstract One of the most important issues in molecular biology is to understand regulatory mechanisms that control gene expression. Gene expression is often regulated by proteins, called transcription factors which bind to short (5 to 20 base pairs),degenerate segments of DNA. Experimental efforts towards understanding the sequence specificity of transcription factors is laborious and expensive, but can be substantially accelerated with the use of computational predictions. This thesis describes the use of algorithms and resources for transcriptionfactor binding site analysis in addressing quantitative modelling, where probabilitic models are built to represent binding properties of a transcription factor and can be used to find new functional binding sites in genomes. Initially, an open-access database(HTPSELEX) was created, holding high quality binding sequences for two eukaryotic families of transcription factors namely CTF/NF1 and LEFT/TCF. The binding sequences were elucidated using a recently described experimental procedure called HTP-SELEX, that allows generation of large number (> 1000) of binding sites using mass sequencing technology. For each HTP-SELEX experiments we also provide accurate primary experimental information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, and assembled clone sequences of binding sequences. The database also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols.The database is available at http://wwwisrec.isb-sib.ch/htpselex/ and and ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex. The Expectation-Maximisation(EM) algorithm is one the frequently used methods to estimate probabilistic models to represent the sequence specificity of transcription factors. We present computer simulations in order to estimate the precision of EM estimated models as a function of data set parameters(like length of initial sequences, number of initial sequences, percentage of nonbinding sequences). We observed a remarkable robustness of the EM algorithm with regard to length of training sequences and the degree of contamination. The HTPSELEX database and the benchmarked results of the EM algorithm formed part of the foundation for the subsequent project, where a statistical framework called hidden Markov model has been developed to represent sequence specificity of the transcription factors CTF/NF1 and LEF1/TCF using the HTP-SELEX experiment data. The hidden Markov model framework is capable of both predicting and classifying CTF/NF1 and LEF1/TCF binding sites. A covariance analysis of the binding sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism. We next tested the LEF1/TCF model by computing binding scores for a set of LEF1/TCF binding sequences for which relative affinities were determined experimentally using non-linear regression. The predicted and experimentally determined binding affinities were in good correlation.

RNA and DNA aptamers as potential tools to prevent cell adhesion in disease

Recent research has shown that receptor-ligand interactions between surfaces of communicating cells are necessary prerequisites for cell proliferation, cell differentiation and immune defense. Cell-adhesion events have also been proposed for pathological conditions such as cancer growth, metastasis, and host-cell invasion by parasites such as Trypanosoma cruzi. RNA and DNA aptamers (aptus = Latin, fit) that have been selected from combinatorial nucleic acid libraries are capable of binding to cell-adhesion receptors leading to a halt in cellular processes induced by outside signals as a consequence of blockage of receptor-ligand interactions. We outline here a novel approach using RNA aptamers that bind to T. cruzi receptors and interrupt host-cell invasion in analogy to existing procedures of blocking selectin adhesion and function in vitro and in vivo.

Design and applications of modified oligonucleotides

Oligonucleotides have a wide range of applications in fields such as biotechnology, molecular biology, diagnosis and therapy. However, the spectrum of uses can be broadened by introducing chemical modifications into their structures. The most prolific field in the search for new oligonucleotide analogs is the antisense strategy, where chemical modifications confer appropriate characteristics such as hybridization, resistance to nucleases, cellular uptake, selectivity and, basically, good pharmacokinetic and pharmacodynamic properties. Combinatorial technology is another research area where oligonucleotides and their analogs are extensively employed. Aptamers, new catalytic ribozymes and deoxyribozymes are RNA or DNA molecules individualized from a randomly synthesized library on the basis of a particular property. They are identified by repeated cycles of selection and amplification, using PCR technologies. Modified nucleotides can be introduced either during the amplification procedure or after selection.

The role of Ken&Barbie in JAK/STAT signalling during development of Drosophila melanogaster

Cell-cell interactions during embryonic development are crucial in the co-ordination of growth, differentiation and maintenance of many different cell types. To achieve this co-ordination each cell must properly translate signals received from neighbouring cells, into spatially and temporally appropriate developmental responses. A surprisingly limited number of signal pathways are responsible for the differentiation of enormous variety of cell types. As a result, pathways are frequently 'reused' during development. Thus, in mammals the JAK/STAT pathway is required during early embryogenesis, mammary gland formation, hematopoiesis and, finally, plays a pivotal role in immune response. In the canonical way, the JAK/STAT pathway is represented by a transmembrane receptor associated with a Janus kinase (JAK), which upon stimulation by an extra-cellular ligand, phosphorylates itself, the receptor and, finally, the signal transducer and activator of transcription (STAT) molecules. Phosphorylated STATs dimerise and translocate to the nucleus where they activate transcription of target genes. The JAK/STAT pathway has been conserved throughout evolution, and all known components are present in the genome of Drosophila melanogaster. Besides hematopoietic and immunity functions, the pathway is also required during development for processes including embryonic segmentation, tracheal morphogenesis, posterior spiracle formation etc. This study describes Drosophila Ken&Barbie (Ken) as a selective regulator of JAK/STAT signalling. ken mutations identified in a screen for modulators of an eye overgrowth phenotype, caused by over-expression of the pathway ligand unpaired, also interact genetically with the pathway receptor domeless (dome) and the transcription factor stat92E. Over-expression of Ken can phenocopy developmental defects known to be caused by the loss of JAK/STAT signalling. These genetic interactions suggest that Ken may function as a negative regulator of the pathway. Ken has C-terminal Zn-finger domain, presumably for DNA binding, and N-terminal BTB/POZ domain, often found in transcriptional repressors. Using EGFP-fused construct expressed in vivo revealed nuclear accumulation of Ken. Therefore, it is proposed that Ken may act as a suppresser of STAT92E target genes. An in vitro assay, termed SELEX, determined that Ken specifically binds to a DNA sequence, with the essential for DNA recognition core overlapping that of STAT92E. This interesting observation suggests that not all STAT92E sites may also allow Ken binding. Strikingly, when effects of ectopic Ken on the expression of putative JAK/STAT pathway target genes were examined, only a subset of the genes tested, namely vvl, trh and kni, were down-regulated by Ken, whereas some others, such as eve and fj, appeared to be unresponsive. Further analysis of vvl, one of the genes susceptible to ectopic Ken, was undertaken. In the developing hindgut, expression of vvl is JAK/STAT pathway dependent, but remains repressed in the posterior spiracles, despite the stimulation of STAT92E by Upd in their primordia. Importantly, ken is also expressed in the developing posterior spiracles. Strikingly, up-regulation of vvl is observed in these tissues in ken mutant embryos. These imply that while ectopic Ken is sufficient to repress the expression of vvl in the hindgut, endogenous Ken is also necessary to prevent its activation in the posterior spiracles. It is therefore conceivable that ectopic vvl expression in the posterior spiracles of the ken mutants may be the result of de-repression of endogenous STAT92E activity. Another consequence of these observations is a fine balance that must exist between STAT92E and Ken activities. Apparently, endogenous level of Ken is sufficient to repress vvl, but not other, as yet unidentified, JAK/STAT pathway targets, whose presumable activation by STAT92E is required for posterior spiracle development as the embryos mutant for dome, the receptor of the pathway, show severe spiracle defects. These defects are also observed in the embryos mis-expressing Ken. Though it is possible that the posterior spiracle phenotype caused by higher levels of Ken results from a JAK/STAT pathway independent activity, it seems to be more likely that Ken acts in a dosage dependent manner, and extra Ken is able to further antagonise JAK/STAT pathway target genes. While STAT92E binding sites required for target gene expression have been poorly characterised, the existence of genome data allows the prediction of candidate STAT92E sites present in target genes promoters to be attempted. When a 6kb region containing the putative regulatory domains flanking the vvl locus are examined, only a single potential STAT92E binding site located 825bp upstream of the translational start can be detected. Strikingly, this site also includes a perfect Ken binding sequence. Such an in silico observation, though consistent with both Ken DNA binding assay in vitro and regulation of STAT92E target genes in vivo, however, requires further analysis. The JAK/STAT pathway is implicated in a variety of processes during embryonic and larval development as well as in imago. In each case, stimulation of the same transcription factor results in different developmental outcomes. While many potential mechanisms have been proposed and demonstrated to explain such pleiotropy, the present study indicates that Ken may represent another mechanism, with which signal transduction pathways are controlled. Ken selectively down-regulates a subset of potential target genes and so modifies the transcriptional profile generated by activated STAT92E - a mechanism, which may be partially responsible for differences in the morphogenetic processes elicited by JAK/STAT signalling during development.

An RNA molecule that specifically inhibits G-protein-coupled receptor kinase 2 in vitro

G-protein-coupled receptors are desensitized by a two-step process. In a first step, G-protein-coupled receptor kinases (GRKs) phosphorylate agonist-activated receptors that subsequently bind to a second class of proteins, the arrestins. GRKs can be classified into three subfamilies, which have been implicated in various diseases. The physiological role(s) of GRKs have been difficult to study as selective inhibitors are not available. We have used SELEX (systematic evolution of ligands by exponential enrichment) to develop RNA aptamers that potently and selectively inhibit GRK2. This process has yielded an aptamer, C13, which bound to GRK2 with a high affinity and inhibited GRK2-catalyzed rhodopsin phosphorylation with an IC50 of 4.1 nM. Phosphorylation of rhodopsin catalyzed by GRK5 was also inhibited, albeit with 20-fold lower potency (IC50 of 79 nM). Furthermore, C13 reveals significant specificity, since almost no inhibitory activity was detectable testing it against a panel of 14 other kinases. The aptamer is two orders of magnitude more potent than the best GRK2 inhibitors described previously and shows high selectivity for the GRK family of protein kinases.

Nuclear dependence of charm production

Using data taken by SELEX during the 1996-1997 fixed target run at Fermilab, we study the production of charmed hadrons on copper and carbon targets with Sigma(-), p, pi(-), and pi(+) beams. Parametrizing the dependence of the inclusive production cross section on the atomic number A as A(alpha), we determine alpha for D(+), D(0), D(s)(+), D(+)(2010), Lambda(+)(c), and their respective anti-particles, as a function of their transverse momentum p(t) and scaled longitudinal momentum x(F). Within our statistics there is no dependence of alpha on x(F) for any charm species for the interval 0.1 < x(F) < 1.0. The average value of alpha for charm production by pion beams is alpha(meson) = 0.850 +/- 0.028. This is somewhat larger than the corresponding average alpha(baryon) = 0.755 +/- 0.016 for charm production by baryon beams (Sigma(-), p).

First observation of the Cabibbo-suppressed decays Xi(+)(c)->Sigma(+)pi(-)pi(+) and Xi(+)(c)->Sigma(-)pi(+)pi(+) and measurement of their branching ratios

We report the first observation of two Cabibbo-suppressed ecay modes, Xi(+)(c) -> Sigma(+)pi(-)pi(+) and Xi c+ -> Sigma(-)pi(+)pi(+). We observe 59 +/- 14 over a background of 87, and 22 +/- 8 over a background of 13 events, respectively, for the signals. The data were accumulated using the SELEX spectrometer during the 1996-1997 fixed target run at Fermilab, chiefly from a 600 Gev/c Sigma(-) beam. The branching ratios of the decays relative to the Cabibbo-favored Xi c+ -> Xi(-)pi(+)pi(+) are measured to be B(Xi(+)(c) -> Sigma(+)pi(-)pi(+))/B(Xi(+)(c) -> Xi(-)pi(+)pi(+)) = 0.48 +/- 0.20, and B(Xi(+)(c) -> Sigma(-)pi(+)pi(+))/B(Xi(+)(c) -> Sigma(-)pi(+)pi(+)) = 0.18 +/- 0.09, respectively. We also report branching ratios for the same decay modes of the Delta(+)(c) relative to Delta(+)(c) -> pK(-)pi(+.) (C) 2008 Elsevier B.V. All rights reserved.

First measurement of π-e→π-eγ pion virtual compton scattering

Pion virtual compton scattering (VCS) via the reaction π-e→π-eγ was observed in the Fermilab E781 SELEX experiment. SELEX used a 600 GeV/c π- beam incident on target atomic electrons, detecting the incident π- and the final state π-, electron and γ. Theoretical predictions based on chiral perturbation theory are incorporated into a Monte Carlo simulation of the experiment and are compared to the data. The number of reconstructed events (=9) and their distribution with respect to the kinematic variables (for the kinematic region studied) are in reasonable accord with the predictions. The corresponding π- VCS experimental cross section is σ=38.8±13 nb, in agreement with the theoretical expectation of σ=34.7 nb.

Investigations on DNA methylation by Dnmt2 and impact of tRNA modifications on TLR7 stimulation

The incorporation of modified nucleotides into ribonucleic acids (RNAs) is important for their structure and proper function. These modifications are inserted by distinct catalytic macromolecules one of them being Dnmt2. It methylates the Cytidine (C) at position 38 in tRNA to 5-methylcytidine (m5C). Dnmt2 has been a paradigm in this respect, because all of its nearest neighbors in evolution are DNA-cytosine C5-methyltransferases and methylate DNA, while its (own) DNA methyltransferase activity is the subject of controversial reports with rates varying between zero and very weak. This work determines whether the biochemical potential for DNA methylation is present in the enzyme. It was discovered that DNA fragments, when presented as covalent RNA:DNA hybrids in the structural context of a tRNA, can be more efficiently methylated than the corresponding natural tRNA substrate. Additional minor deviations from a native tRNA structure that were seen to be tolerated by Dnmt2 were used for a stepwise development of a composite system of guide RNAs that enable the enzyme to perform cytidine methylation on single stranded DNA in vitro. Furthermore, a proof-of-principle is presented for utilizing the S-adenosyl methionine-analog cofactor SeAdoYn with Dnmt2 to search for new possible substrates in a SELEX-like approach.rnIn innate immunity, nucleic acids can function as pathogen associated molecular patterns (PAMPs) recognized by pattern recognition receptors (PRRs). The modification pattern of RNA is the discriminating factor for toll-like receptor 7 (TLR7) to distinguish between self and non-self RNA of invading pathogens. It was found that a 2'-O-methylated guanosine (Gm) at position18, naturally occurring at this position in some tRNAs, antagonizes recognition by TLR7. In the second part of this work it is pointed out, that recognition extends to the next downstream nucleotide and the effectively recognized molecular detail is actually a methylated dinucleotide. The immune silencing effect of the ribose methylation is most pronounced if the dinucleotide motif is composed of purin nucleobases whereas pyrimidines diminish the effect. Similar results were obtained when the Gm modification was transposed into other tRNA domains. Point mutations abolishing base pairings important for a proper tertiary structure had no effect on the immune stimulatory potential of a Gm modified tRNA. Taken together these results suggest a processive type of RNA inspection by TLR7.rn

Progetto e sviluppo di un banco prova per articolazioni protesiche attive ad un grado di libertà basato su programmable automation controller e motion capture

L'elaborato descrive il lavoro svolto in cinque mesi presso il Centro Protesi INAIL di Budrio (BO), che ha portato allo sviluppo di un banco prova per testare articolazioni elettromeccaniche. I dispositivi target, in particolare, sono stati due gomiti mioelettrici: il primo di produzione interna INAIL e il secondo prodotto da Selex ES in collaborazione con il Centro Protesi stesso. Per il controllo del movimento e l'acquisizione dei segnali elettrici si è scelto il PAC CompactRIO (National Instruments), mentre per le acquisizioni cinematiche di angolo e velocità angolare si è sfruttata la stereofotogrammetria.