Network simulations and inferring network connectivity, Simulation in MATLAB

The objective of this problem is to simulate reaction networks for different connectivity. Another objective is to guess the network connectivity from given experimental data. Matlab code has been provided to you (attached) for a large number of reaction networks. Your first task is to just run the matlab code for a variety of networks and examine the trends in the results. We have split the matlab code provided to you into two parts (i) function files that contain the differential equations describing a network e.g. "reactionsABCDE.m" and (ii) corresponding script files that run the function for a given set of parameters e.g. run_reactionsABCDE.m

Tasks to perform

(a) Simulation of Linear Reaction Networks:

Run the matlab code provided to you for the following linear reaction networks (A → B, A → B → C, A→ B → C → D, and A→ B → C → D → E). Comment on the trends for the time dependent profiles for the starting reactant in the network, the intermediates and the end-points in the network. Develop an intuition for the time dependent behavior of various components?

Can you guess the sequence of a linear reaction network if you are given the data for individual species just by looking at the trends? You will use the heuristic information on such trends to reduce the search space in the network connectivity inference problem (part c) on more complicated branched networks with inhibition and activation.

(b) Simulation of Branched Reaction Networks

Run the matalb code provided to you for the following branched reaction networks (i) ABC_AD which is linear ABC with a branching from A to D (ii) ABC_BD which is linear ABC with a branching from B to D. Examine the branched reaction network shown below for the species A through G. Write down the corresponding differential equations and create a function file.

Create script files with different parameter values and simulate the branched network below.

Provide the simulated results for the case where all ki have the value of 1.0

(c) Inferring Network Connectivity from Experimental Data. Our objective is to guess the network connectivity from available experimental data.

The following experimental data are provided to you (as plots and as tables). Based on experience from forward simulation of known networks, can you guess the connectivity of the network? There is a species in one network that is inhibiting another species in the second network. Assume all k values to be 1. Guess network connectivity, simulate the network and try to figure out which metabolite is inhibiting which reaction. Assume that the coefficient of inhibition follows the form: 1/(1+[C]^n) where [C] is the concentration of the inhibiting metabolite and 'n' is the Hill coefficient. The function xlsread can be used to input the data from the EXCEL file into matlab if desired.

Posted Date: 2/15/2013 4:40:13 AM | Location : United States







Related Discussions:- Network simulations and inferring network connectivity, Assignment Help, Ask Question on Network simulations and inferring network connectivity, Get Answer, Expert's Help, Network simulations and inferring network connectivity Discussions

Write discussion on Network simulations and inferring network connectivity
Your posts are moderated
Related Questions
Anybody who can help me with a project using MATLAB on the financial data ( pe ratios etc of stocks? basically do simulation and MATLAB? or neone of u know any person who can do th

Q. Using MATLAB, plot the gain and phase response of the transfer function

implemnt a sliding window model in matlab for mining data stream(e.g transactional database) to find out closed frequent patterns. Mining is performed only over latest W transactio

The advertising director for Diversity Paint and Supply, a chain of four retail stores on Chicago’s North Side, is considering two media possibilities. One plan is for a series of


This question examines the addition of a (synthetic) echo to the voice recording, such as would be found in a reverberant room. To do this, use a difference equation of the form

Create a Matlab function called question3 which finds the root of an equation numerically by using the following recursive formula: xn+1 = x n – f(x n)/g(x n) , for n =1,2,3,……

Q. A periodic sequence of exponential wave forms forms a pulse train whose first cycle is represented by v(t) = [u(t) - u(t - T 0 )]vAe -t/TC Use MATLAB to find Vrms of the

The objective of this problem is to simulate reaction networks for different connectivity. Another objective is to guess the network connectivity from given experimental data. Matl

Using the audiorecorder() function in MATLAB, make a recording of your voice. Use a sample rate of 8kHz, 16 bits per sample, and aim to record around 2 seconds worth. Quantize t