Probing the Unfolded Configurations of a β-Hairpin Using Sketch-Map

This work examines the conformational ensemble involved in β-hairpin folding by means of advanced molecular dynamics simulations and dimensionality reduction. A fully atomistic description of the protein and the surrounding solvent molecules is used, and this complex energy landscape is sampled by means of parallel tempering metadynamics simulations. The ensemble of configurations explored is analyzed using the recently proposed sketch-map algorithm. Further simulations allow us to probe how mutations affect the structures adopted by this protein. We find that many of the configurations adopted by a mutant are the same as those adopted by the wild-type protein. Furthermore, certain mutations destabilize secondary-structure-containing configurations by preventing the formation of hydrogen bonds or by promoting the formation of new intramolecular contacts. Our analysis demonstrates that machine-learning techniques can be used to study the energy landscapes of complex molecules and that the visualizations that are generated in this way provide a natural basis for examining how the stabilities of particular configurations of the molecule are affected by factors such as temperature or structural mutations.

Estratégias de alteração da funcionalidade de proteínas de soja

Este trabalho teve como principal objetivo estudar e modificar as propriedades funcionais das proteínas de soja de forma a otimizar e diversificar a sua aplicação industrial. Para tal, foram propostas e estudadas quatro estratégias: i) extração do isolado de proteínas de soja (IPS) a partir de diferentes matérias-primas, ii) adição de galactomananas (GM) com graus de ramificação e massas moleculares diferentes, iii) hidrólise enzimática controlada das proteínas de soja, iv) processamento por alta pressão hidrostática. O estudo e a interpretação da influência destas estratégias sobre as propriedades funcionais das proteínas de soja, nomeadamente, na capacidade gelificante e emulsionante, foram realizados recorrendo fundamentalmente a ensaios reológicos dinâmicos a baixas deformação, espectroscopia de infravermelho, electroforeses, calorimetria diferencial de varrimento e ensaios de microscopia confocal de varrimento laser. O estudo da extração e caracterização dos isolados de proteínas de soja obtidos a partir de diferentes matérias-primas permitiu concluir que as caraterísticas físico-químicas dos isolados são dependentes da origem da matéria-prima de extração e da severidade dos tratamentos industriais prévios à extração do isolado. Contudo, as propriedades viscoelásticas dos géis obtidos por aquecimento controlado não foram significativamente distintas embora tenha sido possível relacionar o grau de agregação com a diminuição da temperatura de gelificação e com o aumento inicial dos módulos viscoelásticos. As alterações sofridas pelos isolados de origem comercial mostraram ser irreversíveis resultando em géis menos rígidos e com maior caráter viscoso. A adição de galactomanana alterou significativamente o mecanismo de gelificação induzido termicamente das proteínas de soja, bem como as propriedades viscoelásticas dos géis e a microestrutura dos géis, demonstrando-se a ocorrência de separação de fases, em virtude da incompatibilidade termodinâmica entre os biopolímeros, resultando em géis mais rígidos e no decréscimo da temperatura de gelificação. A extensão destas alterações foi dependente da massa molecular, grau de ramificação e da razão IPS/GM. O efeito da hidrólise enzimática por ação da bromelina, nas propriedades gelificantes e emulsionantes das proteínas de soja, mostrou ser dependente do grau de hidrólise (GH). Valores de GH inferiores a 15 % melhoraram as propriedades gelificantes das proteínas de soja. Por outro lado, o aumento do GH teve um efeito negativo nas propriedades emulsionantes, o qual foi atenuado por adição da goma de alfarroba, com efeito positivo na gelificação das proteínas de soja. A concentração crítica limite de compatibilidade entre os hidrolisados de proteína de soja e a goma de alfarroba aumentou com o decréscimo do GH e da massa molecular do polissacacrídeo. O efeito da AP sobre as propriedades físico-químicas e funcionais dos IPS foi influenciado pela origem do isolado e pelas condições de tratamento. O processamento até 100 MPa desencadeou um aumento da atividade emulsionante e considerável melhoria da capacidade gelificante. Contudo, valores de pressão superiores promoveram a desnaturação das proteínas constituintes dos isolados, resultando no decréscimo da temperatura de gelificação e numa re-associação das subunidades proteicas, diminuindo a elasticidade dos géis finais. Os resultados sugeriram que as alterações nas proteínas de soja promovidas durante o tratamento por AP constituem um fator limitante para o desdobramento e re-associação durante o aquecimento térmico, necessários para a formação e fortalecimento de gel formado. O processamento por AP influenciou a estrutura secundária e a microestrutura das amostras. A presença de GA teve um papel baroprotetor. Assim, com este trabalho demonstrou-se que com as estratégias seguidas para manipulação das propriedades funcionais de proteínas de soja, nomeadamente através da adição de um polissacarídeo com propriedades estruturais controladas, da adequada combinação da adição de um polissacarídeo neutro com a hidrólise controlada das proteínas ou com tratamento por alta pressão, é possível a criação de novas funcionalidades, com utilidade no desenvolvimento de novas formulações alimentares, permitindo expandir a aplicação destas proteínas vegetais.

Metal ions modulate the folding and stability of the tumor suppressor protein S100A2

Febs Journal (2009)276:1776-1786

Dimensionamento de Madres Enformadas a Frio Travadas por Painéis do Tipo Sandwich

O presente relatório de estágio foi elaborado no âmbito da Unidade Curricular de DIPRE, Dissertação/Projeto/Estágio, do 2.º ano de Mestrado em Engenharia Civil do Instituto Superior de Engenharia do Porto, do ramo de estruturas, tendo como principal foco, a análise e dimensionamento de madres de aço enformado a frio e sua utilização com painéis do tipo sandwich em coberturas. Seções enformadas a frio são cada vez mais utilizadas em construções modernas, especificamente como estrutura secundária em coberturas, onde geralmente são fixas a painéis através de ligações aparafusadas. A presença de painéis e fixações através de parafusos permitem a estabilização lateral e torsional de madres aumentando desta forma a capacidade resistente, mas por serem elementos estruturais de espessura reduzida, os enformados a frio são suscetíveis a fenómenos de instabilidade associados. Desta forma, a norma EN 1993-1-3 [4] permite a análise e dimensionamento deste tipo de elementos através das disposições regulamentares preconizadas nas partes 1-1 [3] e 1-5 [5] da mesma norma. Num primeiro estudo, o presente trabalho tem como objetivo o dimensionamento e verificação de segurança de elementos enformados a frio com base na seção efetiva determinada com o auxílio das normas EN 1993-1-1 (regras gerais) e EN 1993-1-5 (regras para elementos estruturais constituídos por placas). Numa segunda fase, este trabalho pretende apresentar um estudo do comportamento de interação entre os sistemas madres-painéis. Para tal, são quantificadas as rigidezes das conexões dos sistemas e dos painéis para se realizarem a análise relativamente à restrição lateral e restrição torsional de madres. Neste contexto, concluiu-se que os painéis, quando fixos de forma adequada às madres, contribuem para a estabilidade.

A minimal version of a Protein Tyrosine Phosphatase

Phosphatase and tensin homologue (PTEN) protein belongs to the family of protein tyrosine phos-phatase. Mutations on the phosphatase and tensin homologue (PTEN) protein are highly observed in diverse types of human tumors, being mostly identified on the phosphatase domain of the protein. Although PTEN is a modular protein composed by a phosphatase domain and a C2 domain for mem-brane anchoring, this work aimed at developing a minimal version of PTEN´s phosphatase domain. The minimal version (Small Domain) comprises a 28 residue peptide, with the PTEN 8-mer catalytic peptide accommodated between a α-helix and β-turn as observed in PTEN native structure. Firstly, a de novo prediction of the Small Domain´s secondary structure was carried out by molecular modeling tools. The stability of the predicted structures were then evaluated by Molecular Dynamics. Automated molecular docking of PTEN natural substrate PIP3, its analogue (Inositol) and a PTEN inhibitor (L-tar-tare) were performed with the modeled structure, and PTEN used as a positive control. The gene en-coding for Small Domain was designed and cloned into an expression vector at N-terminal of Green Fluorescence Protein (GFP) encoding gene. The fusion protein was then expressed in Escherichia coli cells. Different expression conditions have been explored for the production of the fusion protein to minimize the formation of inclusion bodies.

Extracting Ramachandran torsional angle distributions from 2D NMR data using Tikhonov regularization /

Solid state nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for studying structural and dynamical properties of disordered and partially ordered materials, such as glasses, polymers, liquid crystals, and biological materials. In particular, twodimensional( 2D) NMR methods such as ^^C-^^C correlation spectroscopy under the magicangle- spinning (MAS) conditions have been used to measure structural constraints on the secondary structure of proteins and polypeptides. Amyloid fibrils implicated in a broad class of diseases such as Alzheimer's are known to contain a particular repeating structural motif, called a /5-sheet. However, the details of such structures are poorly understood, primarily because the structural constraints extracted from the 2D NMR data in the form of the so-called Ramachandran (backbone torsion) angle distributions, g{^,'4)), are strongly model-dependent. Inverse theory methods are used to extract Ramachandran angle distributions from a set of 2D MAS and constant-time double-quantum-filtered dipolar recoupling (CTDQFD) data. This is a vastly underdetermined problem, and the stability of the inverse mapping is problematic. Tikhonov regularization is a well-known method of improving the stability of the inverse; in this work it is extended to use a new regularization functional based on the Laplacian rather than on the norm of the function itself. In this way, one makes use of the inherently two-dimensional nature of the underlying Ramachandran maps. In addition, a modification of the existing numerical procedure is performed, as appropriate for an underdetermined inverse problem. Stability of the algorithm with respect to the signal-to-noise (S/N) ratio is examined using a simulated data set. The results show excellent convergence to the true angle distribution function g{(j),ii) for the S/N ratio above 100.

Locating and sequencing of the E3 and neighbouring regions in bovine advenovirus type 2

Adenoviruses are nonenveloped icosahedral shaped particles. The double stranded DNA viral genome is divided into 5 major early transcription units, designated E1 A, E1 B, and E2 to E4, which are expressed in a regulated manner soon after infection. The gene products of the early region 3 (E3), shown to be nonessential for viral replication in vitro, are believed to be involved in counteracting host immunosurveillance. In order to sequence the E3 region of Bovine adenovirus type 2 (BAV2) it was necessary to determine the restriction map for the plasmid pEA48. A physical restriction endonuclease map for BamHl, Clal, Eco RI, Hindlll, Kpnl, Pstt, Sail, and Xbal was constructed. The DNA insert in pEA48 was determined to be viral in origin using Southern hybridization. A human adenovirus type 5 recombinant plasmid, containing partial DNA fragments of the two transcription units L4 and L5 that lie just outside the E3, was used to localize this region. The recombinant plasmid pEA was subcloned to facilitate sequencing. The DNA sequences between 74.8 and 90.5 map units containing the E3, the hexon associated protein (pVIII), and the fibre gene were determined. Homology comparison revealed that the genes for the hexon associated pV11I and the fibre protein are conserved. The last 70 amino acids of the BAV2 pV11I were the most conserved, showing a similarity of 87 percent with Ad2 pV1I1. A comparison between the predicted amino acid sequences of BAV2 and Ad40, Ad41 , Ad2 and AdS, revealed that they have an identical secondary structure consisting of a tail, a shaft and a knob. The shaft is composed of 22, 15 amino acid motifs, with periodic glycines and hydrophobic residues. The E3 region was found to consist of about 2.3 Kbp and to encode four proteins that were greater than 60 amino acids. However, these four open reading frames did not show significant homology to any other known adenovirus DNA or protein sequence.

Représentation et recherche de motifs cycliques et structuraux d’ARN connus dans les structures secondaires

L'acide désoxyribonucléique (ADN) et l'acide ribonucléique (ARN) sont des polymères de nucléotides essentiels à la cellule. À l'inverse de l'ADN qui sert principalement à stocker l'information génétique, les ARN sont impliqués dans plusieurs processus métaboliques. Par exemple, ils transmettent l’information génétique codée dans l’ADN. Ils sont essentiels pour la maturation des autres ARN, la régulation de l’expression génétique, la prévention de la dégradation des chromosomes et le ciblage des protéines dans la cellule. La polyvalence fonctionnelle de l'ARN résulte de sa plus grande diversité structurale. Notre laboratoire a développé MC-Fold, un algorithme pour prédire la structure des ARN qu'on représente avec des graphes d'interactions inter-nucléotidiques. Les sommets de ces graphes représentent les nucléotides et les arêtes leurs interactions. Notre laboratoire a aussi observé qu'un petit ensemble de cycles d'interactions à lui seul définit la structure de n'importe quel motif d'ARN. La formation de ces cycles dépend de la séquence de nucléotides et MC-Fold détermine les cycles les plus probables étant donnée cette séquence. Mon projet de maîtrise a été, dans un premier temps, de définir une base de données des motifs structuraux et fonctionnels d'ARN, bdMotifs, en terme de ces cycles. Par la suite, j’ai implanté un algorithme, MC-Motifs, qui recherche ces motifs dans des graphes d'interactions et, entre autres, ceux générés par MC-Fold. Finalement, j’ai validé mon algorithme sur des ARN dont la structure est connue, tels que les ARN ribosomaux (ARNr) 5S, 16S et 23S, et l'ARN utilisé pour prédire la structure des riborégulateurs. Le mémoire est divisé en cinq chapitres. Le premier chapitre présente la structure chimique, les fonctions cellulaires de l'ARN et le repliement structural du polymère. Dans le deuxième chapitre, je décris la base de données bdMotifs. Dans le troisième chapitre, l’algorithme de recherche MC-Motifs est introduit. Le quatrième chapitre présente les résultats de la validation et des prédictions. Finalement, le dernier chapitre porte sur la discussion des résultats suivis d’une conclusion sur le travail.

Étude d'un variant de la toxine STb produite par Escherichia coli

Les E. coli entérotoxinogènes (ETEC) sont souvent la cause de diarrhée post-sevrage chez le porc. Deux types d’entérotoxines sont retrouvées chez les ETEC, soit les thermolabiles, comme la toxine LT, et les thermostables, comme EAST-1, STa et STb. Cette dernière est composée de 48 acides aminés et est impliquée dans la pathologie causée par les ETEC. Pour la première fois un variant de la toxine STb fut découvert dans une étude. Nous avons alors émis l’hypothèse qu’il y a présence de variants dans la population de souches ETEC du Québec. Dans les 100 souches STb+ analysées, 23 possédaient le gène de la toxine avec une variation dans la séquence génétique : l’asparagine était présente en position 12 remplaçant ainsi l’histidine. Une corrélation entre la présence du variant et la présence de facteurs de virulence retrouvés dans ces 100 souches ETEC étudiées a été effectuée. Ce variant semble fortement associé à la toxine STa puisque toutes les souches variantes ont hybridé avec le gène codant pour cette dernière. Étant donné sa présence répandue dans la population de souches ETEC du Québec, nous avons de plus émis l’hypothèse que ce variant a des caractéristiques biologiques altérées par rapport à la toxine sauvage. L’analyse par dichroïsme circulaire a montré que le variant et la toxine sauvage ont une structure secondaire ainsi qu’une stabilité similaires. Par la suite, l’attachement au récepteur de la toxine, le sulfatide, a été étudié par résonnance plasmonique de surface (biacore). Le variant a une affinité au sulfatide légèrement réduite comparativement à la toxine sauvage. Puisque l’internalisation de la toxine fut observée dans une étude précédente et qu’elle semble liée à la toxicité, nous avons comparé l’internalisation du variant et de la toxine sauvage à l’intérieur des cellules IPEC-J2. L’internalisation du variant dans les cellules est légèrement supérieure à l’internalisation de la toxine sauvage. Ces résultats suggèrent que le variant est biochimiquement et structurellement comparable à la toxine sauvage.

Synthèse et caractérisation physicochimique de peptides de polyglutamines

Neuf maladies neurodégénératives sont le produit de l’expression de gènes mutés, dans lesquels le codon CAG est répété au-delà d’un seuil pathologique. Ceci produit des protéines mutantes dans lesquelles sont insérés des segments de polyglutamines (polyGln), qui perdent leur activité et acquièrent une nouvelle fonction, ce qui est toxique pour le neurone. Ces altérations sont attribuables aux propriétés particulières de la polyGln. En effet, ces dernières possèdent la capacité de s’assembler pour former des corps d’inclusion intracellulaires. Cette propension à l’agrégation de la polyGln rend difficile l’étude de ces pathologies. C’est ainsi que l’utilisation de peptides peut s’avérer une approche avantageuse. Toutefois, la synthèse de polyGln est associée à de nombreuses délétions et nécessite l’ajout de groupements chargés afin de permettre leur purification. Cependant, ce prérequis donne lieu à des interactions électrostatiques qui biaisent la structure et la cinétique d’agrégation de ces peptides, en plus d’interférer avec l’évaluation d’éventuels agents thérapeutiques. L’objectif du projet est de développer un système permettant l’étude de la polyGln en s’affranchissant des effets de charges. Pour ce faire, deux approches ont été explorées, la première utilise la polyGln non chargée et la seconde utilise une structure polyGln-morpholine ayant des charges labiles en fonction du pH. Ces peptides ont été produits en utilisant une approche linéaire de synthèse peptidique sur support solide avec protection maximale des chaînes latérales. La purification a été effectuée par chromatographie de haute performance en phase inverse en milieu acide. Ces stratégies ont permis de produire des peptides de polyGln de grande pureté avec des rendements acceptables. Une procédure de solubilisation des peptides alliant sonication et lyophilisation a été développée afin d’étudier chacun de ces peptides à l’aide de diverses techniques physicochimiques, telles que la diffusion de la lumière, la spectroscopie de résonance magnétique nucléaire, Raman et UV-visible, le dichroïsme circulaire et la microscopie optique polarisée. La polyGln non chargée solubilisée dans le trifluoroéthanol-eau a montré que la taille des particules et la vitesse d’agrégation sont proportionnelles à la fraction volumique en eau. De plus, la structure secondaire en solution est à prédominance alpha et semble être peu sensible à la fraction d’eau jusqu’à un certain seuil (25%) après lequel la structure aléatoire prédomine. L’analyse des agrégats à l’état solide montre des structures hélicoïdales > aléatoires et ont les caractéristiques des fibrilles amyloïdes. Le peptide de polyGln-morpholines a un pKa de 7,3 en milieu aqueux. Il demeure en solution lorsque le pH < pKa et à faible force ionique, alors qu’il s’autoassemble lorsque ces conditions ne sont pas respectées. Ceci suggère que la répulsion électrostatique est responsable de la stabilisation du peptide en solution. La dimension fractale nous indique que le peptide forme des agrégats compacts dont les constituants ont une taille de 2,5 nm, compatibles avec une conformation aléatoire compacte, en coude bêta ou hélicoïdale. Ceci est en accord avec l’étude structurale des peptides en solution qui a montré des espèces aléatoires > bêta > alpha. De plus, en RMN, l’élargissement des signaux du 1Hγ en cours d’agrégation suggère une interaction via les chaînes latérales. Les analyses en phase solide ont plutôt montré une prédominance de structures bêta et alpha. L’inhibition de l’agrégation à pH 8 varie selon rouge de Congo > tréhalose, alors que le peptide liant la polyGln 1 et la thioflavine T ne semble pas avoir d’effet. Ces approches ont donc permis pour la première fois de s’affranchir des effets de charges auparavant inhérents à l’étude de la polyGln en solution et par conséquent d’obtenir des informations inédites quant à la solubilité, la structure et la cinétique d’agrégation. Enfin, le dispositif à charges labiles permet d’évaluer l’efficacité d’éventuels agents thérapeutiques à pH quasi physiologique.

Riboswitches : le cas des atténuateurs de la transcription du type terminateur/antiterminateur chez les bactéries

Il est essentiel pour chaque organisme d’avoir la possibilité de réguler ses fonctions afin de permettre sa survie et d’améliorer sa capacité de se reproduire en divers habitats. Avec l’information disponible, il semble que les organismes consacrent une partie assez importante de leur matériel génétique à des fonctions de régulation. On peut envisager que certains mécanismes de régulation ont persisté dans le temps parce qu’ils remplissent bien leurs rôles. Les premières études sur les procaryotes ont indiqué qu’il y avait peu de mécanismes de régulation exerçant le contrôle des gènes, mais il a été démontré par la suite qu’une variété de ces mécanismes est utilisée pour la régulation de gènes et d’opérons. En particulier, les opérons bactériens impliqués dans la biosynthèse des acides aminés, l’ARNt synthétase, la dégradation des acides aminés, les protéines ribosomales et l’ARN ribosomal font l’objet d’un contrôle par l’atténuation de la transcription. Ce mécanisme d’atténuation de la transcription diffère d’autres mécanismes pour la génération de deux structures différentes de l’ARNm, où l’une de ces structures réprime le gène en aval, et l’autre permet de continuer la transcription/traduction. Dans le cadre de cette recherche, nous nous sommes intéressé au mécanisme d’atténuation de la transcription chez les procaryotes où aucune molécule ne semble intervenir comme facteur de régulation, en me concentrant sur la régulation des opérons bactériens. Le but principal de ce travail est de présenter une nouvelle méthode de recherche des riborégulateurs qui combine la recherche traditionnelle des riborégulateurs avec la recherche structurale. En incorporant l’étude du repliement de l’ARNm, nous pouvons mieux identifier les atténuateurs répondant à ce type de mécanisme d’atténuation. Ce mémoire est divisé en quatre chapitres. Le premier chapitre présente une revue de la littérature sur l’ARN et un survol sur les mécanismes de régulation de l’expression génétique chez les procaryotes. Les chapitres 2 et 3 sont consacrés à la méthodologie utilisée dans cette recherche et à l’implémentation du logiciel TA-Search. Enfin, le chapitre 4 expose les conclusions et les applications potentielles de la méthode.

Développement de peptidomimétiques antagonistes du récepteur de l’interleukine-1β

Dans ce mémoire, je présente mes études sur la synthèse, la caractérisation et l’évaluation biologique de différentes séries d’analogues du D-heptapeptide appelé 101.10, un modulateur négatif allostérique du récepteur de l’interleukine-1β (IL-1β). Sachant que les peptides ont généralement de faibles propriétés pharmacologiques, le but de ce projet portait sur l’examen des structures nécessaires à la bioactivité, la conformation tridimensionnelle de ces derniers afin d’améliorer la droguabilité du peptide parent. Les stratégies d’optimisation du 101.10 utilisées furent : la coupure N- et C-terminale; la substitution par la proline, α-amino-γ-lactame (Agl), β-amino-γ-lactame (Bgl) et α-amino-β-hydroxy-γ-lactame (Hgl); et la rigidification du squelette à l’aide d’un bicycle, l’indolozidin-2-one (I2aa). Afin de clarifier certaines relations de structure-activité, quelques modifications furent apportées au peptide, incluant l’échange de la thréonine pour la valine, la permutation de la stéréochimie de certains résidus clés ainsi que le remplacement de certaines chaînes latérales par un méthyle. Pour pallier aux difficultés de reproductibilité des résultats avec des échantillons provenant de différentes sources, des études sur l’identité du contre-anion et la pureté du peptide furent conduites. Afin d’évaluer l’effet des modifications sur la conformation aqueuse et l’activité biologique du peptide, des analyses de dichroïsme circulaire et des tests in vitro mesurant l’inhibition de certains effets de l’IL-1β furent effectués. Ces essais cellulaires comportaient l’inhibition de la prolifération de cellules immunes et de l’activation des voies de signalisation inflammatoires du facteur nucléaire κB (NF-κB) et de la protéine kinase activée par mitogène (MAPK), toutes deux stimulées par l’IL-1β. La compilation de ces données a permis de déceler certaines tendances entre la structure, la conformation et l’activité anti-IL-1β des peptidomimétiques.

La RNase P mitochondriale chez Neurospora crassa

Résumé La Ribonucléase P (RNase P) est une enzyme principalement reconnue pour sa participation à la maturation en 5’des ARN de transfert (ARNt). Cependant, d’autres substrats sont reconnus par l’enzyme. En général, la RNase P est composée d’une sous-unité ARN (le P-ARN, codé par le gène rnpB) qui porte le centre actif de l’enzyme et d’une ou de plusieurs sous-unités protéiques (la P-protéine). Les P-ARN chez toutes les bactéries, la majorité des archéobactéries et dans le génome nucléaire de la plupart des eucaryotes, possèdent généralement une structure secondaire très conservée qui inclut le noyau (P1-P4); l’hélice P4 constitue le site catalytique de l’enzyme et l’hélice P1 apparie les extrémités du P-ARN en stabilisant sa structure globale. Les P-ARN mitochondriaux sont souvent moins conservés et difficiles à découvrir. Dans certains cas, les seules régions de structure primaire qui restent conservées sont celles qui définissent le P4 et le P1. Pour la détection des gènes rnpB, un outil de recherche bioinformatique, basé sur la séquence et le profil de structure secondaire, a été développé dans le laboratoire. Cet outil permet le dépistage de toutes les séquences eucaryotes (nucléaires et mitochondriales) du gène avec une très grande confiance (basée sur une valeur statistique, E-value). Chez les champignons, plusieurs ascomycètes encodent un gène rnpB dans leur génome mitochondrial y compris tous les membres du genre d’Aspergillus. Cependant, chez les espèces voisines, Neurospora crassa, Podospora anserina et Sordaria macrospora, une version mitochondriale de ce gène n’existe pas. Au lieu de cela, elles contiennent deux copies nucléaires du gène, légèrement différentes en taille et en contenu nucléotidique. Mon projet a été établi dans le but d’éclaircir l’évolution de la RNase P mitochondriale (mtRNase P) chez ces trois espèces voisines d’Aspergillus. En ce qui concerne les résultats, des modèles de structures secondaires pour les transcrits de ces gènes ont été construits en se basant sur la structure consensus universelle de la sous-unité ARN de la RNase P. Pour les trois espèces, par la comparaison de ces modèles, nous avons établi que les deux copies nucléaires du gène rnpB sont assez distinctes en séquence et en structure pour pouvoir y penser à une spécialisation de fonction de la RNase P. Chez N. crassa, les deux P-ARN sont modifiés probablement par une coiffe et les extrémités 5’, 3’ sont conformes à nos modèles, ayant un P1 allongé. Encore chez N. crassa, nous avons constaté que les deux copies sont transcrites au même niveau dans le cytoplasme et que la plus petite et la plus stable d’entre elles (Nc1) se retrouve dans l’extrait matriciel mitochondrial. Lors du suivi du P-ARN dans diverses sous-fractions provenant de la matrice mitochondriale soluble, Nc1 est associée avec l’activité de la RNase P. La caractérisation du complexe protéique, isolé à partir de la fraction active sur un gel non dénaturant, révèle qu’il contient au moins 87 protéines, 73 d’entre elles ayant déjà une localisation mitochondriale connue. Comme chez la levure, les protéines de ce complexe sont impliquées dans plusieurs fonctions cellulaires comme le processing de l’ADN/ARN, le métabolisme, dans la traduction et d’autres (par exemple : la protéolyse et le repliement des protéines, ainsi que la maintenance du génome mitochondrial). Pour trois protéines, leur fonction est non déterminée.

Une signature du polymorphisme structural d’acides ribonucléiques non-codants permettant de comparer leurs niveaux d’activités biochimiques

Des évidences expérimentales récentes indiquent que les ARN changent de structures au fil du temps, parfois très rapidement, et que ces changements sont nécessaires à leurs activités biochimiques. La structure de ces ARN est donc dynamique. Ces mêmes évidences notent également que les structures clés impliquées sont prédites par le logiciel de prédiction de structure secondaire MC-Fold. En comparant les prédictions de structures du logiciel MC-Fold, nous avons constaté un lien clair entre les structures presque optimales (en termes de stabilité prédites par ce logiciel) et les variations d’activités biochimiques conséquentes à des changements ponctuels dans la séquence. Nous avons comparé les séquences d’ARN du point de vue de leurs structures dynamiques afin d’investiguer la similarité de leurs fonctions biologiques. Ceci a nécessité une accélération notable du logiciel MC-Fold. L’approche algorithmique est décrite au chapitre 1. Au chapitre 2 nous classons les impacts de légères variations de séquences des microARN sur la fonction naturelle de ceux-ci. Au chapitre 3 nous identifions des fenêtres dans de longs ARN dont les structures dynamiques occupent possiblement des rôles dans les désordres du spectre autistique et dans la polarisation des œufs de certains batraciens (Xenopus spp.).

Identification de caractéristiques communes et rares dans les ARN structurés dans la base de données Rfam

Les ARN non codants (ARNnc) sont des transcrits d'ARN qui ne sont pas traduits en protéines et qui pourtant ont des fonctions clés et variées dans la cellule telles que la régulation des gènes, la transcription et la traduction. Parmi les nombreuses catégories d'ARNnc qui ont été découvertes, on trouve des ARN bien connus tels que les ARN ribosomiques (ARNr), les ARN de transfert (ARNt), les snoARN et les microARN (miARN). Les fonctions des ARNnc sont étroitement liées à leurs structures d’où l’importance de développer des outils de prédiction de structure et des méthodes de recherche de nouveaux ARNnc. Les progrès technologiques ont mis à la disposition des chercheurs des informations abondantes sur les séquences d'ARN. Ces informations sont accessibles dans des bases de données telles que Rfam, qui fournit des alignements et des informations structurelles sur de nombreuses familles d'ARNnc. Dans ce travail, nous avons récupéré toutes les séquences des structures secondaires annotées dans Rfam, telles que les boucles en épingle à cheveux, les boucles internes, les renflements « bulge », etc. dans toutes les familles d'ARNnc. Une base de données locale, RNAstem, a été créée pour faciliter la manipulation et la compilation des données sur les motifs de structure secondaire. Nous avons analysé toutes les boucles terminales et internes ainsi que les « bulges » et nous avons calculé un score d’abondance qui nous a permis d’étudier la fréquence de ces motifs. Tout en minimisant le biais de la surreprésentation de certaines classes d’ARN telles que l’ARN ribosomal, l’analyse des scores a permis de caractériser les motifs rares pour chacune des catégories d’ARN en plus de confirmer des motifs communs comme les boucles de type GNRA ou UNCG. Nous avons identifié des motifs abondants qui n’ont pas été étudiés auparavant tels que la « tetraloop » UUUU. En analysant le contenu de ces motifs en nucléotides, nous avons remarqué que ces régions simples brins contiennent beaucoup plus de nucléotides A et U. Enfin, nous avons exploré la possibilité d’utiliser ces scores pour la conception d’un filtre qui permettrait d’accélérer la recherche de nouveaux ARN non-codants. Nous avons développé un système de scores, RNAscore, qui permet d’évaluer un ARN en se basant sur son contenu en motifs et nous avons testé son applicabilité avec différents types de contrôles.