Team:USTC Software/Notebook

Overall Construction (Wei Pan, Bo Ding)

Week 2

“What can we do for the experimentalists?” I saw my team members. It’s the third time for me to take part in iGEM. Interestingly, in this time I am the only person with the background of biology, while nearly everyone’s background is solo. It is really a mixture. I hope there will be chemical reactions among us.

“Design”, the key word floated in everyone’s mind during the discussion. The network inference from the microarray expression profile, for protein concentration profiles, from promoter sequence data, etc, can be done in the DREAM (Database for Reverse Engineering Analysis and Method) project. Meanwhile, the network design can also be done to the required profile. So we will start from here!

So the aims for our project: 1.Present a structure frame for network design 2.Redesign the known networks to test the strategy 3.In silico, design a network with a interesting behavior 4.Realize the network in the cell

And the problems are really hard for our projects. Actually, everything we want to do is we have never done before. So we will try our best.

Week 4

The problems we faced may be harder than our imagination. First, the search space is too large to traverse. The number of potential topologies is , while K is the number of interaction forms and n is the number of species. It is an astronomy number if n >=3. Second, the non-linear system identification is still an unsolved problem in the present. The parameters can not be recognized exactly even the topologies are meeting the requirements. Some approximation should be done mathematically. First, the networks we designed are restricted to the simple network, which means the reactors are at most two kinds in one reaction. The searching space can be shrunk in this way, while all the dynamic behaviors can be simulated with this simple model. Second, we may need some restrictions to the interactions forms and the parameters. The more information is, the better system we can get by the design process.

Week 6

The frame of the software has been constructed as indicated in right figure.

week 9

From the analysis, we should get the information of a reaction from three aspects: topology, formula and the parameters. The optimized algorithms are required with the three different layers problems. For the bottom layer, the parameters are optimized with the known forms of formula. For the upper layer, we need genetic algorithm to look for the formula forms that meet the requirement. And for the top layer, genetic algorithm is also employed to optimize the topology. And I will work with the selection of the formula and the topology.

week 16

The genetic algorithm has been implemented with c++ and tested with some benchmarks. Next, we discussed about the data structures for the program, including the input form and input data, the output form and output data, the description of a reaction, and the items considering in formula. According to the discussion results, I began to implement the whole frame for the program. There are three classes in the program.

week 20

The frame of the program has been implemented. The oscillator system is tested.

week 22

information is integrated. And the adaptation system is tested

week 23

PSO Algorithm is proposed and implemented.

week 25

The multiple initial values system identification is integrated. And the toggle switch system is tested.

week 27

Wikiwikiwiki~ Fast learning needed to digest so much info...

G/LSA and PSO (Yuwei Cui)

Feb. 12, 2009

GUI (Xiaomo Yao)

Jun. 15, 2009

Start GUI design

Jul. 26, 2009 Start GUI project Vedio Commu.

Jul. 27, 2009 Form created, with MDI & dock windows GUI version 090727 released

Aug. 4, 2009 XML editor finished

Aug. 13, 2009 Start developing layout drawing canvas

Aug. 22, 2009 Layout drawing canvas developed; Reactions and modifiers enabled

Aug. 26, 2009 Layout drawing canvas completed

Aug. 27, 2009 GUI version 090827 released

Aug. 31, 2009 Start programming interface with identification/analysis algorithms

Oct. 5, 2009 Interface developed

Oct. 16, 2009 Interface completed GUI version 091017 released

I/O SBML (Jiahao Li)

Feb. 12, 2009

No matter what kind of software we will make, it is sure that one of the most important work is to find out the structure and functions of existing software. Preparing start from Athena and Tinkercell, two multi-function synthetic biological software. OMG!! these software are great! They can do not only some simulating works, but also other create job! Their GUI also performs well. It seems hard to beyond or even achieve them in simulating functions. We have to find other highlights.

Mar. 14, 2009

Be aware of SBML, a standard structured file contains a certain chemical/biological system. The support of SBML related functions in Athena; Tinkercell and other software convince me the necessity of providing these functions in our software. It means our output can supply a standard interface to others.

May. 2, 2009

Start learning SBML and LibSBML. Well~~ download, unzip, install, OMG!! Quite a lot of codes, my teammates advise me to read the tutorial first. Yes, after reading, I realize the structure of SBML file and working process of LibSBML approximately. In terms of our software, some basic information of the SBML file is needed, such as amount of spceies and parameters, their names, information of each reaction, including the mathematical expressions of each kineticlaw. All these data need to extract from given SBML file. So it looks learning to use LibSBML is inevitable.

Jun. 6, 2009

Go on working with LibSBML, Successfully output the species names, parameters and other informations.

Need to face some problems when deal with the mathematical expressions, SBML uses ASTNode to store the formula, although it is not difficult to manage the "tree" data structure, I still need to make clear the functions in "class ASTNode". My teammates need me to supply the information orderly, which means :

1. The information of one reaction should be stored by categories.

2. The terms(the formula related to one certain species) should be put together.

3. The whole information need to be packaged for next step.

Aug. 15, 2009

The readSBMLFile part finished, in this part, I pick each species list in the file orderly, and search all reactions, if the species appeared in "ListOfReactants"; "ListOfProducts" or "ListOfModifiers", i'll push the informations of the reaction into the vertor<term>, the data structure of vector<term> has been decided as "class val_func". After the traversing, i'll package these vectors<term>, and sent to the next step.

I/O Database and Miscellaneous(Yu He)

Apr. 3, 2009

Release of bmptodat plug-in, in order to realize 'desire in, solution out' Source Code File:FuncExtractor.cpp Function: Extract function info from a .bmp file Input: .bmp file path and name (less than 200 bytes) | .dat file path and name (less than 200 bytes) Output: .dat file with range x[0,1], y[0,1], and step length of 0.01 (alterable). Demo:

1.Draw the picture D:\Activities\iGEM\Codes\weigh.bmp like this:

2.Input like this: 3.Get a .dat file at under folder D:\Activities\iGEM\Codes, named weigh.dat like this:

4.Testify the output weigh.dat file by Origin8.0:

To be improved: 1.Rescale-able range of x-axis and y-axis 2.User designated step length 3.Function of ‘smoothing’

Apr. 15, 2009

Release of matrixtoode plug-in
Source Code File:Matrix_ODEArray_conv_1.0.cpp
Function:Convert a interaction-style matrix into text-ified ODE Array
Input: a coefficient matrix through console window
Output: .txt file describing the ODE Array in matlab readable text.
Demo:

1.Input the coefficient matrix like this:

2.Output completed like this:

Note: Depreciated later on.

May 15, 2009

Today the core algorithm of reverse engineering is proposed. We discussed possible approaches to realize these system identification requirements.

The very first thing emerge in my mind is Metropolis method and Monte Carol sampling. If we define an error function which serves as 'judging criteria', and a random number generator used for 'correction' or 'exploration', then we can basically get it running from pure theoretical view.

May 20, 2009

Fourth order multi-dimensional Runge Kutta Method is realized. Tested also on Mathematica7.0 to precision of 6th. digits.

Jul. 15, 2009

Hooray!!! Raw program of system identification on an colony-evolution system, whose behavior is damped oscillation realized.(3 functions, 6 parameters)

Aug. 5, 2009

Raw simulation of an oscillation example on substential level. (6 functions, 5 parameters, nonlinear model).

Sept. 1, 2009

Realizes multiple (13 functions, 11 parameters) targets identification with basic metropilis method and gradient decrease method. Yet, it is still a linear func array.

Sept. 3, 2009

Database in development. It's really a challenge to adopt SQLite now. I don't have that much confidence, since there's little time left.

Sept. 6, 2009

Today I finally make up my mind to realize the database from 'ab-initio' on. Since the schedule is so tight, I still feel safer to work with familiar C language.

Sept. 20, 2009

Okay, the very basic, essential, elementary database in .dat form is finished...

Sept. 28, 2009

Inter-database search and sort function completed.

Oct. 8, 2009

Thanks to Grapeot's timely suggestion and help. The new version of bmptodat plug-in is successfully upgraded by taking color-depth into consideration. This eliminates the x-axis 'triple folding' bug due to differnt color-depth of 8 bits and 24 bits.

Automatic Biological Circuits Design

Team Logo: wanna know more about the hinding metaphors and inspirations in this little red square? Click to check out how much fun this year's iGEM has brought us!