Structure and dynamics of a gene network model incorporating small RNAs

As advances in molecular biology continue to reveal additional layers of complexity in gene regulation, computational models need to incorporate additional features to explore the implications of new theories and hypotheses. It has recently been suggested that eukaryotic organisms owe their phenotypic complexity and diversity to the exploitation of small RNAs as signalling molecules. Previous models of genetic systems are, for several reasons, inadequate to investigate this theory. In this study, we present an artificial genome model of genetic regulatory networks based upon previous work by Torsten Reil, and demonstrate how this model generates networks with biologically plausible structural and dynamic properties. We also extend the model to explore the implications of incorporating regulation by small RNA molecules in a gene network. We demonstrate how, using these signals, highly connected networks can display dynamics that are more stable than expected given their level of connectivity.

Phenotypic effect of mutations in evolving populations of RNA molecules

Background. The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework. Results. We investigate the effect of microscopic mutations on the phenotype of RNA molecules during their in silico evolution and adaptation. We calculate the distribution of the effects of mutations on fitness, the relative fractions of beneficial and deleterious mutations and the corresponding selection coefficients for populations evolving under different mutation rates. Three different situations are explored: the mutation-selection equilibrium (optimized population) in three different fitness landscapes, the dynamics during adaptation towards a goal structure (adapting population), and the behavior under periodic population bottlenecks (perturbed population). Conclusions. The ratio between the number of beneficial and deleterious mutations experienced by a population of RNA sequences increases with the value of the mutation rate µ at which evolution proceeds. In contrast, the selective value of mutations remains almost constant, independent of µ, indicating that adaptation occurs through an increase in the amount of beneficial mutations, with little variations in the average effect they have on fitness. Statistical analyses of the distribution of fitness effects reveal that small effects, either beneficial or deleterious, are well described by a Pareto distribution. These results are robust under changes in the fitness landscape, remarkably when, in addition to selecting a target secondary structure, specific subsequences or low-energy folds are required. A population perturbed by bottlenecks behaves similarly to an adapting population, struggling to return to the optimized state. Whether it can survive in the long run or whether it goes extinct depends critically on the length of the time interval between bottlenecks. © 2010 Stich et al; licensee BioMed Central Ltd.

Binding of small molecules to DNA junctions

Previous results in our laboratory suggest that the (CG) 4 segments whether present in a right-handed or a left-handed conformation form distinctive junctions with adjacent random sequences. These junctions and their associated sequences have unique structural and thermodynamic properties that may be recognized by DNA-binding molecules. This study probes these sequences by using the following small ligands: actinomycin D, 1,4-bis(((di(aminoethyl)amino)ethyl)amino)anthracene-9,10-dione, ametantrone, and tris(phenanthroline)ruthenium (II). These ligands may recognize the distinctive features associated to the (CG)4 segment and its junctions and thus interact preferentially near these sequences. Restriction enzyme inhibition assays were used to determine whether or not binding interactions took place, and to approximate locations of these interactions. These binding studies are first carried out using two small synthetic oligomers BZ-III and BZ-IV. The (5meCG)4 segment present in BZ-III adopts the Z-conformation in the presence of 50 m M Co(NH3)63+. In BZ-IV, the unmethylated (CG)4 segment changes to a non-B conformation in the presence of 50 m M Co(NH3)63+. BZ-IV, containing the (CG)4 segment, was inserted into a clone plasmid then digested with the restriction enzyme Hinf I to produce a larger fragment that contains the (CG)4 segment. The results obtained on the small oligomers and on the larger fragment for restriction enzyme Mbo I indicate that 1,4-bis(((di(aminoethyl)amino)ethyl)amino)anthracene-9,10-dione binds more efficiently at or near the (CG)4 segment. Restriction enzymes EcoRV, Sac I and Not I with cleavage sites upstream and downstream of the (CG)4 insert were used to further localize binding interactions in the vicinity of the (CG)4 insert. RNA polymerase activity was studied in a plasmid which contained the (CG)4 insert downstream from the promoter sites of SP6 and T7 RNA polymerases. Activities of these two polymerases were studied in the presence of each one of the ligands used throughout the study. Only actinomycin D and spider, which bind at or near the (CG)4 segment, alter the activities of SP6 and T7 RNA polymerases. Surprisingly, enhancement of polymerase activity was observed in the presence of very low concentrations of actinomycin D. These results suggest that the conformational features of (CG) segments may serve in regulatory functions of DNA. ^

The role of splicing factors and small nuclear RNAs in spliceosomal formation

Protein coding genes are comprised of protein-coding exons and non-protein-coding introns. The process of splicing involves removal of the introns and joining of the exons to form a mature messenger RNA, which subsequently undergoes translation into polypeptide. The spliceosome is a large, RNA/protein assembly of five small nuclear RNAs as well as over 300 proteins, which catalyzes intron removal and exon ligation. The selection of specific exons for inclusion in the mature messenger RNA is spatiotemporally regulated and results in production of an enormous diversity of polypeptides from a single gene locus. This phenomenon, known as alternative splicing, is regulated, in part, by protein splicing factors, which target the spliceosome to exon/intron boundaries. The first part of my dissertation (Chapters II and III) focuses on the discovery and characterization of the 45 kilodalton FK506 binding protein (FKBP45), which I discovered in the silk moth, Bombyx mori, as a U1 small nuclear RNA binding protein. This protein family binds the immunosuppressants FK506 and rapamycin and contains peptidyl-prolyl cis-trans isomerase activity, which converts polypeptides from cis to trans about a proline residue. This is the first time that an FKBP has been identified in the spliceosome. The second section of my dissertation (Chapters IV, V, VI and VII) is an investigation of the potential role of small nuclear RNA sequence variants in the control of splicing. I identified 46 copies of small nuclear RNAs in the 6X whole genome shotgun of the Bombyx mori p50T strain. These variants may play a role in differential binding of specific proteins that mediate alternative splicing. Along these lines, further investigation of U2 snRNA sequence variants in Bombyx mori demonstrated that some U2 snRNAs preferentially assemble into high molecular weight spliceosomal complexes over others. Expression of snRNA variants may represent another mechanism by which the cell is able to fine tune the splicing process.

The Role of Splicing Factors and Small Nuclear RNAS in Spliceosomal Formation

Protein coding genes are comprised of protein-coding exons and non-protein-coding introns. The process of splicing involves removal of the introns and joining of the exons to form a mature messenger RNA, which subsequently undergoes translation into polypeptide. The spliceosome is a large, RNA/protein assembly of five small nuclear RNAs as well as over 300 proteins, which catalyzes intron removal and exon ligation. The selection of specific exons for inclusion in the mature messenger RNA is spatio-temporally regulated and results in production of an enormous diversity of polypeptides from a single gene locus. This phenomenon, known as alternative splicing, is regulated, in part, by protein splicing factors, which target the spliceosome to exon/intron boundaries. The first part of my dissertation (Chapters II and III) focuses on the discovery and characterization of the 45 kilodalton FK506 binding protein (FKBP45), which I discovered in the silk moth, Bombyx mori, as a U1 small nuclear RNA binding protein. This protein family binds the immunosuppressants FK506 and rapamycin and contains peptidyl-prolyl cis-trans isomerase activity, which converts polypeptides from cis to trans about a proline residue. This is the first time that an FKBP has been identified in the spliceosome. The second section of my dissertation (Chapters IV, V, VI and VII) is an investigation of the potential role of small nuclear RNA sequence variants in the control of splicing. I identified 46 copies of small nuclear RNAs in the 6X whole genome shotgun of the Bombyx mori p50T strain. These variants may play a role in differential binding of specific proteins that mediate alternative splicing. Along these lines, further investigation of U2 snRNA sequence variants in Bombyx mori demonstrated that some U2 snRNAs preferentially assemble into high molecular weight spliceosomal complexes over others. Expression of snRNA variants may represent another mechanism by which the cell is able to fine tune the splicing process.

FLIP: a targetable mediator of resistance to radiation in non-small cell lung cancer

Resistance to radiotherapy due to insufficient cancer cell death is a significant cause of treatment failure in non-small cell lung cancer (NSCLC). The endogenous caspase-8 inhibitor, FLIP, is a critical regulator of cell death that is frequently overexpressed in NSCLC and is an established inhibitor of apoptotic cell death induced via the extrinsic death receptor pathway. Apoptosis induced by ionizing radiation (IR) has been considered to be mediated predominantly via the intrinsic apoptotic pathway; however, we found that IR-induced apoptosis was significantly attenuated in NSCLC cells when caspase-8 was depleted using RNA interference (RNAi), suggesting involvement of the extrinsic apoptosis pathway. Moreover, overexpression of wild-type FLIP, but not a mutant form that cannot bind the critical death receptor adaptor protein FADD, also attenuated IR-induced apoptosis, confirming the importance of the extrinsic apoptotic pathway as a determinant of response to IR in NSCLC. Importantly, when FLIP protein levels were down-regulated by RNAi, IR-induced cell death was significantly enhanced. The clinically relevant histone deacetylase (HDAC) inhibitors vorinostat and entinostat were subsequently found to sensitize a subset of NSCLC cell lines to IR in a manner that was dependent on their ability to suppress FLIP expression and promote activation of caspase-8. Entinostat also enhanced the anti-tumor activity of IR in vivo. Therefore, FLIP down-regulation induced by HDAC inhibitors is a potential clinical strategy to radio-sensitize NSCLC and thereby improve response to radiotherapy. Overall, this study provides the first evidence that pharmacological inhibition of FLIP may improve response of NCSLC to IR.

First complete genome sequence of a capsicum chlorosis tospovirus isolate from Australia with an unusually large S RNA intergenic region

The first complete genome sequence of capsicum chlorosis virus (CaCV) from Australia was determined using a combination of Illumina HiSeq RNA and Sanger sequencing technologies. Australian CaCV had a tripartite genome structure like other CaCV isolates. The large (L) RNA was 8913 nucleotides (nt) in length and contained a single open reading frame (ORF) of 8634 nt encoding a predicted RNA-dependent RNA polymerase (RdRp) in the viral-complementary (vc) sense. The medium (M) and small (S) RNA segments were 4846 and 3944 nt in length, respectively, each containing two non-overlapping ORFs in ambisense orientation, separated by intergenic regions (IGR). The M segment contained ORFs encoding the predicted non-structural movement protein (NSm; 927 nt) and precursor of glycoproteins (GP; 3366 nt) in the viral sense (v) and vc strand, respectively, separated by a 449-nt IGR. The S segment coded for the predicted nucleocapsid (N) protein (828 nt) and non-structural suppressor of silencing protein (NSs; 1320 nt) in the vc and v strand, respectively. The S RNA contained an IGR of 1663 nt, being the largest IGR of all CaCV isolates sequenced so far. Comparison of the Australian CaCV genome with complete CaCV genome sequences from other geographic regions showed highest sequence identity with a Taiwanese isolate. Genome sequence comparisons and phylogeny of all available CaCV isolates provided evidence for at least two highly diverged groups of CaCV isolates that may warrant re-classification of AIT-Thailand and CP-China isolates as unique tospoviruses, separate from CaCV.

Identification and Characterization of New Small Molecule Inhibitors of Picornavirus Replication

The Picornaviridae family consists of positive-strand RNA viruses that are the causative agents of a variety of diseases in humans and animals. Few drugs targeting picornaviruses are available, making the discovery of new antivirals a high priority. Here, we identified and characterized three compounds from a library of kinase inhibitors that block replication of poliovirus, coxsackievirus B3, and encephalomyocarditis virus. The antiviral effect of these compounds is not likely related to their known cellular targets because other inhibitors targeting the same pathways did not inhibit viral replication. Using an in vitro translation-replication system, we showed that these drugs inhibit different stages of the poliovirus life cycle. A4(1) inhibited the formation of a functional replication complex, while E5(1) and E7(2) affected replication after the replication complex had formed. A4(1) demonstrated partial protection from paralysis in a murine model of poliomyelitis. Poliovirus resistant to E7(2) had a single mutation in the 3A protein. This mutation was previously found to confer resistance to enviroxime-like compounds, which target either PI4KIIIβ (major enviroxime-like compounds) or OSBP (minor enviroxime-like compounds), cellular factors involved in lipid metabolism and shown to be important for replication of diverse positive-strand RNA viruses. We classified E7(2) as a minor enviroxime-like compound, because the localization of OSBP changed in the presence of this inhibitor. Interestingly, both E7(2) and major enviroxime-like compound GW5074 interfered with the viral polyprotein processing. Multiple attempts to isolate resistant mutants in the presence of A4(1) or E5(1) were unsuccessful, showing that effective broad-spectrum antivirals could be developed on the basis of these compounds. Studies with these compounds shed light on pathways shared by diverse picornaviruses that could be potential targets for the development of broad-spectrum antiviral drugs.

STRUCTURAL ANALYSIS OF RIBOSOME BINDING ELEMENTS OF THE 3´ UTR IN TWO POSITIVE-STRAND PLANT RNA VIRUSES

Turnip crinkle virus (TCV) and Pea enation mosaic virus (PEMV) are two positive (+)-strand RNA viruses that are used to investigate the regulation of translation and replication due to their small size and simple genomes. Both viruses contain cap-independent translation elements (CITEs) within their 3´ untranslated regions (UTRs) that fold into tRNA-shaped structures (TSS) according to nuclear magnetic resonance and small angle x-ray scattering analysis (TCV) and computational prediction (PEMV). Specifically, the TCV TSS can directly associate with ribosomes and participates in RNA-dependent RNA polymerase (RdRp) binding. The PEMV kissing-loop TSS (kl-TSS) can simultaneously bind to ribosomes and associate with the 5´ UTR of the viral genome. Mutational analysis and chemical structure probing methods provide great insight into the function and secondary structure of the two 3´ CITEs. However, lack of 3-D structural information has limited our understanding of their functional dynamics. Here, I report the folding dynamics for the TCV TSS using optical tweezers (OT), a single molecule technique. My study of the unfolding/folding pathways for the TCV TSS has provided an unexpected unfolding pathway, confirmed the presence of Ψ3 and hairpin elements, and suggested an interconnection between the hairpins and pseudoknots. In addition, this study has demonstrated the importance of the adjacent upstream adenylate-rich sequence for the formation of H4a/Ψ3 along with the contribution of magnesium to the stability of the TCV TSS. In my second project, I report on the structural analysis of the PEMV kl-TSS using NMR and SAXS. This study has re-confirmed the base-pair pattern for the PEMV kl-TSS and the proposed interaction of the PEMV kl-TSS with its interacting partner, hairpin 5H2. The molecular envelope of the kl-TSS built from SAXS analysis suggests the kl-TSS has two functional conformations, one of which has a different shape from the previously predicted tRNA-shaped form. Along with applying biophysical methods to study the structural folding dynamics of RNAs, I have also developed a technique that improves the production of large quantities of recombinant RNAs in vivo for NMR study. In this project, I report using the wild-type and mutant E.coli strains to produce cost-effective, site-specific labeled, recombinant RNAs. This technique was validated with four representative RNAs of different sizes and complexity to produce milligram amounts of RNAs. The benefit of using site-specific labeled RNAs made from E.coli was demonstrated with several NMR techniques.

Rola białka Hfq w oddziaływaniu niekodujących RNA bakterii z sekwencją kodującą mRNA

Wydział Biologii: Instytut Biologii Molekularnej i Biotechnologii

Spliceosomal small nuclear ribonucleoprotein particle trafficking and maturation in the nuclear compartment

We show for the first time that upon injection into the cytoplasm of the oocyte, fluorescein-labeled spliceosomal snRNAs, in the context of functional snRNPs, are targeted to elongating pre-mRNAs. This finding presents us with a novel assay with which to dissect the mechanism by which snRNPs are targeted to nascent pre-mRNA transcripts. Two critical advantages offered by this system are immediately evident. First, it allows us to investigate the mechanisms employed to recruit snRNPs as it actually transpires within the realm of the cell nucleus. Second, it allows a genome-wide analysis of snRNP recruitment to nascent transcripts, and, hence, the conclusions drawn from these studies do not depend on the sequence of any particular promoter or pre-mRNA. Indeed, it is with this assay that we have stumbled upon a most unanticipated discovery: Contrary to the current paradigm, the co-transcriptional recruitment of splicing snRNPs to nascent transcripts is not contingent on their role in splicing in vivo. Based on these and other data, we have constructed a two-step recruitment-loading model wherein snRNPs are first recruited to pre-mRNA transcripts and only then loaded directly onto cis-acting sequences on nascent pre-mRNA. While conducting studies on snRNP trafficking, a new discovery was made. We found that the lampbrush chromosomes could be visualized by light microscopy in vivo, and that these chromosomes have an architecture that is identical with those in formaldehyde treated nuclear spread preparations. Importantly, we now have the first system with which we can examine the dynamic interactions of macromolecules with specific RNA polymerase II transcriptional units in the live nucleus.

Alternate Conformations Regulate Ribosomal Recoding in a Positive-sense RNA Virus

Positive-sense RNA viruses are important animal, plant, insect and bacteria pathogens and constitute the largest group of RNA viruses. Due to the relatively small size of their genomes, these viruses have evolved a variety of non-canonical translation mechanisms to optimize coding capacity expanding their proteome diversity. One such strategy is codon redefinition or recoding. First described in viruses, recoding is a programmed translation event in which codon alterations are context dependent. Recoding takes place in a subset of messenger RNA (mRNAs) with some products reflecting new, and some reflecting standard, meanings. The ratio between the two is both critical and highly regulated. While a variety of recoding mechanisms have been documented, (ribosome shunting, stop-carry on, termination-reinitiation, and translational bypassing), the two most extensively employed by RNA viruses are Programmed Ribosomal Frameshifting (PRF) and Programmed Ribosomal Readthrough (PRT). While both PRT and PRF subvert normal decoding for expression of C-terminal extension products, the former involves an alteration of reading frame, and the latter requires decoding of a non-sense codon. Both processes occur at a low but defined frequency, and both require Recoding Stimulatory Elements (RSE) for regulation and optimum functionality. These stimulatory signals can be embedded in the RNA in the form of sequence or secondary structure, or trans-acting factors outside the mRNA such as proteins or micro RNAs (miRNA). Despite 40+ years of study, the precise mechanisms by which viral RSE mediate ribosome recoding for the synthesis of their proteins, or how the ratio of these products is maintained, is poorly defined. This study reveals that in addition to a long distance RNA:RNA interaction, three alternate conformations and a phylogenetically conserved pseudoknot regulate PRT in the carmovirus Turnip crinkle virus (TCV).

Characterisation of RNase Y in H. pylori: a non-essential enzyme involved in the processing of cag-PAI non coding RNA 1 (CncR1) sRNA

The importance of Helicobacter pylori as a human pathogen is underlined by the plethora of diseases it is responsible for. The capacity of H. pylori to adapt to the restricted host-associated environment andto evade the host immune response largely depends on a streamlined signalling network. The peculiar H. pylori small genome size combined with its paucity of transcriptional regulators highlights the relevance of post-transcriptional regulatory mechanisms as small non-coding RNAs (sRNAs). However, among the 8 RNases represented in H. pylori genome, a regulator guiding sRNAs metabolism is still not well studied. We investigated for the first time the physiological role in H. pylori G27 strain of the RNase Y enzyme. In the first line of research we provide a comprehensive characterization of the RNase Y activity by analysing its genomic organization and the factors that orchestrate its expression. Then, based on bioinformatic prediction models, we depict the most relevant determinants of RNase Y function, demonstrating a correlation of both structure and domain organization with orthologues represented in Gram-positive bacteria. To unveil the post-transcriptional regulatory effect exerted by the RNase Y, we compared the transcriptome of an RNase Y knock-out mutant to the parental wild type strain by RNA-seq approach. In the second line of research we characterized the activity of this single strand specific endoribonuclease on cag-PAI non coding RNA 1 (CncR1) sRNA. We found that deletion or inactivation of RNase Y led to the accumulation of a 3’-extended CncR1 (CncR1-L) transcript over time. Moreover, beneath its increased half-life, CncR1-L resembled a CncR1 inactive phenotype. Finally, we focused on the characterization of the in vivo interactome of CncR1. We set up a preliminary MS2-affinity purification coupled with RNA-sequencing (MAPS) approach and we evaluated the enrichment of specific targets, demonstrating the suitability of the technique in the H. pylori G27 strain.

Egfr Activating Mutations And Their Association With Response To Platinum-doublet Chemotherapy In Brazilian Non-small Cell Lung Cancer Patients.

The aim of the study was to analyze the frequency of epidermal growth factor receptor (EGFR) mutations in Brazilian non-small cell lung cancer patients and to correlate these mutations with response to benefit of platinum-based chemotherapy in non-small cell lung cancer (NSCLC). Our cohort consisted of prospective patients with NSCLCs who received chemotherapy (platinum derivates plus paclitaxel) at the [UNICAMP], Brazil. EGFR exons 18-21 were analyzed in tumor-derived DNA. Fifty patients were included in the study (25 with adenocarcinoma). EGFR mutations were identified in 6/50 (12 %) NSCLCs and in 6/25 (24 %) adenocarcinomas; representing the frequency of EGFR mutations in a mostly self-reported White (82.0 %) southeastern Brazilian population of NSCLCs. Patients with NSCLCs harboring EGFR exon 19 deletions or the exon 21 L858R mutation were found to have a higher chance of response to platinum-paclitaxel (OR 9.67 [95 % CI 1.03-90.41], p = 0.047). We report the frequency of EGFR activating mutations in a typical southeastern Brazilian population with NSCLC, which are similar to that of other countries with Western European ethnicity. EGFR mutations seem to be predictive of a response to platinum-paclitaxel, and additional studies are needed to confirm or refute this relationship.

Neurodevelopment Of Small-for-gestational Age Infants: Behavioral Aspects In First Year.

The purpose of this study was to compare the behavior of full-term small-for-gestational age (SGA) with full-term appropriate-for gestational age (AGA) infants in the first year of life. We prospectively evaluated 68 infants in the 2nd month, 67 in the 6th month and 69 in the 12th month. The Bayley Scales of Infant Development-II were used, with emphasis on the Behavior Rating Scale (BRS). The groups were similar concerning the item interest in test materials and stimuli; there was a trend toward differences in the items negative affect, hypersensitivity to test materials and adaptation to change in test materials. The mean of Raw Score was significantly lower for the SGA group in the items predominant state, liability of state of arousal, positive affect, soothability when upset, energy, exploration of objects and surroundings, orientation toward examiner. A lower BRS score was associated with the SGA group in the 2nd month.