search

ash211 / gtfold

GTfold is a fast, scalable multicore code for predicting RNA secondary structure and is one to two orders of magnitude faster than the de facto standard programs and achieves comparable accuracy of prediction.
http://gtfold.sourceforge.net/
5 stars 18 forks source link

readme

GTfold: A Scalable Multicore Code for RNA Secondary Structure Prediction

(C) 2007-2010, David A. Bader Amrita Mathuriya College of Computing, Georgia Institute of Technology Christine E. Heitsch School of Mathematics, Georgia Institute of Technology Stephen C. Harvey School of Biology, Georgia Institute of Technology Sainath Mallidi College of Computing, Georgia Institute of Technology

The prediction of the correct secondary structures of large RNAs is one of the unsolved challenges of computational molecular biology. Among the major obstacles is the fact that accurate calculations scale as O(n^4), so the computational requirements become prohibitive as the length increases. Existing folding programs implement heuristics and approximations to overcome these limitations. GTfold, a parallel multicore and scalable program, is one to two orders of magnitude faster than the de facto standard programs and achieves comparable accuracy of prediction. Development of GTfold opens up a new path for the algorithmic improvements and application of an improved thermodynamic model to increase the prediction accuracy.

This work was supported in part by National Institutes of Health grant NIH NIGMS R01 GM083621, and by NSF Grants CNS-0614915 and DBI-04-20513. Additionally, the research of Christine E. Heitsch, Ph.D., is supported in part by a Career Award at the Scientific Interface (CASI) from the Burroughs Wellcome Fund (BWF).

Additional contributors to GTfold's implementation include Georgia Tech visiting students Gregory Nou, Sonny Hernandez, and Manoj Soni.

INSTALLATION - Please execute the following commands

gtar zxvf gtfold-version.tar.gz cd gtfold-version

1) Installing with the root user b) ./configure c) make d) make install

2) Installing in a local directory a) ./configure --prefix=(Full path of the directory where gtfold will be installed.) b) make c) make install

Also, you can specify C and C++ compiler using the CC and CXX variables. For example - ./configure CC=gcc CXX=g++ --prefix=/home/username/gtfold

Understanding Source Code - 1) Please use a program called CSCOPE to get understanding of the source code. It may be preinstalled in your Linux system or get it from http://cscope.sourceforge.net/ 2) Please refer Chapter 3 of Mirela Stefania Andronescu's Masters thesis understanding of calculating the energy of various types of loops. http://www.cs.ubc.ca/grads/resources/thesis/Nov03/Mirela_Andronescu.pdf 3) You can look Chapter 4 to understand the thermodynamic recurssion formulas. Also, pesudocode given in the paper should be helpful. 4) To understand traceback, you may want to look at the following publication - Biopolymers, Volume 49 Issue 2, Pages 145 - 165 Complete suboptimal folding of RNA and the stability of secondary structures Stefan Wuchty, Walter Fontana, Ivo L. Hofacker, Peter Schuster 5) Please look at the paper here - A. Mathuriya, D.A. Bader, C.E. Heitsch, and S.C. Harvey, ``GTfold: A Scalable Multicore Code for RNA Secondary Structure Prediction,''
24th Annual ACM Symposium on Applied Computing (SAC), Computational Sciences Track, Honolulu, HI, March 8-12, 2008.

Version History:

1.0. 16 August 2008 Base Version 1.1. 17 August 2008 Converted DP tables from static to dynamic allocation 1.2. 18 August 2008 Fixed a bug with chPair rnap2/continue introduced in gtfold-1.1 Changed to 0-based indexing for dangle, stack, tskth, tskti, iloop22, iloop21, iloop11. 1.3. 20 August 2008 Changed RNA mapping to 0 based 1.4 23 August 2008 Flattened arrays stack, tskth, tskti. 1.5 29 August 2008 Added OpenMP parallelization 1.6 4 October 2008 Changed calcW function in algorithms.c to have Wim1 = MIN (W[i-1],0) in place of Wim1 = W[i-1]. Initialized W[i] values to INFINITY_ in function initTables in place of 0. loader.cc file is modified to open the data files in read mode instead of read/write mode. The sequence file is opened in read mode instead of write mode in file main.cc in main function. The version directly produces ct files from the sequence in the working directory. 1.7 22 October 2008 Corrected a bug in file loader.cc, function - initInt22Values 1.8 03 January 2009 (1) Adding functionality of Internal Loop Speedup Algorithm ( ILSA )
(2) Corrected one invalid memory access bug, in function calcWM() of file algorithms.c (3) MFE will be printed in correct floating point format - change in file main.c (4) Converted constant 30, 31, 32 of function calcVBI() in terms of MAXLOOP constant (5) Modified gtfold command line options. Run "gtfold -help" to see the usage. (6) Now, multConst array contains the multiloop penalty constants. Functions for multicloop energy calculations are modified to read the penalties from this array. 1.9 15 January 2009 Documented GTfold with various comments.
1.10 22 January 2009 More documentation 1.11 30 Spetember 2009 Added Better support for input sequence files (removing intermediate white spaces) 1.12 5 Mar 2010 Added support for constraints, forcing, prohibitng a base pair, specifying single stranded regions. Modified gtfold command line options. Run "gtfold -help". 1.13 25 Mar 2010 Added a toggle switch for preventing isolated base pairs Handling constraints has been changed drastically and it produces much better output now. Modified gtfold command line options. Run "gtfold -help". 1.14 30 Mar 2010 Added support to change data directory on fly. Support for generating binaries which reads the data files from its subdir data/ More robust handling of input file sequences Changed the output print sequence of GTfold, now also prints in Vienna format Modified the command line options. Run "gtfold -help" to see the usage. Changed the output format of traceback to print loop energies instead of running energy. Changed the initialized values of W[i] to 0 to play nice with the constraints. 1.15 18 May 2010 Fixed a nasty bug with reading constraints by increasing the length of string to 100 that reads the lines from sontraints file. 1.16 20 May 2010 Fixed another nasty bug with prohibited constraints handling. Earlier not used to update the WM array if a base pair is prohibited which is causing GTfold to trace partial structure. Some documentation added.