
****************************************************************************************************************************************

Comments are included in all scripts as a guidance through the code. Equations 1-12 are implemented in the scripts in
different order comparing to the paper, namely, first active component of GTPases (Ca, Ra, Rhos), then inactive component of 
GTPases (Ci, Ri, Rhoi), followed by Par (Pc, Pun) and Scribble (Sc, Sun) complexes and active and inactive integrin (Ia, Ii) at the end.


Open source code include several Matlab scripts:


1) Polarization.m - Calculate response of GTPases and Par and Scribble complexes for input values defined in Table 1 and 2. 
It allows activation by integrin at left and right side of the cell close to the membrane (equation 13) by defining l=1 or r=1.
Effect of the amount of activated integrin on polarization can be analysed with this code.
The code produces figures 4 and 7 from the paper.


2) EMT.m - Takes into account degradation or upregulation (options: downregulation=1 or upregulation=1, respectively; equation 16) 
of each active component of GTPases (option 1 for Rac1, option 2 for Rho and option 3 for Cdc42). Responses of GTPases and
Par and Scribble complexes for each case are plotted with specific color maps.
The code produces figures 3, 8 and 9 from the paper.


3) MPSA_PDE.m, RunSim.m and Plot_figures.m are used for sensitivity analysis.
MPSA_PDE.m is used for generation of different sets of selected 21 input parameters (range of variation -
five times smaller to five times higher values in comparison to the reference set defined in Table 1 and 2).
The script creates sets of dummy parameters as well and defines the criteria of acceptable and unacceptable 
output values.
To perform Monte Carlo simulations, the code calls RunSim.m script.
For plotting figures of the obtained results Plot_figures.m script is called.
Set of results obtained with these scripts for 3000 simulations is included and can be used to plot results
form Plot_figures.m without need to run new simulations.
The code produces figures 5 and 6 from the paper.

*****************************************************************************************************************************************


