Web: http://www.ima.umn.edu | Email: ima-staff@ima.umn.edu | Telephone: (612) 624-6066 | Fax: (612) 626-7370
Additional newsletters available at http://www.ima.umn.edu/newsletters

IMA Newsletter #372

October 2007

2007-2008 Program

Mathematics of Molecular and Cellular Biology

See http://www.ima.umn.edu/2007-2008 for a full description of the 2007-2008 program on Mathematics of Molecular and Cellular Biology.

News and Notes

Math video is a winner in Science & Visualization Challenge

Möbius Transformations Revealed, a short video by IMA director Douglas Arnold and colleague Jonathan Rogness, was awarded honorable mention in the National Science Foundation's 2007 Science and Visualization Challenge. It is featured in the current issue of the journal Science along with the other winning entries. The video depicts the beauty of Möbius transformations and shows how moving to a higher dimension reveals their essential unity. The video has already been watched by more than 35,000 viewers on YouTube. It can be found at www.ima.umn.edu/~arnold/moebius/ .

IMA is seeking a new associate director:

The IMA is looking for a new associate director to begin September 1, 2008.

Math Matters public lecture series starts another season

Berkeley professor and former IMA postdoc Bernd Sturmfels will open the 2007–8 season of the IMA Math Matters public lecture series with a talk entitled Algebra, statistics, computation, and biology.

NSF Math Institutes facilitating CDI proposals

The National Science Foundation has announced a major new initiative beginning in 2008 on Cyber-enabled Discovery and Innovation (CDI). The NSF Mathematical Sciences Institutes are making a coordinated effort to inform the mathematical community about this initiative and to facilitate the development of research proposals. As part of this effort the IMA is supporting visits for project planning.

IMA Events

IMA Annual Program Year Workshop

RNA in Biology, Bioengineering and Nanotechnology

October 29 - November 2, 2007

Organizers: Tamar Schlick (New York University), Eric Westhof (Université de Strasbourg I (Louis Pasteur))
Schedule

Tuesday, October 2

11:15a-12:15pFast food, square magic, and polyhedraBenjamin Nill (Freie Universität Berlin)Lind Hall 409

Wednesday, October 3

11:15a-12:15pThe shape of knotsKenneth C. Millett (University of California)Lind Hall 409 MMCB

Tuesday, October 9

11:15a-12:15pModeling multicellularity: from cell rheology to gastrulationTimothy Newman (Arizona State University)Lind Hall 409 PS
5:00p-6:30pReceptionLind Hall 400 PUB10.9.07
7:00p-8:15pAlgebra, statistics, computation, and biologyBernd Sturmfels (University of California)Willey Hall 125 PUB10.9.07

Wednesday, October 10

11:15a-12:15pCombinatorics of RNA secondary structuresChristine E. Heitsch (Georgia Institute of Technology)Lind Hall 409 MMCB

Tuesday, October 16

11:15a-12:15pEstimating connectivity in neuronal and other networksDuane Nykamp (University of Minnesota)Lind Hall 409 PS

Wednesday, October 17

11:15a-12:15pRNA Folding Energy Landscapes in mRNA Splicing Shi-Jie Chen (University of Missouri)Lind Hall 409 MMCB

Tuesday, October 23

10:00a-11:00aStart Seeing Mathematics!Douglas N. Arnold (University of Minnesota)Hubert H. Humphrey Center's Cowles Auditorium
11:00a-12:00pWho Wants to Be a Mathematician—a math contest with Minnesota high school students Douglas N. Arnold (University of Minnesota)Hubert H. Humphrey Center's Cowles Auditorium
11:15a-12:15pOptimization in mathematical biologyAndrew Stein (University of Minnesota)Lind Hall 409 PS

Wednesday, October 24

11:15a-12:15pThe role of anionic residues on the salt dependence of the binding of aminoacyl-tRNA synthetases to tRNA: A Poisson-Boltzmann analysisMarcia O. Fenley (Florida State University)Lind Hall 409 MMCB

Friday, October 26

11:15a-12:15pModeling DNA unknotting by type II topoisomerasesMariel Vazquez (San Francisco State University)Lind Hall 409 MMCB
1:25p-2:25pModeling HCV antiviralsJeffrey Saltzman (Merck & Co., Inc.)Vincent Hall 570 IPS
3:30p-4:30pChromatin modeling and simulation: A tale of histone tailsTamar Schlick (New York University)101 Walter Library

Monday, October 29

9:30a-10:00aRegistration and coffeeEE/CS 3-176 W10.29-11.2.07
10:00a-10:15aWelcome to the IMADouglas N. Arnold (University of Minnesota)EE/CS 3-180 W10.29-11.2.07
10:15a-10:30aIntroductionTamar Schlick (New York University)
Eric Westhof (Université de Strasbourg I (Louis Pasteur))
EE/CS 3-180 W10.29-11.2.07
10:30a-11:00aComputational methods for RNA secondary structure determinationMichael Zuker (Rensselaer Polytechnic Institute)EE/CS 3-180 W10.29-11.2.07
11:05a-11:35aPrediction of the secondary structure common to two sequences: Free energy minimization and comparative analysisDavid H. Mathews (University of Rochester)EE/CS 3-180 W10.29-11.2.07
11:35a-11:50aDiscussionsEE/CS 3-180 W10.29-11.2.07
11:50a-1:30pLunch W10.29-11.2.07
1:30p-2:00pAnalysis and design of nucleic acid devicesNiles A. Pierce (California Institute of Technology)EE/CS 3-180 W10.29-11.2.07
2:05p-2:20pA continuous probabilistic model of local RNA 3-D structureJes Frellsen (University of Copenhagen)EE/CS 3-180 W10.29-11.2.07
2:25p-2:40pRNA matrices and RNA secondary structuresAsamoah Nkwanta (Morgan State University)EE/CS 3-180 W10.29-11.2.07
2:45p-3:00pRNA folding during transcription facilitated by non-native structuresTao Pan (University of Chicago)EE/CS 3-180 W10.29-11.2.07
3:05p-3:20pDiscussionsEE/CS 3-180 W10.29-11.2.07
3:30p-3:45pGroup Photo W10.29-11.2.07
4:00p-5:30pReception and Poster SessionLind Hall 400 W10.29-11.2.07
Efficient parameter estimation for RNA secondary structure predictionMirela Andronescu (University of British Columbia)
SwS: a solvation web service for nucleic acidsPascal Auffinger (Université de Strasbourg I (Louis Pasteur))
The Rfam database: we need youAlex Bateman (Wellcome Trust Sanger Institute)
Paul Gardner (Wellcome Trust Sanger Institute)
Utilizing the RNAJunction database for the design of RNA nanostructuresEckart Bindewald (SAIC-Frederick, Inc.)
Riboswitches, RNA conformational switches and prokaryotic gene regulation (with Eva Freyhulta Vincent Moultonb) Peter Clote (Boston College)
Efficient algorithms for probing the RNA mutation landscape (with Waldispuehl, Devadas, Berger) Peter Clote (Boston College)
Binding of aminoglycosidic antibiotics to the oligonucleotide A-site model Maciej Dlugosz (University of Warsaw)
RNA dinucleotide step parametersMauricio Esguerra (Rutgers University)
A continuous probabilistic model of local RNA 3-D structureJes Frellsen (University of Copenhagen)
Designing structured RNA pools for in vitro selection of RNAsHin Hark Gan (New York University)
Functional classification of all non-coding microbial sequences through phylogenetic profilingDaniel Gautheret (Université de Paris XI (Paris-Sud))
Multi scale simulation of RNA catalytic activityGeorge Madalin Giambasu (University of Minnesota)
Local pairwise structural RNA alignments by pruning of the dynamical programming matrixJan Gorodkin (University of Copenhagen)
RNA pseudoknotted secondary structure prediction using hierarchical foldingHosna Jabbari (University of British Columbia)
Determining functional conformations of two HDV III strainsWojciech (Voytek) Kasprzak (National Cancer Institute)
Finding additional functional elements in essential RNA sites: not conserved, but not unimportantRob Knight (University of Colorado)
Annotated tertiary interaction motifs in RNA structuresChristian Laing (New York University)
MASTR: Simultaneous multiple alignment and structure prediction of non-coding RNAs using simulated annealingStinus Lindgreen (University of Copenhagen)
Computational models for RNA silencing pathways under time-dependent transgene transcription ratesRoderick Melnik (Wilfrid Laurier University)
RNA matrices and RNA secondary structuresAsamoah Nkwanta (Morgan State University)
RNA folding during transcription facilitated by non-native structuresTao Pan (University of Chicago)
Genomic identification of structural RNAs using phylo-SCFGs Jakob Skou Pedersen (University of Copenhagen)
Locomotif: from graphical motif description to RNA motif searchJens Reeder (Universität Bielefeld)
Feyman diagrams, RNA folding, and the transition polynomialYongwu Rong ()
Computational comparative genomics for discovery of cis-regulatory RNAs in bacteria Walter Larry Ruzzo (University of Washington)
Collective properties of evolving populations of RNA moleculesMichael Stich (Instituto Nacional de Tecnica Aeroespacial)
On the design of oligos for gene synthesisChris Thachuk (University of British Columbia)
Comparative genomics beyond sequence based alignments: RNA structures in the ENCODE regionsElfar Torarinsson (University of Copenhagen)
Efficient algorithms for probing the RNA mutation landscape and prediction of deleterious mutationsJerome Waldispuhl (Massachusetts Institute of Technology)
Improved RNA gene predictions through dinucleotide controlled randomization of multiple sequence alignmentsStefan Washietl (Universität Wien)
Estimating the fraction of non-coding RNAs in mammalian transcriptomesYurong Xin (New York University)

Tuesday, October 30

10:45a-11:45aCoffee and TreatsLind Hall 400 W10.29-11.2.07
1:00p-1:30pCoffeeEE/CS 3-176 W10.29-11.2.07
1:30p-2:00pAnalysis, prediction, and design of viral RNA secondary structuresChristine E. Heitsch (Georgia Institute of Technology)EE/CS 3-180 W10.29-11.2.07
2:05p-2:35pARTS and DARTS: A method and database for exploring RNA tertiary structuresRuth Nussinov (National Cancer Institute)EE/CS 3-180 W10.29-11.2.07
2:35p-2:50pDiscussionsEE/CS 3-180 W10.29-11.2.07
2:50p-3:20pImaging RNA structures and folding intermediates using electron cryo-microscopyTobin R. Sosnick (University of Chicago)EE/CS 3-180 W10.29-11.2.07
3:20p-3:50pCoffeeEE/CS 3-176 W10.29-11.2.07
3:50p-4:20pPrediction of RNA-RNA interactionsIvo L. Hofacker (Universität Wien)EE/CS 3-180 W10.29-11.2.07
4:20p-4:35pDiscussions W10.29-11.2.07
4:35p-4:50pGenomic identification of structural RNAs using phylo-SCFGs Jakob Skou Pedersen (University of Copenhagen)EE/CS 3-180 W10.29-11.2.07
4:55p-5:10pLocomotif: from graphical motif description to RNA motif searchJens Reeder (Universität Bielefeld)EE/CS 3-180 W10.29-11.2.07
5:15p-5:30pDesigning Structured RNA Pools for RNA In Vitro Selection-- Hin Hark Gan, Namhee Kim and Tamar SchlickTamar Schlick (New York University)EE/CS 3-180 W10.29-11.2.07
5:30p-5:45pDiscussionsEE/CS 3-180 W10.29-11.2.07

Wednesday, October 31

9:30a-10:00aCoffeeEE/CS 3-176 W10.29-11.2.07
10:00a-10:30aTheory and application of a novel RNA folding approach based on nucleotide cyclic motifsFrancois Major (University of Montreal)EE/CS 3-180 W10.29-11.2.07
10:35a-11:05aStructure-neutral RNA substitutions from 3D structure alignments and 3D motif searchNeocles B. Leontis (Bowling Green State University)EE/CS 3-180 W10.29-11.2.07
11:05a-11:20aDiscussionsEE/CS 3-180 W10.29-11.2.07
11:20a-11:50aCoffeeEE/CS 3-176 W10.29-11.2.07
11:50a-12:20pUsing RNA 3D structure data in SCFG/MRF models to do sequence alignment and motif inferenceCraig L. Zirbel (Bowling Green State University)EE/CS 3-180 W10.29-11.2.07
12:25p-12:40pFeynman diagrams, RNA folding, and the transition polynomialYongwu Rong (George Washington University)EE/CS 3-180 W10.29-11.2.07
12:45p-1:00pUltraconserved nonsense: Pervasive unproductive splicing of SR proteins associated with exceptionally conserved DNA elements -- a bizarre prevalent mode of gene regulationSteven Brenner (University of California)EE/CS 3-180 W10.29-11.2.07
1:05p-1:20p Architecture and reactivity of RNAEric Westhof (Université de Strasbourg I (Louis Pasteur))EE/CS 3-180 W10.29-11.2.07
1:20p-1:30pDiscussionsEE/CS 3-180 W10.29-11.2.07
7:00p-9:00pGroup Event Guthrie Theater818 South 2nd Street, Minneapolis W10.29-11.2.07

Thursday, November 1

10:45a-11:45aCoffee and TreatsLind Hall 400 W10.29-11.2.07
1:00p-1:30pCoffeeEE/CS 3-176 W10.29-11.2.07
1:30p-2:00pHow RNA tells right from wrong: base pairs, tertiary interactions, and counterions in RNA foldingSarah Woodson (Johns Hopkins University)EE/CS 3-180 W10.29-11.2.07
2:05p-2:35pStructure, dynamics and catalytic mechanisms of two ribozymesDavid M.J. Lilley (University of Dundee)EE/CS 3-180 W10.29-11.2.07
2:35p-2:50pDiscussionsEE/CS 3-180 W10.29-11.2.07
2:50p-3:20pCoffeeEE/CS 3-176 W10.29-11.2.07
3:20p-3:50pA topological classification of RNA foldsHenri Orland (Commissariat à l'Énergie Atomique Saclay (CEA))EE/CS 3-180 W10.29-11.2.07
3:55p-4:25pExploring the energy landscape of RNADevarajan Thirumalai (University of Maryland)EE/CS 3-180 W10.29-11.2.07
4:25p-4:40pDiscussionsEE/CS 3-180 W10.29-11.2.07
5:00p-6:30pMath Matters public lecture receptionLind Hall 400 W10.29-11.2.07
7:00p-8:15pMath Matters public lecture: U.S. premier screening of the film "Achieving the unachievable" with the film's writer/director Jean Bergeron (Laval University)Willey Hall 125 W10.29-11.2.07

Friday, November 2

9:30a-10:00aCoffee W10.29-11.2.07
10:00a-10:30aComputational approaches to RNA nanodesignBruce A. Shapiro (National Cancer Institute)EE/CS 3-180 W10.29-11.2.07
10:35a-11:05aRNA tertiary structure as a proto-language for nano-constructionLuc Jaeger (University of California)EE/CS 3-180 W10.29-11.2.07
11:05a-11:35aRNA / RNP synthetic biology Tan Inoue (Kyoto University)EE/CS 3-180 W10.29-11.2.07
11:35a-11:50aDiscussionsEE/CS 3-180 W10.29-11.2.07
11:50a-12:00pConcluding remarksTamar Schlick (New York University)
Eric Westhof (Université de Strasbourg I (Louis Pasteur))
EE/CS 3-180 W10.29-11.2.07
Abstracts
Mirela Andronescu (University of British Columbia) Efficient parameter estimation for RNA secondary structure prediction
Abstract: Joint work with Anne Condon, Holger H. Hoos, David H. Mathews, and Kevin P. Murphy. Motivation: Accurate prediction of RNA secondary structure from the base sequence is an unsolved computational challenge. The accuracy of predictions made by free energy minimization is limited by the quality of the energy parameters in the underlying free energy model. The most widely used model, the Turner99 model, has hundreds of parameters, and so a robust parameter estimation scheme should efficiently handle large data sets with thousands of structures. Moreover, the estimation scheme should also be trained using available experimental free energy data in addition to structural data. Results: In this work, we present constraint generation (CG), the first computational approach to RNA free energy parameter estimation that can be efficiently trained on large sets of structural as well as thermodynamic data. Our constraint generation approach employs a novel iterative scheme, whereby the energy values are first computed as the solution to a constrained optimization problem. Then the newly-computed energy parameters are used to update the constraints on the optimization function, so as to better optimize the energy parameters in the next iteration. Using our method on biologically sound data, we obtain revised parameters for the Turner99 energy model. We show that by using our new parameters, we obtain significant improvements in prediction accuracy over current state-of-the-art methods. Reference:
Mirela Andronescu, Anne Condon, Holger H. Hoos, David H. Mathews, and Kevin P. Murphy, Efficient parameter estimation for RNA secondary structure prediction, Bioinformatics. 2007 Jul 1;23(13):i19-28.
Douglas N. Arnold (University of Minnesota) Who Wants to Be a Mathematician—a math contest with Minnesota high school students
Abstract: The American Mathematical Society is sponsoring Who Wants to Be a Mathematician—a math contest with Minnesota high school students competing for cash and prizes. The event will be hosted by the Institute for Mathematics and its Applications (IMA) at the University of Minnesota's Hubert H. Humphrey Center's Cowles Auditorium. The student contestants can win up to $2000 in cash by answering multiple choice mathematics questions. Contestants for the game are: * Rohit Agrawal (Wayzata) * Nathan Bertram (Becker) * Alex Fish (Bloomington Jefferson) * Kyle Husebye (Stillwater) * Paul Kravik (Minnehaha Academy) * Perihan MacDonald (St. Paul Central) * Calli Thompson (Stillwater) * Justin Wheeler (Hill-Murray) At 10:00 a.m., Douglas N. Arnold, Director of the Institute for Mathematics and its Applications, will give a talk entitled Start Seeing Math! about visualization of complex numbers and functions. Who Wants to be a Mathematician is a game of skill and excitement patterned on the TV show Who Wants To Be A Millionaire. The game rewards math performance while at the same time having fun. The American Mathematical Society has been performing the game for over six years at locations from Boston to Honolulu, but this is the first chance for Minnesota high school students to play. The contest, which is free and open to the public, will last about one hour. Contact either Douglas N. Arnold at 162-264-6066 or Mike Breen, AMS Public Awareness Officer at 800-321-4267 for more information.
Douglas N. Arnold (University of Minnesota) Start Seeing Mathematics!
Abstract: This special lecture is about visualization of complex numbers and functions.
Pascal Auffinger (Université de Strasbourg I (Louis Pasteur)) SwS: a solvation web service for nucleic acids
Abstract: Modeling accurately the solvation of nucleic acid systems is an important issue since it has been shown that water, together with the surrounding ionic atmosphere, is an essential component of RNA and DNA structure. A new web service, called SwS ("Solvation web Service" for nucleic acids), will be presented. This web service, based on the nucleic acid structures contained in the NDB, is devoted to the statistical analysis of the first solvation shell of important structural fragments and has been developed to allow accurate comparisons between theoretical (molecular dynamics simulations) and experimental (x-ray) data and to better understand molecular recognition phenomena involving water and ions. Such data will also, at a more subtle level, improve our views on assembly rules of tertiary structural motifs of nucleic acids.
Alex Bateman (Wellcome Trust Sanger Institute), Paul Gardner (Wellcome Trust Sanger Institute) The Rfam database: we need you
Abstract: The Rfam database is a collection of multiple sequence alignments and covariance models representing many common non-coding RNA gene (ncRNA) families. Rfam aims to facilitate the identification and classification of new members of known sequence families, and distributes annotation of ncRNAs in over 200 complete genome sequences. Rfam release 8.0 contains 574 ncRNA families (including 427 bona fide RNA genes, and 145 regulatory elements). For each family we provide predicted secondary structures, multiple sequence alignments, species distribution, annotation and links to other external specialised resources. All our data is available and searchable online or for download and local installation.
Jean Bergeron (Laval University) Math Matters public lecture: U.S. premier screening of the film "Achieving the unachievable" with the film's writer/director
Abstract: M.C. Escher is among the most mathematical of artists. In 1956 he challenged the laws of perspective with his graphic Print Gallery, and found himself trapped by an impossible barrier. His uncompleted master-piece quickly became the most puzzling enigma of modern art, for both artists and scientists. Half a century later, mathematician Hendrik Lenstra took everyone by surprise by drawing a fantastic bridge between the intuition of the artist and his own, and completed Escher's work mathematically. This story is presented in the 52 minute film Achieving the Unachievable by documentary filmmaker Jean Bergeron. After the screening, the film's U.S. premier, Bergeron will be available to answer questions.
Eckart Bindewald (SAIC-Frederick, Inc.) Utilizing the RNAJunction database for the design of RNA nanostructures
Abstract: Joint work with Wojciech Kasprzak1, Mary O’Connor2, Brett Boyle2 and Bruce A. Shapiro2. 1 Basic Research Program, SAIC-Frederick, Inc., NCI Frederick, Frederick, Maryland, USA
2 Center for Cancer Research Nanobiology Program, NCI Frederick, Frederick, Maryland, USA. We are presenting RNAJunction, which is a database containing extracted and annotated 3D coordinate data of RNA junctions, kissing loops, internal loops and bulges. The database contains more than 12000 structural elements and allows web-based querying by sequence, type and PDB information. The database allows searching by geometric constraints (inter-helix angles); this is useful for the design of RNA nanostructures. We show how these structural elements can be utilized to generate ring structures and other complexes using the NanoTiler software. We present five different approaches for assembling RNA complexes from building blocks. Several examples of automatically generated computational RNA models are presented. Funded in part by DHHS #N01-CO-12400.
Steven Brenner (University of California) Ultraconserved nonsense: Pervasive unproductive splicing of SR proteins associated with exceptionally conserved DNA elements -- a bizarre prevalent mode of gene regulation
Abstract: Nonsense-mediated mRNA decay (NMD) is a cellular RNA surveillance system that recognizes transcripts with premature termination codons and degrades them. We previously discovered large numbers of natural alternative splice forms that appear to be targets for NMD, and we speculated that this might be a mode of gene regulation which we termed RUST (regulated unproductive splicing and translation). This seems to be confirmed by our finding that all conserved members of the SR family of splice regulators have an unproductive alternative mRNA isoform targeted for NMD. Strikingly, the splice pattern for each is conserved in mouse and always associated with an ultraconserved or highly-conserved region of ~100 or more nucleotides of perfect identity between human and mouse. Remarkably, this seems to have evolved independently in every one of the genes, suggesting that this is a natural mode of regulation.
Shi-Jie Chen (University of Missouri) RNA Folding Energy Landscapes in mRNA Splicing
Abstract: We develop a statistical mechanical model for RNA/RNA complexes with both intramolecular and intermolecular interactions. As an application of the model, we compute the free energy landscapes, which give the full distribution for all the possible conformations for spliceosomal snRNA complexes. Different snRNA experiments found contrasting structures, our free energy landscape theory shows why these structures emerge and how they compete with each other. In addition, the energy landscapes suggest possible mechanisms for the conformational switches in splicing. The change of the energy landscape shape gives information about the conformational changes. We find multiple (native-like and misfolded) intermediates formed through base-pairing rearrangements in snRNA complexes. Furthermore, the energy landscape gives the stabilities of all the possible (functional) intermediates and such information is directly related to splicing efficiency.
Peter Clote (Boston College) Riboswitches, RNA conformational switches and prokaryotic gene regulation (with Eva Freyhulta Vincent Moultonb)
Abstract: Linnaeus Centre for Bioinformatics, Uppsala University, 75124 Uppsala, Sweden, eva.freyhult@lcb.uu.se, b School of Computing Sciences, University of East Anglia, Norwich, NR4 7TJ, UK, vincent.moulton@cmp.uea.ac.uk, c Department of Biology, Boston College, Chestnut Hill, MA 02467, USA, clote@bc.edu. This work is funded in part by NSF DBI-0543506. Metabolite-sensing 5 -UTR (untranlated regions) of certain mRNAs, called riboswitches, have been discovered to undergo a conformational change upon ligand-binding, which thereby can up- or down-regulate the corresponding protein product. For instance, upon the binding of nucleotide guanine, the G-box riboswitch in the 5 UTR of the XPT gene of Bacillus subtillis un- dergoes a conformational change to create a terminator loop, thereby pre- maturely terminating transcription of the XPT gene. Since XPT is involved in guanine metabolism, this is an example of negative autoregulation by a riboswitch. Although riboswitches have been postulated to be an ancient genetic regulatory system, first developed in bacteria, the remarkable dis- covery of Cheah et al. in Nature 2007 suggests that eukaryotes may have co-opted riboswitches to control alternative splicing of genes. Here we describe a new algorithm RNAbor (Freyhult, Moulton, Clote Bioinformatics 2007) which gives information on possible conformational switches by computing the Boltzmann probability of structural neighbors of a given RNA secondary structure. A secondary structure T of a given RNA sequence s is called a δ-neighbor of S if T and S differ by exactly δ base pairs. RNAbor computes the number (Nδ ), the Boltzmann partition function (Zδ ) and the minimum free energy (MFEδ ) and corresponding structure over the collection of all δ-neighbors of S. This computation is done simultaneously for all δ ≤ m, in run time O(m2 n3 ) and memory O(mn2 ), where n is the sequence length. We apply RNAbor for the detection of possible RNA con- formational switches, and compare RNAbor with an existent switch detection method. We also provide examples of how RNAbor can at times improve the accuracy of secondary structure prediction.
Peter Clote (Boston College) Efficient algorithms for probing the RNA mutation landscape (with Waldispuehl, Devadas, Berger)
Abstract: The diversity and importance of the role played by RNAs in the regulation and development of the cell has now been demonstrated. This broad range of functions is achieved through specific structures which have been (presumably) optimized through evolution. The existence of a well-founded energy function for RNA has enabled accurate ab-initio secondary structure prediction. State-of-the-art methods such as McCaskill, use a statistical mechanics framework based on the computation of the partition function over the canonical ensemble of all possible secondary structures on a given sequence. Unfortunately, these techniques do not permit any modification of the input sequence during their execution and thus cannot investigate the mutation landscape of this sequence.
Maciej Dlugosz (University of Warsaw) Binding of aminoglycosidic antibiotics to the oligonucleotide A-site model
Abstract: Coauthors J. M. Antosiewicz and J. Trylska. Aminoglycosidic antibiotics are anti-bacterial molecules which target the A-site of the small ribosomal subunit. Using Brownian dynamics we simulated the encounter of four different aminoglycosidic antibiotics with their RNA binding site on the ribosome. The considered antibiotics include neamine, neomycin, paromomycin and ribostamycin. They are amine sugar derivatives, composed of 2 to 4 rings, with a positive total charge of +4 to +6e. The influence of structural, electrostatic and hydrodynamic properties of antibiotics on the kinetics of their association with the ribosomal A-site is discussed. Diffusion limited rates of association are computed and their dependence on ionic strength of the surrounding is examined. The mechanism of diffusion towards the RNA and the formation of the encounter complex is analyzed.
Mauricio Esguerra (Rutgers University) RNA dinucleotide step parameters
Abstract: We present a first view of the space of conformations adopted by RNA in the currently best-resolved structure of the large ribosomal subunit using the dinucleotide ‘step’ parameters computed with the 3DNA software. We have explored how the base-step parameters for the 16 possible nucleotide steps of RNA vary in helical vs. non-helical regions.
Jes Frellsen (University of Copenhagen) A continuous probabilistic model of local RNA 3-D structure
Abstract: Joint work with Ida Moltke, Martin Thiim and Thomas Hamelryck (The Bioinformatics Center, University of Copenhagen) So far, the most common approach to modeling local RNA 3-D structure has been to describe the local conformational space as discrete in a non-probabilistic framework. We present an original approach to modeling local RNA 3-D structure, namely a probabilistic model that treats the conformational space as continuous. In our model the backbone dihedral angles and the base dihedral angles are modeled with a Dynamic Bayesian Network using directional statistics. The model assigns a probability distribution to the conformational space and therefore it has numerous applications. It allows for fast probabilistic sampling of locally RNA-like structures and it can therefore be used in RNA 3-D structure prediction, where one of the problems is how to efficiently search through the space of plausible RNA structures. Today, the state-of-the-art method for suggesting plausible RNA structures is based on assembling fragments from libraries. Further, the model can also be used for deriving probabilities of seeing different local structures and it can therefore be used for quality validation of experimentally determined structures.
Jes Frellsen (University of Copenhagen) A continuous probabilistic model of local RNA 3-D structure
Abstract: Joint work with Ida Moltke, Martin Thiim and Thomas Hamelryck (The Bioinformatics Center, University of Copenhagen) So far, the most common approach to modeling local RNA 3-D structure has been to describe the local conformational space as discrete in a non-probabilistic framework. We present an original approach to modeling local RNA 3-D structure, namely a probabilistic model that treats the conformational space as continuous. In our model the backbone dihedral angles and the base dihedral angles are modeled with a Dynamic Bayesian Network using directional statistics. The model assigns a probability distribution to the conformational space and therefore it has numerous applications. It allows for fast probabilistic sampling of locally RNA-like structures and it can therefore be used in RNA 3-D structure prediction, where one of the problems is how to efficiently search through the space of plausible RNA structures. Today, the state-of-the-art method for suggesting plausible RNA structures is based on assembling fragments from libraries. Further, the model can also be used for deriving probabilities of seeing different local structures and it can therefore be used for quality validation of experimentally determined structures.
Hin Hark Gan (New York University) Designing structured RNA pools for in vitro selection of RNAs
Abstract: In vitro selection of RNAs is a versatile experimental technology for discovering novel RNA molecules from randon sequence pools. However, finding complex RNA molecules is difficult because simple motifs dominate in random pools. Thus, engineering sequence pools possessing complex structures could increase the probability of discovering novel RNAs. The mathematical problem of designing structured RNA pools is to optimize the sequence/structure space to yield the structural characteristics of the target pool. We represent experimental pool generation as nucleotide mixing (transition) matrix applied to a starting sequence, and pool structures as RNA graphs. These tools allow us to map regions of RNA sequence space using mixing matrix and their structural distributions. The target structured pool corresponds to an optimal combination of mixing matrices, starting sequences, and associated pool fractions. We show that our pool design approach allows generation of pools with user-defined characteristics, such as proportions of specific target motifs, starting functional sequences, and sequence length. Our pool design method has been automated and made available through the webserver RAGPOOLS (http://rubin2.biomath.nyu.edu) that offers a theoretical companion tool for RNA in vitro selection and related problems. Thus, RAGPOOLS can serve as a guide to researchers who aim to synthesize RNA pools with desired properties and/or perform in silico experiments. References: Kim N, Shin JS, Elmetwaly S, Gan HH, and Schlick T, RAGPOOLS: RNA-As-Graph-Pools A web server for assisting the design of structured RNA pools for i n vitro s election. Bioinformatics 2007 (In Press). Kim N, Gan HH, Schlick T, A computational proposal for designing structured RNA pools for in vitro selection of RNAs.RNA 2007, 13(4):478-92.
Daniel Gautheret (Université de Paris XI (Paris-Sud)) Functional classification of all non-coding microbial sequences through phylogenetic profiling
Abstract: Joint work with Antonin Marchais and Magali Naville (IGM. Bât 400 - Université Paris-Sud - 91405 Orsay cedex – France). Although comparative genomics has been instrumental in the identification of novel non-coding RNA (ncRNA) in model genomes, this technique cannot, in the form it is currently practised, keep up with the pace of genome sequencing. As a result, hundreds of microbial genomes, including entire families of important pathogens, have been left out of the picture in terms of ncRNA function analysis. As ncRNAs play major regulatory and adaptive roles in bacteria, there is an urgent need for innovative computational methods that would permit a quick and efficient detection of ncRNAs in any genome of interest. Here we propose a protocol that exploits the depth of phylogenetic information in all available genomes (with virtually no limitation in the number of species) to produce a functional classification of all ncRNA candidates and other non-coding conserved elements in any target bacterial genome. Our protocol involves a low-stringency screening for intergenic conserved elements (ICEs) in the target genome, followed by the construction of the presence/absence profile of each ICE across the complete bacterial genome collection. All ICEs are then clustered according to the distance of their phylogenetic profiles, as done by Pellegrini et al. [1] for classifying protein genes. A simultaneous clustering of ICEs and ORFs produces a complete classification of coding and non-coding elements in the target genome. We ran this pipeline on E. coli and B. subtilis. In both species, known small RNAs and riboswitches were significantly concentrated (P~10-6) in two or three clusters containing as well many orthologous E. coliand B. subtilis ORFs and ~200 undefined ICEs in each species. Phylogenetic profile clustering is independent of sequence similarity and appears to predict functional ncRNAs with a much higher specificity than comparative sequence analysis. Furthermore, some clusters that lack known ncRNAs show very interesting phylogenetic presence/absence patterns that indicate either horizontal transfers or the emergence of common adaptive non-coding elements in distant bacterial species. Finally, the co-occurrence of ICEs and protein coding genes in the same clusters may constitute an important source of information on ICE/ORF functional relationships. A complete run of our ICE classification pipeline on a bacterial genome only requires a few hours. [1] Pellegrini M, Marcotte EM, Thompson MJ, Eisenberg D, Yeates TO. (1999) Assigning protein functions by comparative genome analysis: protein phylogenetic profiles. Proc. Natl. Acad. Sci. U.S.A. 96:4285-4288
George Madalin Giambasu (University of Minnesota) Multi scale simulation of RNA catalytic activity
Abstract: Joint work with Taisung Lee and Darrin M. York (Department of Chemistry, University of Minnesota). We present a series of multi-scale simulation studies on RNA catalysis. The results of several series of molecular dynamics (MD) and QM/MM simulations on of the full-length hammerhead ribozyme and the L1 Ligase ribozyme are presented. For the hammerhead ribozyme we have used simulations to investigate the role of metal ions and the possible solvent structure in the crystal, and study/predict the mutation effects at the C3 and G8 sites. For the L1 Ligase we have studied the details of a major conformational change prior to the reaction and possible conformations of the ligation site in the reactant state.These simulations (each with a length of 50 to 100 ns, with a total of more than 1.5 ms) are at least one to two orders longer than any previous reported simulations and significant amount of unrevealed insights have been found through our simulations.
Jan Gorodkin (University of Copenhagen) Local pairwise structural RNA alignments by pruning of the dynamical programming matrix
Abstract: Joint with Jakob H. Havgaard and Elfar Torarinsson (Division of Genetics and Bioinformatics, IBHV, University of Copenhagen, Frederiksberg, Denmark). The Sankoff algorithm for simultaneously folding and aligning RNA sequences is computationally very heavy. Recently a number of groups have applied various constraints to lower the computational requirements to reasonable levels. Whereas the original Sankoff algorithm as well as many of the implementations, only conduct global alignments, the FOLDALIGN implementation makes both local and global structural alignment. The most recent version of FOLDALIGN introduces pruning of the dynamical programming matrix as a simple and effective heuristic which lowers the time and memory requirements significantly without lowering the predictive performance. FOLDALIGN is currently one of few Sankoff alogorithms capable of conductiong local alignments while being a practical tool. It has also been used in genome-wide screen for putative RNA structures in corresponding, but unaligned regions between human and mouse. In addition to the pairwise version of FOLDALIGN we have also made a multiple alignment method which either takes the pairwise alignments or McCaskill basepair probability matrices as input. References - Fast pairwise structural RNA alignments by pruning of the dynamical programming matrix. J. H. Havgaard, E. Torarinsson and J. Gorodkin PLoS Computational Biology, in press - Multiple structural alignment and clustering of RNA sequences. E. Torarinsson, J. H. Havgaard and J. Gorodkin Bioinformatics, 23:926-932, 2007. - Thousands of corresponding human and mouse genomic regions unalignable in primary sequence contain common RNA structure. E. Torarinsson, M. Sawera, J. H. Havgaard, M. Fredholm and J. Gorodkin Genome Research, 16:885-889, 2006.
Christine E. Heitsch (Georgia Institute of Technology) Analysis, prediction, and design of viral RNA secondary structures
Abstract: Understanding how biological sequences encode structural and functional information is a fundamental scientific challenge. For RNA viral genomes, the information encoded in the sequence extends well-beyond their protein coding role to the role of intra-sequence base pairing in viral packaging, replication, and gene expression. Working with the Pariacoto virus as a model sequence, we investigate the compatibility of predicted base pairings with the dodecahedral cage known from crystallographic studies. To build a putative secondary structure, we first analyze different possible configurations using a combinatorial model of RNA folding. We give results on the trade-offs among types of loop structures, the asymptotic degree of branching in typical configurations, and the characteristics of stems in "well-determined" substructures. These mathematical results yield insights into the interaction of local and global constraints in RNA secondary structures, and suggest new directions in understanding the folding of RNA viral genomes.
Christine E. Heitsch (Georgia Institute of Technology) Combinatorics of RNA secondary structures
Abstract: Under a suitable abstraction, complex biological problems can reveal surprising mathematical structure. For instance, the planar self-bonding, or nested "secondary structure," of an RNA molecule is nicely represented by a plane tree. Under our plane tree model of RNA folding, we prove combinatorial theorems which yield insight into the coding of structural and functional information in RNA sequences. One result demonstrates the importance of local constraints in specifying a global structure while another characterizes the degree of branching in minimal energy configurations. Additionally, by transitioning between RNA configurations by an appropriate local move, we obtain an isomorphism with the lattice of noncrossing partitions. Thus, the interaction between combinatorics and molecular biology motivates new combinatorial theorems as well as advancing biological applications.
Ivo L. Hofacker (Universität Wien) Prediction of RNA-RNA interactions
Abstract: Most noncoding RNAs exert their function by interacting with other RNA molecules, as in the case of microRNA and their mRNA targets. Recent bioinformatical as well as experimental approaches produce thousands of novel RNA transcripts most of which cannot be annotated. Our best hope for elucidating the function of these novel transcripts is through identification of their interaction partners. In this talk I will review several approaches to predict the structure of two RNAs upon hybridization. I will present a fast method to compute the probability that some region of an RNA is unpaired and thus accessible for inter-molecular interactions (or, equivalently, the free energy needed to open up the site), as well as a method to quickly search for possible hybridization sites. The combination of these two approaches yields a promising approach to identify the RNA interaction partners of noncoding RNAs. The site accessibility is also a potent predictor of siRNA efficacy and can be used to improve microRNA target predictions.
Tan Inoue (Kyoto University) RNA / RNP synthetic biology
Abstract: In general, molecular design of RNA is difficult at the 3D level because of its highly complicated folding process. In the 1990s, biochemical and structural analyses revealed that many functional noncoding natural RNAs are organized into modules and fold into defined 3D structures. Moreover, several commonly used RNA–RNA binding motifs in these RNAs were identified by phylogenetic comparison and high-resolution structural analyses. Consequently, it has become possible to design self-folding RNAs precisely by employing such motifs and mimicking the modular organization of natural RNAs. As one such example, we have investigated the design and construction of a self-folding RNA scaffold consisting of standard doublestranded helices connected by the two RNA–RNA binding motifs. Results indicated that the constructed RNA folds compactly into the designed 3D structure. We have also reported the synthesis and development of an artificial RNA enzyme by installing a reaction site and a catalytic site into the designed RNA scaffold. For medical and biological applications, the goals of our current project are 1) to establish multifunctional RNP molecules with tumor seeking sensors, imaging agents and toxins that kill target cells, and 2) to establish artificial signal transduction systems for regulating function of a cell by employing designed RNA and RNP molecules. The strategy may be applicable to the synthesis and development of a variety of nonnatural functional RNAs with defined 3D structures.
Hosna Jabbari (University of British Columbia) RNA pseudoknotted secondary structure prediction using hierarchical folding
Abstract: Improving the accuracy and efficiency of computational RNA secondary structure prediction is an important challenge, particularly for pseudoknotted secondary structures. We propose a new approach for prediction of pseudoknotted structures, motivated by the hypothesis that RNA structures fold hierarchically, with pseudoknot free pairs forming initially, and pseudoknots forming later so as to minimize energy relative to the initial pseudoknot free structure. Our HFold (Hierarchical Fold) algorithm has O(n3) running time, and can handle a wide range of biological structures, including nested kissing hairpins, which have previously required O(n6) time using traditional minimum free energy approaches. We also report on an experimental evaluation of HFold.
Luc Jaeger (University of California) RNA tertiary structure as a proto-language for nano-construction
Abstract: Common occurrence of many small structural motifs in natural RNA molecules suggests that nature utilize s a vocabulary of sequence patterns to compose structural molecules with sophisticated topologies such as the ribosome and large ribozymes. By careful analysis of sequences and tertiary structures of natural RNAs, 3D RNA modules and their folding and assembly principles are presently gathered for generating the syntax of a proto-language for rational design and prediction of RNA 3D shapes. RNA architectonics refers to the creation of this proto-language and to its use to build new RNAs with self-assembly properties. Recently, RNA architectonics led to the reliable prediction and design of the tertiary structure of several artificial RNA building blocks able to form programmable filaments and 2D RNA arrays at the nano-scale level. As a proof of concept, we also demonstrated that structurally complex RNAs based on a syntax involving a repertoire of several different RNA motifs can self-assemble into complex supra-molecular 3D nano-particles. This studies show that RNA architectonics can be used as a tool to explore and compare the biophysical properties of various RNA tertiary structure motifs that would be otherwise more difficult to investigate in isolation or within their natural context. It also demonstrates that RNA is an ideal medium for sculpting addressable and responsive self-assembling architectures of any desired shapes in the 20 to 50 nm scale. Moreover, it suggests that RNA supra-molecular assembly can potentially lead to the development of highly sophisticated therapeutic nano-devices for biological and medical applications. References 1. Jaeger, L. & Chworos, A. (2006) The Architectonics of Programmable RNA and DNA Nanostructures. Current Opinion in Structural Biology, 16, 531-543. 2. Chworos, A, Severcan, I., Koyfman, A. Y., Wienkam, P., Oroudjev, E., Hansma, H. G. & Jaeger, L. (2004). Building programmable jigsaw puzzles with RNA. Science 306, 2068-2072. 3. Nasalean, L., Baudrey, S., Leontis, N.B. & Jaeger, L. (2006) Controlling RNA self-assembly to form filaments. Nucleic Acids Res. 34, 1381-1392 4. Bates, A.D., Callen B.P., Cooper J.M., Cosstick, R., Geary, C., Glidle, A., Jaeger, L., Pearson, J.L., Proupín-Pérez, M., Xu, C., & Cumming, D.R. S. (2006) Construction and characterization of a gold nanoparticle wire assembled using Mg2+-dependent RNA-RNA interactions. Nanoletters 6, 445-448.
Wojciech (Voytek) Kasprzak (National Cancer Institute) Determining functional conformations of two HDV III strains
Abstract: Joint work with Sarah D. Linnstaedt2, John L. Casey2, and Bruce A. Shapiro3. 1Basic Research Program, SAIC-Frederick, Inc., NCI Frederick, Frederick, MD
2Department of Microbiology and Immunology, Georgetown University Medical Center, Washington, DC
3Center for Cancer Research Nanobiology Program, National Cancer Institute, Frederick, MD Hepatitis Delta virus (HDV) is a sub-viral human pathogen aggravating Hepatitis B virus (HBV) liver infections. The short HDV genome (~1680 nt) is a single stranded, circular RNA encoding only one protein, the hepatitis delta antigen (HDAg). The host enzyme ADAR1 edits the HDV stop codon (UAG) into a tryptophan (W) codon (UGG) enabling expression of the two forms of the protein, short and long, from the same open reading frame. HDAg-S is required for replication, while HDAg-L enables viral particle formation and inhibits replication. The balance between the two forms is crucial and editing must be regulated. We have applied our programs, MPGAfold and StructureLab, to predict and examine the folding coformations/states of an HDV III construct. This construct includes the editing site (amber/W) and has the editing capabilities of the full HDV III. The predicted secondary structure folding dynamics indicates that the HDV III RNA forms a meta-stable branched structure and a stable rod structure. Both were observed in vitro, and the branched structure was identified as the one enabling editing. Computational predictions and the experimental data also indicate that an Ecuadorian strain folds into the editing-capable structures more readily than a Peruvian strain, and we indicate the reasons for the difference. Thus the folding dynamics of HDV III strains appears to strongly influence their RNA editing levels. Funded in part by NCI Contract N01-CO-1240.
Rob Knight (University of Colorado) Finding additional functional elements in essential RNA sites: not conserved, but not unimportant
Abstract: Joint work with Vikas Malaiya, Jana Chocholousova, Matthew Iyer, Irene Majerfeld, and Michael Yarus. Evolutionary conservation has often been used to recover the essential pieces of RNA sites, yet can only reveal elements that are necessary, rather than sufficient, for function. Biochemical studies in several systems, including the hammerhead ribozyme and the purine riboswitch, indicate additional regions, such as loop-loop interactions, that are required for function yet are not phylogenetically conserved. Here we use a minimal motif for binding the amino acid tryptophan to ask the ultimate question of an RNA motif: do we know the essential elements well enough to embed the motif in a random-sequence background and obtain functional molecules? We show the utility of this technique for discovering additional sequence requirements for the motif, in this case the requirement for an unpaired G in a specific range of locations and structures relative to the main loop identified by SELEX, and discuss its implications for calculating the probability of obtaining functional RNAs from random-sequence pools.
Christian Laing (New York University) Annotated tertiary interaction motifs in RNA structures
Abstract: RNA tertiary motifs play an important role in RNA folding. To understand the complex organization of RNA tertiary interactions, we compiled a dataset containing 54 high-resolution RNA crystal structures. Seven RNA tertiary motifs (coaxial helix, A-minor, ribose zipper, pseudoknot, kissing hairpin, tRNA D-loop:T-loop and tetraloop-tetraloop receptor) were searched by different computer programs. For the non-redundant RNA dataset, 605 RNA tertiary interactions were found. Most of these 3D interactions occur in the 16S and 23S rRNAs. Exhaustive search of these motifs reveals diversity of interaction. Correlation between motifs (e.g. pseudoknot or coaxial helix with A-minor) shows that they can form "composite" motifs. These findings may lead to tertiary structure constraints useful for RNA 3D prediction.
Neocles B. Leontis (Bowling Green State University) Structure-neutral RNA substitutions from 3D structure alignments and 3D motif search
Abstract: The function of structured RNA molecules depends on forming the correct 3D structure, so the most significant constraints on their sequences are structural. Structure-disrupting substitutions are selected against during evolution while structurally neutral substitutions can accumulate as populations evolve. The relevant interactions include basepairs, base-stacking, and base-phosphate interactions, all of which can be disrupted by certain substitutions. Further constraints are imposed by interactions with other molecules. The general question we address is how to determine which substitutions are structure-neutral in RNA molecules, at the level of individual bases and base-pairs and at the level of 3D motifs and molecular architectures. The availability of two or more 3D structures of large RNA molecules such as the 16S and 23S rRNAs presents opportunities for exploring this question empirically, once the two structures are appropriately aligned. Detailed examination and comparison of nucleotide-nucleotide interaction geometries provides another avenue for addressing the same question. Finally, some RNA motifs occur multiple times in single structures and in non-homologous positions in other molecules, giving another way to study the neutrality of base substitutions. These three approaches will be described and their results compared. Along the way, we will briefly describe FR3D (“Find RNA 3D”), a set of Matlab programs we have developed to annotate RNA structures and to carry out searches for recurrent RNA motifs.
David M.J. Lilley (University of Dundee) Structure, dynamics and catalytic mechanisms of two ribozymes
Abstract: CR-UK Nucleic Acid Structure Group, MSI/WTB complex, University of Dundee, Dundee DD1 5EH, UK d.m.j.lilley@dundee.ac.uk The nucleolytic ribozymes are catalytic RNA molecules that generate site-specific cleavage by means of a transesterification reaction involving the 2’ and 5’ O atoms. We have made a study of two of these, the hairpin and VS ribozymes. The hairpin ribozyme folds to generate an intimate loop-loop interaction to create the local environment in which catalysis can proceed. By means of FRET we can observe individual hairpin ribozyme molecules as they undergo multiple cycles of cleavage and ligation, and measure the rates of the internal reactions. On average, the cleaved ribozyme undergoes several docking-undocking events before a ligation reaction occurs. On the basis of these experiments, we have explored the role of the nucleobases G8 and A38 in the catalysis. Both cleavage and ligation reactions are pH dependent, corresponding to the titration of a group with pKA = 6.2. We have used a novel ribonucleoside in which these bases are replaced by imidazole to investigate the role of acid-base catalysis in this ribozyme. We observe significant rates of cleavage and ligation, and a bell-shaped pH dependence for both. The VS ribozyme is the largest of the nucleolytic ribozymes, and the only one for which there is no crystal structure. The ribozyme consists of five helical sections organised by two three-way junctions, each of which undergo metal ion-induced folding. Using a ‘divide and conquer’ approach based principally on the analysis of component junctions by FRET, we deduced the global structure of the ribozyme. We have now solved the structure of the complete ribozyme at low resolution using small-angle X-ray scattering in solution. The binding of the substrate stem-loop generates a catalytically-productive interaction with the A730 loop active site. We have identified two critical nucleotides in the catalytic process; A756 within the A730 loop, and G638 in the substrate internal loop. Mutation or functional group substitution of either nucleobase leads to > 1,000-fold impairment of catalytic activity, while leaving the structure and binding to the ribozyme unaltered. The pH dependencies of the rate of cleavage of substrate with guanine, adenine, 2,6-diaminopurine or inosine at position 638 are fully consistent with a mechanism in which G638 and A756 act in concert in general acid-base catalysis. The proposed mechanism of the VS and hairpin ribozymes, together with the manner of the generation of the active sitea and their topology, are strikingly similar each other. This has probably arisen by convergent evolution. M.K. Nahas, T.J. Wilson, S. Hohng, K. Jarvie, D.M.J. Lilley and T. Ha Observation of internal cleavage and ligation reactions of a ribozyme Nature Struct. Molec. Biol. 11, 1107-1113 (2004). Z. Zhao, A. McLeod, S. Harusawa, L. Araki, M. Yamaguchi, T. Kurihara and D. M. J. Lilley Nucleobase participation in ribozyme catalysis. J. Amer. Chem. Soc. 127, 5026-5027 (2005). T. J. Wilson, J. Ouellet, Z. Zhao, S. Harusawa, L. Araki, T. Kurihara and D. M. J. Lilley Nucleobase catalysis in the hairpin ribozyme. RNA 12, 980-987 (2006). T. J. Wilson, A. C. McLeod and D. M. J. Lilley A guanine nucleobase important for catalysis by the VS ribozyme EMBO J. 26, 2489-2500 (2007).
Stinus Lindgreen (University of Copenhagen) MASTR: Simultaneous multiple alignment and structure prediction of non-coding RNAs using simulated annealing
Abstract: Joint work with Anders Krogh, University of Copenhagen and Paul Gardner, Wellcome Trust Sanger Institute. With the growing interest in non-coding RNAs and their function, there is also a growing need for computational tools that can be used to analyze new sequences and predict their secondary structure. For a single RNA sequence, a common approach is to find the minimum free energy conformation. However, by looking at many related sequences, it is possible to incorporate evolutionary information and perform a comparative analysis which can improve the prediction. Preferably, one should perform the multiple alignment of the sequences and the prediction of their common secondary structure simultaneously. We present a novel heuristic method that uses simulated annealing to iteratively improve the sequence alignment and a common secondary structure. This is done in the context of a sampling approach using fairly simple moves, that either change the sequence alignment by moving gaps around or update the structure on the base pair level. The prediction is evaluated using a cost function that combines the log-likelihood of the alignment, the base pair probabilities and a covariation term. The method, implemented in the C++ program MASTR (Multiple Alignment of STructural RNAs), is competitive to other current programs, both in terms of speed, alignment quality and structure quality. References:
Lindgreen S, Gardner PP, Krogh A (2007): "MASTR: Multiple alignment and structure prediction of non-coding RNAs using simulated annealing", Bioinformatics (accepted)
Francois Major (University of Montreal) Theory and application of a novel RNA folding approach based on nucleotide cyclic motifs
Abstract: Joint work with Marc Parisien (Institute for Research in Immunology and Cancer, Department of Computer Science, University of Montreal). We change the classical rationale underlying RNA structure prediction by incorporating the contributions of the non-Watson-Crick base pairs. To do so, we define a new first-order object for representing nucleotide relationships in structured RNAs, which we call nucleotide cyclic motif (NCM) (1). In comparison to the classical stacks of Watson‑Crick base pairs, the properties that make NCMs appealing for structure determination are the facts that: i) the same algorithm can be employed for predicting secondary, tertiary, and 3-D structures; ii) the RNA structural motifs are either made of one or more NCMs (2); iii) the NCMs embrace indistinctly both canonical and non‑canonical base pairs; and, iv) the NCMs precisely designate how any nucleotide in a sequence relates to the others. A structure generator and scoring function has been developed: MC-Fold. We show how MC-Fold, combined to MC-Sym (3), builds RNA 3-D structures from sequence data and, combined to MC-Cons, clusters and aligns RNA family sequences. We show how low-resolution data can be incorporated in the modeling to reach conformational states that are difficult to access by sequence data alone. 1. Lemieux, S. and Major, F. (2006) Automated extraction and classification of RNA tertiary structure cyclic motifs. Nucleic Acids Res., 34, 2340-2346. 2. St-Onge, K., Thibault, P., Hamel, S. and Major, F. (2007) Modeling RNA tertiary structure motifs by graph-grammars. Nucleic Acids Res, 35, 1726-1736. 3. Major, F. (2003) Building Three-Dimensional Ribonucleic Acid Structures. IEEE Comp Science Eng 5:44-53.
David H. Mathews (University of Rochester) Prediction of the secondary structure common to two sequences: Free energy minimization and comparative analysis
Abstract: In this talk, I will discuss the dynamic programming methods for simultaneously predicting secondary structure and alignment for two sequences. This approach was first suggested by Sankoff in 1985. Recently, it has been implemented by several groups, using one or more heuristics to reduce the computational cost. Our implementation, Dynalign, finds the lowest free energy common secondary structure. It uses single sequence secondary structure prediction and sequence alignment data as input to reduce the search space to make the calculation tractable for long sequences.
Roderick Melnik (Wilfrid Laurier University) Computational models for RNA silencing pathways under time-dependent transgene transcription rates
Abstract: Joint work with Jack Yang and Roy Mahapatra. The synthesis of dsRNA is analyzed using a pathway model with amplifications caused by the aberrant RNAs. The transgene influx rates are assumed time-decaying and Gaussian functions of time. The dynamics of the transgene induced RNA silencing is investigated with a system of coupled non-autonomous nonlinear differential equations describing the process phenomenologically. The silencing phenomena are detected after a period of transcription. Important contributions of several parameters, including those leading to bifurcation patterns, are discussed with a series of numerical examples.
Kenneth C. Millett (University of California) The shape of knots
Abstract: What do the knots that occur in natural macromolecules look like? What about those arising in random equilateral spatial polygons such as those used to model DNA or polymers? The asphericity is a number between 1.0 and 0.0 that measures the degree to which the ellipsoid of inertia is more like a sausage (near 1.0), a rugby ball (near 0.5) or, a soccer ball (near 0.0). An analogous measure, the spatial asphericity, is defined by considering the smallest ellipsoid containing the knot. Computer simulations and examples of proteins will be used to illustrate these measures. Via the simulations, we will consider the influence of constrained topology on these measures in comparison to the unconstrained averages. In a continuation of earlier research, we will look at the equilibrium lengths defined by the asphericity, the lengths at with the constrained and unconstrained have equal values, and compare these with those determined earlier. We present computer simulations to examine the differences between the average asphericity of polygons with constrained and unconstrained topology.
Timothy Newman (Arizona State University) Modeling multicellularity: from cell rheology to gastrulation
Abstract: I will present an overview of recent work focused on constructing computer models of developmental systems. In particular, we have constructed an algorithm, called the Subcellular Element Model, which is able to simulate large numbers (thousands) of deformable cells in three dimensions. I will describe the inner workings of the algorithm, and indicate that the method is capable of modeling developmental systems over a wide range of scales - capturing cell visco-elasticity at small scales, and long-ranged coordinated cell movement at large scales. Modeling of primitive streak extension in the chick embryo will be discussed as a concrete application.
Benjamin Nill (Freie Universität Berlin) Fast food, square magic, and polyhedra
Abstract: When considering the possible ways of ordering 200 Chicken McNuggets or dealing with the more classical challenge of finding the number of magic squares with given magic number we are inevitably led to the problem of counting lattice points in lattice polytopes. We are going to describe two modern gems of lattice point enumeration, namely Brion's theorem and Barvinok's algorithm, from an easily accessible viewpoint.
Asamoah Nkwanta (Morgan State University) RNA matrices and RNA secondary structures
Abstract: Two lower-triangular arrays with entries that count RNA secondary structures of a given length are mentioned in this short talk. The array entries also count specific lattice walks. There is a one-to-one correspondence between RNA structures and a subset of the walks. We will discuss various ways in which the walks can be used as a tool to help predict primary RNA sequences.
Asamoah Nkwanta (Morgan State University) RNA matrices and RNA secondary structures
Abstract: Two lower-triangular arrays with entries that count RNA secondary structures of a given length are mentioned in this short talk. The array entries also count specific lattice walks. There is a one-to-one correspondence between RNA structures and a subset of the walks. We will discuss various ways in which the walks can be used as a tool to help predict primary RNA sequences.
Ruth Nussinov (National Cancer Institute) ARTS and DARTS: A method and database for exploring RNA tertiary structures
Abstract: Joint work with Oranit Dror, Mira Avraham, Haim Wolfson (Tel Aviv University and SAIC, NCI-Frederick). An increasing number of non-coding RNAs have recently been discovered as key players in a variety of cellular pathways and pathological processes. Much like proteins,the function of these active RNAs can be inferred from their tertiary (3D) structures. However, in contrast to proteins, the number of tools and databases for 3D structural analysis of RNA is still limited. With the aim to fill this void, we have developed a computational method, named ARTS, for aligning RNA tertiary structures. Given a pair of RNA structures, the method searches for a-priori unknown common substructures. The search is truly three dimensional and irrespective of the order of the nucleotide chain. The detected common substructures are either large global folds or small local tertiary motifs. The method is highly-efficient and was used in a fully automatic framework for clustering all the currently available RNA structures. The result is a database, named DARTS, which reveals the current fold repertoire of solved RNA structures and provides a hierarchical classification for them. Both the method and the database should be useful for structural and functional analysis of RNA. They may shed new light on the evolutionary relationship between RNAs and reveal possible building blocks and functional properties. This publication has been funded in whole or in part with Federal funds from the National Cancer Institute, National Institutes of Health, under contract # NO1-CO-12400.
Duane Nykamp (University of Minnesota) Estimating connectivity in neuronal and other networks
Abstract: Determining the connectivity structure of complex networks is hindered by one's inability to simultaneously measure the activity of all nodes. In experiments probing networks such as gene regulatory networks, computer networks, or neural networks, many hidden nodes could be interacting with the small set of measured nodes and corrupting estimates of connectivity in unknown ways. For example, if a hidden node had connections onto two measured nodes, this common input could introduce correlations among the measured nodes, which might lead one to erroneously infer a connection between the measured nodes. We present a model-based approach to control for such effects of hidden nodes in networks. We demonstrate the promise of this approach via simulations of small networks of neurons driven by a visual stimulus.
Henri Orland (Commissariat à l'Énergie Atomique Saclay (CEA)) A topological classification of RNA folds
Abstract: After reviewing some elementary properties of RNA, we show how the RNA folding problem can be formulated exactly in terms of an NxN matrix field theory. This formulation introduces a classification of RNA structures according to their topological genus. The large N limit of this theory generates the secondary structures of RNA (planar graphs), whereas 1/N corrections are identified as pseudo-knots. We show how the RNA structures can be analyzed in terms of primitive pseudo knots of low genus and how this concept can be included in Monte Carlo calculations to actually predict RNA folds.
Tao Pan (University of Chicago) RNA folding during transcription facilitated by non-native structures
Abstract: RNA folding in the cell occurs during transcription. Expedient RNA folding must avoid the formation of undesirable structures as the nascent RNA emerges from the RNA polymerase. We show that efficient folding during transcription of three conserved non-coding RNAs from E. coli, RNase P RNA, SRP RNA and tmRNA is facilitated by their cognate polymerase pausing at specific locations. These pause sites are located between the upstream and the downstream portions of all the native long-range helices in these non-coding RNAs. In the paused complexes, the nascent RNAs form labile structures that sequester these upstream portions in a manner as to guide folding. Both the pause sites and the secondary structure of the non-native portions of the paused complexes are phylogenetically conserved. Specific pausing-induced structural formation can be a general strategy to facilitate the folding of long-range helices. This polymerase-based mechanism may result in portions of non-coding RNA sequences to be evolutionarily conserved for efficient folding during transcription.
Tao Pan (University of Chicago) RNA folding during transcription facilitated by non-native structures
Abstract: RNA folding in the cell occurs during transcription. Expedient RNA folding must avoid the formation of undesirable structures as the nascent RNA emerges from the RNA polymerase. We show that efficient folding during transcription of three conserved non-coding RNAs from E. coli, RNase P RNA, SRP RNA and tmRNA is facilitated by their cognate polymerase pausing at specific locations. These pause sites are located between the upstream and the downstream portions of all the native long-range helices in these non-coding RNAs. In the paused complexes, the nascent RNAs form labile structures that sequester these upstream portions in a manner as to guide folding. Both the pause sites and the secondary structure of the non-native portions of the paused complexes are phylogenetically conserved. Specific pausing-induced structural formation can be a general strategy to facilitate the folding of long-range helices. This polymerase-based mechanism may result in portions of non-coding RNA sequences to be evolutionarily conserved for efficient folding during transcription.
Jakob Skou Pedersen (University of Copenhagen) Genomic identification of structural RNAs using phylo-SCFGs
Abstract: RNA structures often evolve with characteristic substitution patterns that preserve base-pairs in spite of changes in primary sequence. With the advent of closely related full-length genomes, it has become possible to exploit this comparative signal for genomic identification of structural RNAs (1). Phylo-SCFGs (2) are attractive models for this problem since they can describe both RNA structure, using stochastic context-free grammars (SCFGs), and sequence evolution, using phylogenetic models. Using variations of classical algorithms, multiple alignments with any number sequences can be handled efficiently. EvoFold implements this approach and has been used to screen multiple-sequence genomic-alignments of both vertebrates and Drosopholids for structural RNAs (1,3). This has resulted in hundreds of high-confidence novel candidates of both ncRNAs and cis-regulatory structures. 1) Identification and Classification of Conserved RNA Secondary Structures in the Human Genome. Pedersen JS , Bejerano G, Siepel G, Rosenbloom K, Lindblad-Toh K, Lander ES, Kent J, Miller W, and Haussler D. PLoS Comput Biol. 2006 Apr;2(4):e33. 2) Using stochastic context free grammars and molecular evolution to predict RNA secondary structure. Knudsen B and Hein JJ. Bioinformatics. 1999; 15 (6): 446-454. 3) Discovery of functional elements in 12 Drosophila genomes using evolutionary signatures. Stark A, Lin MF, Kheradpour P, and Pedersen JS, et al. 2007 (in press).
Jakob Skou Pedersen (University of Copenhagen) Genomic identification of structural RNAs using phylo-SCFGs
Abstract: RNA structures often evolve with characteristic substitution patterns that preserve base-pairs in spite of changes in primary sequence. With the advent of closely related full-length genomes, it has become possible to exploit this comparative signal for genomic identification of structural RNAs (1). Phylo-SCFGs (2) are attractive models for this problem since they can describe both RNA structure, using stochastic context-free grammars (SCFGs), and sequence evolution, using phylogenetic models. Using variations of classical algorithms, multiple alignments with any number sequences can be handled efficiently. EvoFold implements this approach and has been used to screen multiple-sequence genomic-alignments of both vertebrates and Drosopholids for structural RNAs (1,3). This has resulted in hundreds of high-confidence novel candidates of both ncRNAs and cis-regulatory structures. 1) Identification and Classification of Conserved RNA Secondary Structures in the Human Genome. Pedersen JS , Bejerano G, Siepel G, Rosenbloom K, Lindblad-Toh K, Lander ES, Kent J, Miller W, and Haussler D. PLoS Comput Biol. 2006 Apr;2(4):e33. 2) Using stochastic context free grammars and molecular evolution to predict RNA secondary structure. Knudsen B and Hein JJ. Bioinformatics. 1999; 15 (6): 446-454. 3) Discovery of functional elements in 12 Drosophila genomes using evolutionary signatures. Stark A, Lin MF, Kheradpour P, and Pedersen JS, et al. 2007 (in press).
Niles A. Pierce (California Institute of Technology) Analysis and design of nucleic acid devices
Abstract: DNA and RNA are versatile construction materials. By appropriately designing the sequence of bases in each strand, synthetic nucleic acid systems can be programmed to self-assemble into complex structures that implement dynamic mechanical tasks. Motivated by the challenge of encoding arbitrary mechanical function into nucleic acid sequences, we are developing a suite of computational algorithms for analyzing the underlying free energy landscapes that control the behavior of a system. This talk will focus on new algorithms for predicting the equilibrium properties of an entire test tube of interacting nucleic acid strands. The utility of the approach will be demonstrated by elucidating the empirical behavior of hybridization chain reaction mechanisms that are under development with application to biosensing, transport, and therapeutics.
Jens Reeder (Universität Bielefeld) Locomotif: from graphical motif description to RNA motif search
Abstract: Motivated by the recent rise of interest in small regulatory RNAs, we present Locomotif - a new approach for locating RNA motifs that goes beyond the previous ones in three ways: (1) Motif search is based on efficient dynamic programming algorithms, incorporating the established thermodynamic model of RNA secondary structure formation. (2) Motifs are described graphically, using a Java-based editor, and search algorithms are derived from the graphics in a fully automatic way. The editor allows us to draw secondary structures, annotated with size and sequence information. They closely resemble the established, but informal way in which RNA motifs are communicated in the literature. Thus, the learning effort for Locomotif users is minimal. (3) Locomotif employs a client-server approach. Motifs are designed by the user locally. Search programs are generated and compiled on a bioinformatics server. They are made available both for execution on the server, and for download as C source code plus an appropriate make-file. Availability: Locomotif is available at http://bibiserv.techfak.uni-bielefeld.de/locomotif.
Jens Reeder (Universität Bielefeld) Locomotif: from graphical motif description to RNA motif search
Abstract: Motivated by the recent rise of interest in small regulatory RNAs, we present Locomotif - a new approach for locating RNA motifs that goes beyond the previous ones in three ways: (1) Motif search is based on efficient dynamic programming algorithms, incorporating the established thermodynamic model of RNA secondary structure formation. (2) Motifs are described graphically, using a Java-based editor, and search algorithms are derived from the graphics in a fully automatic way. The editor allows us to draw secondary structures, annotated with size and sequence information. They closely resemble the established, but informal way in which RNA motifs are communicated in the literature. Thus, the learning effort for Locomotif users is minimal. (3) Locomotif employs a client-server approach. Motifs are designed by the user locally. Search programs are generated and compiled on a bioinformatics server. They are made available both for execution on the server, and for download as C source code plus an appropriate make-file. Availability: Locomotif is available at http://bibiserv.techfak.uni-bielefeld.de/locomotif.
Yongwu Rong Feyman diagrams, RNA folding, and the transition polynomial
Abstract: Feynman diagrams were introduced by physicists. They arise naturally in mathematics (from knots and singular knots), and in molecular biology (from RNA folding). In particular, work of G. Vernizzi, H. Orland, and A. Zee has shown that the "genus" of Feynman diagrams plays an important role in the prediction of RNA structures. The transition polynomial for 4-regular graphs was defined by Jaeger to unify polynomials given by vertex reconfigurations similar to the skein relations of knots. It is closely related to the Kauffman bracket, Tutte polynomial, and the Penrose polynomial. We define a transition polynomial for Feynman diagrams and discuss its properties. In particular, we show that the genus of a Feynman diagram is encoded in the transition polynomial. This is joint work with Kerry Luse.
Yongwu Rong (George Washington University) Feynman diagrams, RNA folding, and the transition polynomial
Abstract: Feynman diagrams, RNA folding, and the transition polynomial
Walter Larry Ruzzo (University of Washington) Computational comparative genomics for discovery of cis-regulatory RNAs in bacteria
Abstract: Discovery of novel functional noncoding RNA is a multifaceted problem. Motif representation, inference and search are all important, as is incorporation of relevant biological knowledge. With careful attention to all of these, we have developed a comparative genomics "pipeline" for discovery of cis-regulatory RNA elements in bacteria. We represents motifs using covariance models (CMs), as in the Rfam database. Motif inference in unaligned sequences with extraneous flanking regions (i.e., local alignment) relies on CMfinder [2]. We apply it to intergenic regions upstream of homologous genes in different bacteria, since cis-regulatory elements often are found and conserved there [3]. We use Ravenna [1] for efficient, sensitive CM search to identify additional instances. This is critical since (i) a given RNA element often regulates multiple genes in a pathway, not just homologs, and (ii) more examples allow us to refine the model (also via CMfinder), in turn enabling further discoveries. This strategy recovers most known RNAs in Firmicutes [3]. More importantly, we discovered 6 new likely riboswitch families, most experimentally verified, plus over 20 other elements in a wide variety of bacteria [3,4]. Collectively, these RNAs are involved in diverse but individually specific cellular processes, such as ribosome biogenesis, molybdenum cofactor biosynthesis and the citric acid cycle. One of the more surprising finds is a widespread riboswitch that apparently regulates such disparate processes as natural competence in Vibrio cholerae and the use of metal ions as electron acceptors in Geobacter sulfurreducens. These candidate RNAs add to the growing list of RNA motifs involved in multiple cellular processes, and suggest that many additional RNAs remain to be discovered. Many computational challenges also remain. [1] Weinberg and Ruzzo. Sequence-based heuristics for faster annotation of non-coding RNA families. Bioinformatics, 2006, 22(1):35-39 [2] Yao, Weinberg and Ruzzo. CMfinder--A Covariance Model Based RNA Motif Finding Algorithm. Bioinformatics, 2006, 22(4): 445-452. [3] Yao, Barrick, Weinberg, Neph, Breaker, Tompa and Ruzzo. A Computational Pipeline for High Throughput Discovery of cis-Regulatory Noncoding RNA in Prokaryotes. PLoS Computational Biology. 3(7): e126, July 6, 2007. [4] Weinberg, Barrick, Yao, Roth, Kim, Gore, Wang, Lee, Block, Sudarsan, Neph, Tompa, Ruzzo and Breaker. Identification of 22 candidate structured RNAs in bacteria using the CMfinder comparative genomics pipeline. Nucl. Acids Res., July 2007 35: 4809-4819.
Jeffrey Saltzman (Merck & Co., Inc.) Modeling HCV antivirals
Abstract: Hepatitis C affects between 4 and 5 million people in the United States with nearly 75% suffering chronic infection. Further, current estimates indicate 170 million people are infected, worldwide, with the Hepatitis C virus (HCV). Pegylated interferon and ribavirin combination therapy is currently the standard of care for treatment of HCV but is considered less than optimal. For example, sustained virologic response (SVR) for genotype 1 virus is observed in pivotal clinical trials for only 54% of treated patients. SVR is defined as undetectable HCV RNA in plasma on the order of 6 months after cessation of treatment. This combination therapy is also associated with a high incidence of significant side effects suffered over a treatment interval of at least six months. Identification of better treatment alternatives is the goal of a vigorous antiviral program at Merck. In this talk we describe mathematical models explaining HCV infection and response to treatment. The model equations and resultant simulations were chiefly derived from the literature. Simulations have already provided support for the antiviral drug development program at Merck. Mathematical modeling gives evidence of target engagement and drug efficacy. Using accumulated simulation experience gained from basic research, preclinical data, and the literature, design optimizations for early clinical studies may be proposed. As internal modeling expertise is enhanced by experience gained in early development, this mathematical knowledge base may help optimize costly phase III trials by guiding key decisions such as dose selection and length of dosing.
Tamar Schlick (New York University) Designing Structured RNA Pools for RNA In Vitro Selection-- Hin Hark Gan, Namhee Kim and Tamar Schlick
Abstract: In vitro selection is a versatile experimental tool for discovering novel synthetic RNAs. However, most RNAs identified from random sequence pools are small and simple folding motifs. To significantly increase the probability of discovering novel RNAs, we develop an approach for engineering sequence pools that links RNA sequence space regions with corresponding structural distributions via a "mixing matrix" approach combined with a graph theory analysis of RNA 2D structure space. Our pool design method has been automated and made available through the web server RAGPOOLS (http://rubin2.biomath.nyu.edu). The RAGPOOLS serves as a guide to researchers who aim to synthesize RNA pools with desired properties and/or experiment in silico with various designs by our approach.
Tamar Schlick (New York University) Chromatin modeling and simulation: A tale of histone tails
Abstract: Eukaryotic chromatin is the fundamental protein/nucleic acid unit that stores the genetic material. Understanding how chromatin fibers fold and unfold is important for interpreting fundamental biological processes like DNA replication and transcription regulation. Using a mesoscopic model of oligonucleosome chains and tailored sampling protocols, we elucidate the energetics of oligonucleosome folding/unfolding and the role of each histone tail in regulating chromatin structure. Simulations reveal unfolding of oligonucleosome at low salt due to strong electrostatic repulsion between linker DNAs and moderate folding at higher salt due to a balance between the attractive inter-nucleosomal interactions (mediated by the histone tails) and repulsive interactions between the linker DNA. Furthermore, tail-mediated fiber/fiber interactions emerge as the oligonucleosome chain folds into more compact self interactions. Analyses of the tail positional distributions reveal a broad spread of tail positions consistent with the dynamic and flexible nature of the tails and the role of each tail in folding and energetics. These studies open the door to investigations of higher-order structures of compact chromatin and the biochemical modulation of chromatin structure by altered histone tail charge.

Of possible interest:

  • T Schlick, Molecular Modeling: An Interdisciplinary Guide, Springer-Verlag, New York, 2002;
  • D. Beard and T. Schlick, Structure 9: 105--114 (2001);
  • G. Arya et al., Biophys. J. 91: 133–150 (2006);
  • G. Arya & T. Schlick, Proc. Natl. Acad. Sci. 103: 16236–16241 (2006);
  • G. Arya & T. Schlick, J. Chem. Phys. 126: 044107 (2007).
Bruce A. Shapiro (National Cancer Institute) Computational approaches to RNA nanodesign
Abstract: We have developed a number of computational tools that permit a user to design RNA based nano-particles with various functionalities. One of these tools is a newly developed relational database, RNAJunction, which contains structural and sequence information for all known RNA n-way junctions and kissing loop interactions. The database also contains the results from applying molecular mechanics and structural clustering techniques to the motifs. The database of motifs can be searched in a variety of ways and provide a source for further analysis and RNA nano building blocks. Another computational tool, NanoTiler, permits a user to interactively and automatically construct user specified RNA-based nano-scale shapes. The combination of the RNAJunction database, NanoTiler and other computational tools allows the rapid prototyping of designed RNA shapes. We discuss some of the principles involved in these design tools and show how the RNA nanodesign process can be accomplished with the use of these methodologies.
Tobin R. Sosnick (University of Chicago) Imaging RNA structures and folding intermediates using electron cryo-microscopy
Abstract: We investigate the applicability of electron cryomicroscopy (cryo-EM) with single particle reconstruction in the RNA structural studies as small as 154 residues (~50 kD). This size is at least two-fold smaller than the generally conceived limits for single-particle image reconstruction by cryo-EM of macromolecules. For the Specificity and Catalytic domain of bacterial RNase P RNA, single-particle reconstruction of the native structures exhibits good agreement with their respective crystal structures. For the major thermodynamic folding intermediate of the /B. subtilis/ specificity domain, the single-particle reconstruction has considerable similarity to the previously proposed structural models of this intermediate. These results indicate that cryo-EM can directly image conformations of relatively small RNA molecules in different structural and functional states.
Andrew Stein (University of Minnesota) Optimization in mathematical biology
Abstract: A common part of many mathematical biology problems involves optimization. Two major areas where it's important is in finding the optimal set of parameters to fit a model and in minimizing the energy of a biophysical system. In this talk, the speaker will present an overview of a few simple optimization algorithms, including the Nelder-Mead Simplex Method, Newton-Raphson, Conjugate Gradient and provide examples of their use in his own research on tumor invasion and collagen micromechanics. If time permits, the speaker will also touch on global optimization methods. A note to the audience: The author is not an expert on optimization techniques and one purpose of this talk is to improve his own understanding of the field, as well the audience's. Comments and discussion at this seminar will be welcome.
Michael Stich (Instituto Nacional de Tecnica Aeroespacial) Collective properties of evolving populations of RNA molecules
Abstract: RNA molecules, through their dual appearance as sequence and structure, represent a suitable model to study evolutionary properties of quasispecies. The essential ingredient in this model is the differentiation between genotype (molecular sequences which are affected by mutation) and phenotype (molecular structure, affected by selection). This framework allows a quantitative analysis of organizational properties of quasispecies as they adapt to different environments, such as their robustness, the effect of the degeneration of the sequence space, or the adaptation under different mutation rates and the error threshold associated.
Bernd Sturmfels (University of California) Algebra, statistics, computation, and biology
Abstract: Recently, algebraic methods have been developed to unify and advance a variety of techniques of statistical analysis, providing new and improved tools for computational biology. Professor Sturmfels, one of the founders of the new field of algebraic statistics, will introduce the subject and describe its emerging applications to genome science and developmental biology. He will be assisted by a fictional character named DiaNA who plays hopscotch and rolls tetrahedral dice with faces labeled "A," "C," "G," and "T."
Chris Thachuk (University of British Columbia) On the design of oligos for gene synthesis
Abstract: Methods for reliable synthesis of long genes offer great promise for protein synthesis via expression of synthetic genes, with applications to improved analysis of protein structure and function, as well as engineering of novel proteins. Current technologies for gene synthesis use computational methods for design of short oligos, which can then be reliably synthesized and assembled into the desired target gene. For collision-oblivious oligo design -- when mishybridizations between oligos are ignored -- we give a simple and efficient dynamic programming algorithm. We conjecture that the collision-aware oligo design problem is NP-hard and provide evidence that mishybridizations between oligos occur infrequently in the designs from the collision-oblivious algorithm. We extend our dynamic programing algorithm to achieve collision-aware oligo design, when the target gene can be partitioned into independently-assembled short segments. We evaluate our methods on a large biological gene set.
Devarajan Thirumalai (University of Maryland) Exploring the energy landscape of RNA
Abstract: Recent single molecule experiments and high-resolution temperature jump experiments show that the energy landscape of RNA is rugged. As a result, even the formation of a hairpin, exhibits all the signatures of folding (multiple pathways and complex kinetics) usually associated with self-assembly of ribozymes. I will describe the kinetics of hairpin formation, initiated by both temperature and force quench, using computations. The profound differences between the two methods will be illustrated in terms of the pathways to the native state. Analogies to folding of ribozymes will also be given.
Elfar Torarinsson (University of Copenhagen) Comparative genomics beyond sequence based alignments: RNA structures in the ENCODE regions
Abstract: Joint work with Z. Yao, E. D. Wiklund, J. B. Bramsen , C. Hansen, J. Kjems, N. Tommerup, W. L. Ruzzo, and J. Gorodkin. Recent computational scans for non-coding RNAs (ncRNAs) in multiple organisms have relied on existing multiple sequence alignments. However, as sequence similarity drops, a key signal of RNA structure - frequent compensating base changes - is increasingly likely to cause sequence-based alignment methods to misalign, or even refuse to align, homologous ncRNAs, consequently obscuring that structural signal. We have used CMfinder [1], a structure-oriented local alignment tool, to search vertebrate multiple alignments in the ENCODE regions. In agreement with other studies [2], we find a large number of potential RNA structures in the ENCODE regions. We report 6,587 candidates with an estimated false positive rate of 50%. More intriguingly, many of these candidates may be better represented by alignments taking the RNA secondary structure into account than those based on primary sequence alone, often quite dramatically. For example, approximately one quarter of these 6,587 candidates show revisions in more than 50% of their aligned positions. Furthermore, our results are strongly complementary to those discovered by sequence-alignment-based approaches—84% of our candidates are not covered by Washietl et al.[2], increasing the number of ncRNA candidates in the ENCODE region by 32%. In a group of eleven ncRNA candidates that were tested by RT-PCR, 10 were confirmed to be present as RNA transcripts in human tissue. Our results broadly suggest caution in any analysis relying on multiple sequence alignments in less well-conserved regions, clearly support growing appreciation for the biological significance of ncRNAs, and strongly argue for considering RNA structure directly in any searches for these elements. 1. Yao, Z., Weinberg, Z. and Ruzzo, W.L. 2006. CMfinder - A Covariance Model Based RNA Motif Finding Algorithm. Bioinformatics 22: 445-452. 2. Washietl, S., Pedersen, J.S., Korbel, J.O., Gruber, A.R., Hackermuller, J., Hertel, J., Lindemeyer, M., Reiche, K., Stocsits, C., Tanzer, A., et al. 2007. Structured RNAs in the ENCODE Selected Regions of the Human Genome. Genome Research 17: 852-864.
Mariel Vazquez (San Francisco State University) Modeling DNA unknotting by type II topoisomerases
Abstract: Type II topoisomerases simplify DNA knots and links efficiently by performing strand-passage on DNA strands. Experimental studies have shown that these enzymes simplify the topology of DNA below thermodynamical equilibrium, however the key behind their efficiency is yet to be revealed. Motivated by these experimental observations, we study random transitions of knotted polygonal chains of fixed length. We use Monte Carlo computer simulations and computational knot theory methods to model strand-passage, with and without topological biases, on the knotted chains and to identify the resulting knotting probabilities. This project is funded by NIH MBRS SCORE grant S06 GM052588.
Jerome Waldispuhl (Massachusetts Institute of Technology) Efficient algorithms for probing the RNA mutation landscape and prediction of deleterious mutations
Abstract: We develop an efficient algorithm to compute, for a given RNA sequence and simultaneously for each k, the minimum free energy structure MFE_k and the Boltzmann partition function Z_k over all secondary structures of all k-point mutants of the given sequence. Using the partition function, we rigorously sample from the ensemble of low energy k-point mutants in order to explore the mutation landscape. Our algorithm, named RNAmutants, allows us to investigate deleterious mutations (mutations that radically modify secondary structure) in the Hepatitis C virus cis-acting replication (HCV CAR) element and the hairpin of human immunodeficiency virus trans-activation response (HIV-1 TAR) element. More generally, using RNAmutants, we study the resiliance of an RNA molecule to pointwise mutations. By computing the mutation profile of a sequence, a novel graphical representation of the mutational tendency of nucleotide positions, we analyze the deleterious nature of mutating specific nucleotide positions or groups of positions. In particular, we show qualitative agreement between published HIV experimental mutagenesis studies and our analysis of deleterious mutations using RNAmutants. Our work predicts other deleterious mutations, which could be verified experimentally. Work in collaboration with P. Clote, B. Berger and S. Devadas.
Stefan Washietl (Universität Wien) Improved RNA gene predictions through dinucleotide controlled randomization of multiple sequence alignments
Abstract: Tanja Gesell (1) & Stefan Washietl (2,3) 1. Center for Integrative Bioinformatics, Max Perutz Laboratories, Vienna
2. EMBL-European Bioinformatics Institute, Hinxton United Kingdom
3. Department of Theoretical Chemistry, University of Vienna, Austria Most noncoding RNA gene prediction programs are based on the detection of conserved RNA secondary structures in multiple alignments [1]. Although this approach seems to be the most promising, the main problem of current algorithms is the large number of false positive predictions, in particular in large vertebrate genomes. As the available algorithms assume a mononucleotide background model, a major source of erroneously predicted RNA structures is the biased dinucleotide content found e.g. in vertebrate genomes. While there are well known algorithms for randomization of single sequences preserving dinucleotide content, no algorithms exist for multiple alignments. We present a novel algorithm addressing this problem. Our approach which involves in silico evolution along a phylogenetic tree has two key features: (i) We make use of a new evolutionary model that considers site specific and overlapping dependencies [2]. This enables us to simulate alignments with given dinucleotide (or higher order nucleotide) content. (ii) The model includes site-specific rate factors that preserve critical conservation patterns of the original alignments. We developed a time-efficient distance based approximation method to estimate a tree under this complex model which is used as guide for simulating new alignments. Based on this improved null model we have implemented a noncoding RNA gene prediction algorithm called SISSIz, that builds upon the RNAalifold and AlifoldZ programs. The new dinucleotide based program shows significantly improved accuracy over its mononucleotide counterparts on a vertebrate test set. [1] Washietl S., Hofacker I.L., Lukasser M., Huettenhofer A., Stadler P.F. Mapping of conserved RNA secondary structures predicts thousands of functional noncoding RNAs in the human genome. Nat. Biotechnol. (2005), 23:1383-90. [2] Gesell T., von Haeseler A. In silico sequence evolution with site-specific interactions along phylogenetic trees. Bioinformatics. (2006), 22:716-22
Eric Westhof (Université de Strasbourg I (Louis Pasteur)) Architecture and reactivity of RNA
Abstract: RNA architecture results from the hierarchical assembly of preformed double-stranded helices defined by Watson-Crick base pairs and RNA modules maintained by non-Watson-Crick base pairs. Surprisingly, the most common RNA-RNA interaction motif, the A-minor motif, is also the least specific in its local requirements. A-minor motifs are mediated by adenines binding into the shallow/minor groove of any combination of stacked and helical Watson-Crick base pairs. Thus, A-minor motifs are mutationally robust and can accommodate many combinations of neutral mutations. This complicates the search of functional RNAs in genomes and dilutes the links between RNA structure and evolution. The bacterial ribosomal decoding A site exploits this lack of local atomic specificity. There, the adenines A1492 and A1493 of the A site are seen either tucked in within the internal loop or bulging out and poised for interaction. This dynamic equilibrium contributes to the decoding process during recognition of the codon:anticodon Watson-Crick base pairings. In contrast, for RNA folding, where specificity is a requirement, global, positional and orientational, constraints on the native fold must occur upstream in the folding process. Critical parameters are the lengths of the helices, the co-axiality of the helical stacks, and the structure adopted at the junctions of helices. The molecular neutrality present in the local interactions is thus partially compensated by these global topological criteria, much less accessible to sequence analysis since they are attached to the three-dimensional architecture. The search for functional RNAs in genomes is thereby complexified through this dilution of the direct links between sequences and structures. The simultaneous treatment of 3D structures, structural alignments, and annotations of the interactions should allow hopefully to derive some rules of molecular evolution in structured RNAs. Lescoute, A. and Westhof, E. (2006) The interaction networks of structured RNAs. Nucleic Acids Res 34, 6587. Hammann, C. and Westhof, E. (2007) Searching genomes for ribozymes and riboswitches. Genome Biology 8, 210.
Sarah Woodson (Johns Hopkins University) How RNA tells right from wrong: base pairs, tertiary interactions, and counterions in RNA folding
Abstract: RNAs must self-assemble into unique three-dimensional structures in the cell, yet how RNA molecules find their native structure in a short time is not well understood. This is a challenging problem, because RNA secondary structures are thermodynamically stable but not uniquely specified by the sequence, while tertiary interactions are specific but not very stable. Consequently, many RNAs become trapped in metastable, non-native intermediates. Recent footprinting and SAXS experiments on a bacterial ribozyme show that tertiary interactions make helix assembly more specific during the initial collapse transition. Specific collapse increases the flux through folding pathways that lead directly to the native structure. The stability of the folded RNA also depends on the charge density of the counterions; small multivalent counterions stabilize the RNA more than large monovalent ions. In low charge density counterions, the transition state ensemble becomes broader, accelerating the search for the native structure.
Yurong Xin (New York University) Estimating the fraction of non-coding RNAs in mammalian transcriptomes
Abstract: Recent studies of mammalian transcriptomes have identified numerous RNA transcripts that do not code for proteins; their identity, however, is largely unknown. Here we explore an approach based on sequence randomness patterns to discern different RNA classes. The relative z-score we use helps identify the known ncRNA class from the genome, intergene, and intron classes. This leads us to a fractional ncRNA measure of putative ncRNA datasets which we model as a mixture of genuine ncRNAs and other transcripts derived from genomic, intergenic and intronic sequences. We use this model to analyze six representative datasets, identified by the FANTOM3 project and two computational approaches based on comparative analysis (RNAz and EvoFold). Our analysis suggests fewer ncRNAs than estimated by DNA sequencing and comparative analysis, but the verity of our approach and its prediction requires more extensive experimental RNA data.
Craig L. Zirbel (Bowling Green State University) Using RNA 3D structure data in SCFG/MRF models to do sequence alignment and motif inference
Abstract: RNA 3D structure files contain essentially complete information about the interactions that form the 3D structure of an RNA molecule for a given organism. Homologous molecules in other organisms will have very similar 3D structures, but we expect to see sequence variability due to structurally neutral base substitutions, insertions, and deletions, among other things. RNA databases have far more RNA sequences than RNA 3D structures, and this will always be the case. We wish to use the 3D structure data to make inferences about the 3D structure of homologous molecules on the basis of their sequences. We think of homologous RNA sequences as being random variants of the molecule for which we have a 3D structure. If the probabilistic model of this variation is simple enough, we can use it to align the sequences to the 3D structure, and thus infer the structural role of each base in the sequence. Stochastic context free grammars (SCFGs) can account for the nested Watson-Crick basepairs prevalent in RNA, and by choosing appropriate basepair substitution probabilities, they can be used to model structurally neutral basepair substitutions for non-Watson-Crick basepairs as well. We use an SCFG formalism enhanced by production rules based on Markov Random Fields (MRF). This allows us to model base triples such as are found in sarcin motifs, and local crossing interactions such as are found in kink turn internal loops. However, as usual with SCFG, it does not allow us to model longer-range pseudoknots. The SCFG/MRF model can be used for two purposes: First, to make RNA multiple sequence alignments based on the 3D structure of one molecule, without reference to a hand-curated seed alignment. Second, to infer the 3D structure of small motifs such as internal loops from their sequences.
Michael Zuker (Rensselaer Polytechnic Institute) Computational methods for RNA secondary structure determination
Abstract: The talk will begin with a definition of RNA secondary structure, including three different ways to display these structures. Two distinct approaches will be presented for determining secondary structure from sequence data. The comparative method requires a multiple sequence alignment of a collection of homologous RNA sequences. It uses phylogeny to determine common, conserved base pairs that are more likely to be the result of evolution than to exist by chance. On the other hand, recursive algorithms may be used on single sequences to compute minimum free energy structures, partition functions and other biophysical quantities. These algorithms ignore evolution and use empirically derived energy parameters based on physical chemistry. Examples will be given for both methods.
Visitors in Residence
Mirela Andronescu University of British Columbia 10/28/2007 - 11/2/2007
Douglas N. Arnold University of Minnesota 7/15/2001 - 6/30/2008
F. Javier Arsuaga San Francisco State University 9/3/2007 - 12/21/2007
Rolf Backofen Albert-Ludwigs-Universität Freiburg 10/27/2007 - 11/4/2007
Alex Bateman Wellcome Trust Sanger Institute 10/28/2007 - 11/2/2007
Daniel J. Bates University of Minnesota 9/1/2006 - 8/31/2008
Peter W. Bates Michigan State University 9/1/2007 - 12/31/2007
John Baxter University of Minnesota 8/1/2007 - 7/30/2009
Jean Bergeron Laval University 10/31/2007 - 11/2/2007
Yermal Sujeet Bhat University of Minnesota 9/1/2006 - 8/31/2008
Eckart Bindewald SAIC-Frederick, Inc. 10/28/2007 - 11/2/2007
Michael Bon Commissariat à l'Énergie Atomique (CEA)-Centre d'Études Nucléaires de Saclay (CENS) 10/28/2007 - 11/2/2007
Michael P. Brenner Harvard University 10/7/2007 - 10/8/2007
Steven Brenner University of California 10/28/2007 - 10/31/2007
Anke Busch Albert-Ludwigs-Universität Freiburg 10/26/2007 - 11/4/2007
Maria-Carme T. Calderer University of Minnesota 10/29/2007 - 11/2/2007
Hannah Callender University of Minnesota 9/1/2007 - 8/31/2009
David Case Scripps Research Institute 10/28/2007 - 11/2/2007
Bernard Chazelle Princeton University 10/6/2007 - 10/8/2007
Shi-Jie Chen University of Missouri 9/3/2007 - 11/3/2007
Arkadiusz Chworos University of California 10/28/2007 - 11/2/2007
Peter Clote Boston College 10/28/2007 - 11/2/2007
Anne E. Condon University of Wisconsin 10/28/2007 - 11/2/2007
L. Pamela Cook University of Delaware 10/6/2007 - 10/7/2007
Ludovica Cecilia Cotta-Ramusino University of Minnesota 10/1/2007 - 8/30/2009
Irma Cruz-White Chipola College 10/28/2007 - 11/2/2007
Isabel K. Darcy University of Iowa 9/1/2007 - 1/19/2008
Kaushik Dayal University of Minnesota 10/29/2007 - 11/2/2007
Maciej Dlugosz University of Warsaw 10/27/2007 - 11/3/2007
Charles Doering University of Michigan 10/7/2007 - 10/8/2007
Olivier Dubois University of Minnesota 9/3/2007 - 8/31/2009
Mauricio Esguerra Rutgers University 10/28/2007 - 11/2/2007
Marcia O. Fenley Florida State University 10/23/2007 - 11/2/2007
Jes Frellsen University of Copenhagen 10/27/2007 - 11/3/2007
Hin Hark Gan New York University 10/28/2007 - 11/2/2007
Paul Gardner Wellcome Trust Sanger Institute 10/28/2007 - 11/3/2007
Daniel Gautheret Université de Paris XI (Paris-Sud) 10/29/2007 - 11/2/2007
Zhumur Ghosh Indian Association for the Cultivation of Science (IACS) 10/28/2007 - 11/3/2007
George Madalin Giambasu University of Minnesota 10/29/2007 - 11/2/2007
Anna Gilbert University of Michigan 10/6/2007 - 10/8/2007
Jan Gorodkin University of Copenhagen 10/28/2007 - 11/2/2007
Jason E. Gower University of Minnesota 9/1/2006 - 8/31/2008
Christine E. Heitsch Georgia Institute of Technology 9/4/2007 - 11/3/2007
Milena Hering University of Minnesota 9/1/2006 - 8/31/2008
Peter Hinow University of Minnesota 9/1/2007 - 8/31/2009
Ivo L. Hofacker Universität Wien 10/28/2007 - 11/2/2007
Mihaela Iftime Boston University 10/28/2007 - 11/1/2007
Tan Inoue Kyoto University 10/28/2007 - 11/3/2007
Hosna Jabbari University of British Columbia 10/28/2007 - 11/2/2007
Luc Jaeger University of California 10/28/2007 - 11/2/2007
Richard D. James University of Minnesota 9/4/2007 - 6/30/2008
Tiefeng Jiang University of Minnesota 9/1/2007 - 6/30/2008
Hans Kaper National Science Foundation 10/6/2007 - 10/8/2007
Wojciech (Voytek) Kasprzak National Cancer Institute 10/28/2007 - 11/2/2007
Markos A. Katsoulakis University of Massachusetts 10/27/2007 - 10/30/2007
Namhee Kim New York University 10/27/2007 - 11/2/2007
Soojeong Kim University of Iowa 8/30/2007 - 1/20/2008
Rob Knight University of Colorado 10/28/2007 - 11/2/2007
Debra Knisley East Tennessee State University 8/17/2007 - 6/1/2008
Robert V. Kohn New York University 10/6/2007 - 10/8/2007
Christian Laing New York University 10/28/2007 - 11/2/2007
Chang Hyeong Lee Worcester Polytechnic Institute 10/14/2007 - 1/4/2008
Hyunju Lee University of Minnesota 10/29/2007 - 11/2/2007
Neocles B. Leontis Bowling Green State University 10/28/2007 - 11/2/2007
Anton Leykin University of Minnesota 8/16/2006 - 8/15/2008
David M.J. Lilley University of Dundee 10/28/2007 - 11/2/2007
Stinus Lindgreen University of Copenhagen 10/27/2007 - 11/4/2007
Andy Lorenz Boston College 10/28/2007 - 11/2/2007
Roger Lui Worcester Polytechnic Institute 9/1/2007 - 6/30/2008
Laura Lurati University of Minnesota 9/1/2006 - 8/31/2008
Maciej Maciejczyk Cornell University 10/28/2007 - 11/4/2007
Francois Major University of Montreal 10/28/2007 - 11/2/2007
Bibekanand Mallick Indian Association for the Cultivation of Science (IACS) 10/28/2007 - 11/2/2007
David H. Mathews University of Rochester 10/27/2007 - 11/2/2007
Roderick Melnik Wilfrid Laurier University 10/28/2007 - 11/2/2007
Ezra Miller University of Minnesota 9/1/2007 - 6/30/2008
Kenneth C. Millett University of California 9/13/2007 - 10/17/2007
Maria Giovanna Mora International School for Advanced Studies (SISSA/ISAS) 9/1/2007 - 12/31/2007
Alejandro Morales Valencia University of Guadalajara 10/27/2007 - 11/4/2007
Junalyn Navarra-Madsen Texas Woman's University 10/27/2007 - 11/3/2007
Timothy Newman Arizona State University 9/1/2007 - 6/30/2008
Olalla Nieto Faza University of Minnesota 10/29/2007 - 11/2/2007
Benjamin Nill Freie Universität Berlin 9/30/2007 - 10/5/2007
Asamoah Nkwanta Morgan State University 10/28/2007 - 11/2/2007
Ruth Nussinov National Cancer Institute 10/28/2007 - 10/31/2007
Duane Nykamp University of Minnesota 9/1/2007 - 6/30/2008
Wilma K. Olson Rutgers University 10/29/2007 - 11/2/2007
Henri Orland Commissariat à l'Énergie Atomique Saclay (CEA) 10/28/2007 - 11/2/2007
Hans G. Othmer University of Minnesota 9/1/2007 - 6/30/2008
Tao Pan University of Chicago 10/28/2007 - 11/2/2007
Jakob Skou Pedersen University of Copenhagen 10/28/2007 - 11/2/2007
Arlie O. Petters Duke University 10/6/2007 - 10/8/2007
Niles A. Pierce California Institute of Technology 10/28/2007 - 10/30/2007
Petr Plechac University of Tennessee 10/27/2007 - 10/30/2007
Craig T. Poling Lockheed Martin 10/7/2007 - 10/8/2007
Yann Ponty Boston College 10/28/2007 - 11/2/2007
Eric Rawdon University of St. Thomas 10/29/2007 - 11/2/2007
Jens Reeder Universität Bielefeld 10/28/2007 - 11/2/2007
Donald Richards Pennsylvania State University 10/6/2007 - 10/8/2007
Yongwu Rong George Washington University 10/30/2007 - 11/2/2007
Laureen Ross McCalib University of Minnesota 10/7/2007 - 10/8/2007
Ioulia Rouzina University of Minnesota 10/29/2007 - 11/2/2007
Walter Larry Ruzzo University of Washington 10/28/2007 - 11/2/2007
Jeffrey Saltzman Merck & Co., Inc. 10/25/2007 - 10/27/2007
Fadil Santosa University of Minnesota 10/6/2007 - 10/12/2007
Tamar Schlick New York University 10/25/2007 - 11/2/2007
Deena Schmidt University of Minnesota 9/1/2007 - 8/31/2009
Chehrzad Shakiban University of Minnesota 9/1/2006 - 8/31/2008
Bruce A. Shapiro National Cancer Institute 10/28/2007 - 11/2/2007
Tobin R. Sosnick University of Chicago 10/28/2007 - 11/2/2007
Andrew Stein University of Minnesota 9/1/2007 - 8/31/2009
Michael Stich Instituto Nacional de Tecnica Aeroespacial 10/27/2007 - 11/4/2007
Andrew M. Stuart University of Warwick 10/29/2007 - 11/2/2007
Bernd Sturmfels University of California 10/9/2007 - 10/10/2007
Bernd Sturmfels University of California 10/7/2007 - 10/8/2007
De Witt L. Sumners Florida State University 10/24/2007 - 11/8/2007
Vladimir Sverak University of Minnesota 9/1/2007 - 6/30/2008
David Swigon University of Pittsburgh 9/4/2007 - 12/14/2007
Berhanu Tameru Tuskegee University 10/28/2007 - 11/2/2007
Chris Thachuk University of British Columbia 10/28/2007 - 11/2/2007
Devarajan Thirumalai University of Maryland 10/28/2007 - 11/2/2007
Madhan Tirumalai University of Houston 10/28/2007 - 11/2/2007
Elfar Torarinsson University of Copenhagen 10/28/2007 - 11/4/2007
Erkan Tüzel University of Minnesota 9/1/2007 - 8/31/2009
Tanya Vassilevska Lawrence Livermore National Laboratory 10/28/2007 - 11/2/2007
Mariel Vazquez San Francisco State University 9/3/2007 - 12/21/2007
Dionisios G. Vlachos University of Delaware 10/28/2007 - 10/30/2007
Jerome Waldispuhl Massachusetts Institute of Technology 10/28/2007 - 11/3/2007
Haiyan Wang Arizona State University 10/28/2007 - 11/1/2007
Zhian Wang University of Minnesota 9/1/2007 - 8/31/2009
Stefan Washietl Universität Wien 10/28/2007 - 11/2/2007
Annika Wedemeier German Cancer Research Center 9/15/2007 - 10/12/2007
Eric Westhof Université de Strasbourg I (Louis Pasteur) 10/28/2007 - 11/2/2007
Sebastian Will Albert-Ludwigs-Universität Freiburg 10/26/2007 - 11/4/2007
Robert Williams US Air Force Research Laboratory 10/17/2007 - 10/17/2007
Sarah Woodson Johns Hopkins University 10/31/2007 - 11/2/2007
Zhijun Wu Iowa State University 9/4/2007 - 6/1/2008
Yurong Xin New York University 10/27/2007 - 11/2/2007
Yi Xing University of Iowa 10/28/2007 - 11/2/2007
Yanji Xu Minnesota Supercomputing Institute 10/29/2007 - 11/2/2007
Masoud Yari Indiana University 10/29/2007 - 11/2/2007
Darrin M. York University of Minnesota 10/28/2007 - 11/2/2007
Hongchao Zhang University of Minnesota 9/1/2006 - 8/31/2008
Zhichuan Zhang University of British Columbia 10/28/2007 - 11/2/2007
Craig L. Zirbel Bowling Green State University 10/28/2007 - 11/2/2007
Michael Zuker Rensselaer Polytechnic Institute 10/28/2007 - 11/2/2007
Legend: Postdoc or Industrial Postdoc Long-term Visitor

IMA Affiliates:
3M, Boeing, Carnegie Mellon University, Corning, ExxonMobil, Ford, General Electric, General Motors, Georgia Institute of Technology, Honeywell, IBM, Indiana University, Iowa State University, Johnson & Johnson, Kent State University, Lawrence Livermore National Laboratory, Lockheed Martin, Los Alamos National Laboratory, Medtronic, Michigan State University, Michigan Technological University, Microsoft Research, Mississippi State University, Motorola, Northern Illinois University, Ohio State University, Pennsylvania State University, Purdue University, Rice University, Rutgers University, Sandia National Laboratories, Schlumberger-Doll, Schlumberger-Doll Research, Seoul National University, Siemens, Telcordia, Texas A & M University, University of Central Florida, University of Chicago, University of Cincinnati, University of Delaware, University of Houston, University of Illinois at Urbana-Champaign, University of Iowa, University of Kentucky, University of Maryland, University of Michigan, University of Minnesota, University of Notre Dame, University of Pittsburgh, University of Tennessee, University of Texas, University of Wisconsin, University of Wyoming, US Air Force Research Laboratory, Wayne State University, Worcester Polytechnic Institute