Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Messages - azadeh.s

Pages: [1]
1
General Discussion / Re: "Affinity Profile"
« on: February 22, 2019, 10:15:54 am »
Hello

I'am also novice for using your great package REDUCE_Suit.
I'm interested to study the affinity profile for a vary famous gene COL1A1( I downloded COL1A1 from Ensembel database) I tried to understand the dynamic of affinity profile for a transcription factor Sp1 through its binding site for both wild and mutant COL1A1 genes .

To start my study I ran the commande :
AffinityProfile \
   -prefix=aff_ \
   -threshold=0 \
   -psam=../../Desktop/source/REDUCE-Suite-v2.2/examples/AffinityProfile/psam_001.xml \
   -sequence=sequence_COL1A1.fasta

As a Psam matrix, I used the matrix given in your example.
But then I don't understand the result.

cat seq_psam.dat
   psam_001.xml
NM_000088.3   0.0301094
####################
cat AffinityProfile.log
sequence file: sequence_COL1A1.fasta   cut-off: 0
reading in sequences from file [sequence_COL1A1.fasta]
reading in PSAM ../../Desktop/source/REDUCE-Suite-v2.2/examples/AffinityProfile/psam_001.xml [1/1]
[  1]: ../../Desktop/source/REDUCE-Suite-v2.2/examples/AffinityProfile/psam_001.xml (forward); pmax=0.025626

Time used: 00:00:00:00

#################
cat aff_psam_001.xml

NM_000088.3   0.0301094   1:2.78962e-104   2:2.8034e-101   3:2.508e-07   4:2.97891e-94   5:3.73255e-103   6:4.75891e-116   7:6.73541e-70   8:6.52039e-71   9:5.81238e-65   10:4.04049e-53   11:9.0064e-52   12:2.29734e-51   13:1.30952e-52   14:1.14509e-48   15:2.33791e-51   16:6.36074e-41   17:4.78319e-46   18:9.88362e-103   19:1.18759e-104   20:7.56195e-118   21:3.05295e-99   22:1.19914e-103   23:2.68201e-21   24:6.17375e-91   25:6.30612e-151   26:6.59257e-50   27:1.50964e-89   28:1.19137e-86   29:9.42612e-47   30:5.0609e-66   31:3.89768e-149   32:1.41179e-49   33:1.62746e-79   34:2.80803e-108   35:5.26342e-77   36:1.28883e-55   37:4.3514e-24   38:5.40903e-86   39:5.26603e-61   40:2.5533e-83   41:2.91799e-121   42:4.92685e-68   43:1.42043e-72   44:3.66844e-54   45:1.07013e-92   46:2.2805e-56   47:4.31802e-112   48:2.12206e-67   49:4.32273e-65   50:5.36831e-108   51:3.10334e-59   52:5.76347e-71   53:3.189e-74   54:1.16991e-15   55:2.31613e-120   56:4.17858e-69   57:3.11132e-70   58:6.68288e-27   59:1.72033e-105   60:7.81592e-63   61:7.91623e-64   62:9.43495e-63   63:1.05312e-74   64:9.53001e-123   65:2.04295e-66   66:5.71577e-30   67:5.90979e-108   68:8.76473e-54   69:2.51005e-73   70:2.08671e-81   71:9.14769e-38   72:5.51473e-31   73:2.33375e-44   74:1.11588e-46   75:1.80495e-39   76:3.66284e-49   77:1.52272e-44   78:4.28484e-71   79:8.19242e-23   80:7.09162e-57   81:1.31163e-58   82:7.40765e-44   83:1.63155e-65   84:6.85533e-61   85:1.67216e-60   86:3.10819e-77   87:7.99735e-26   88:2.20357e-81   89:2.67345e-53   90:7.94344e-34   91:2.97164e-83   92:5.35325e-70   93:4.00621e-67   94:9.56203e-61   95:8.90908e-50   96:2.06665e-60   97:4.78235e-53   98:2.39696e-50   99:1.73158e-52   100:4.63302e-32   101:1.96568e-84   102:9.10009e-69   103:1.6031e-73   104:8.72453e-44   105:3.37299e-88   106:5.02218e-44   107:1.70837e-107   108:6.30181e-110   109:6.10561e-89   110:6.96046e-35   111:2.86401e-101   112:2.9393e-31   113:4.69209e-129   114:2.34756e-56   115:9.01049e-68   116:3.00809e-54   ...
...
..
...
5906:9.52601e-72   5907:9.93085e-28   5908:1.0446e-76   5909:1.53848e-65   5910:6.02333e-52   5911:2.12289e-43   5912:7.00049e-71   5913:7.89907e-37   5914:8.86224e-60   5915:1.47167e-70


####################
cat   seq_psam.dat
   psam_001.xml
NM_000088.3   0.0301094


#########################

Could you please help me to understand these results ?
What is this matrix  "aff_psam_001.xml" ?
How could I obtain more details on the affinity change and dynamics of binding site ?

Thank you in advance
Best Regards

Azadeh

Pages: [1]
Created and maintained by Dr. Xiang-Jun Lu [律祥俊]. See also http://forum.x3dna.org and http://x3dna.org