This repository has been archived by the owner on Nov 21, 2024. It is now read-only.
-
Notifications
You must be signed in to change notification settings - Fork 11
/
README.txt
136 lines (91 loc) · 6.82 KB
/
README.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
# moFF : A modest Feature Finder (but still robust) to extract apex MS1 itensity directly from Thermo raw file
================================
moFF is written in python and it is based on a Go library that is able to read raw file from Thermo machine
Required library :
Python 2.7
pandas > 0.17.
numpy > 1.9.0
argparse > 1.2.1
scikit-learn > 0.17
moFF is composed by two stand alone modules :
moff_mbr.py : matching between run
moff.py : apex intensity
To run the entire workflow (mbr and apex ) you should use moff_all.py.
moFF uses txic to extract the XiC data from the raw files, so the execute txic must be located in the same folder where you have all moFF scripts.
The txic program is compatibale with the raw file of all the Orbitrap and triple quadrupole Thermo machines.
For the moment it does not work with the Thermo Fusion machine.
The input files that contain the list of the MS2 identified peptides (you can use any search engines) must contains the information showed in moFF_setting..property for each peptide.
moFF_setting.property : it specifies the minimun specificic requirements of the input files tha are :
-- tab delimited file
-- the header of the infput file should contains the following the fields and columnns names :
'peptide' : sequence of the peptide
'prot': protein ID
'rt': retention time of peptide
'mz' : mass over charge
'mass' : mass of the peptide
'charge' : charge of the ionized peptide
see the sample input files in the folder f1_folder.
The retention time must be specified in second
In the folder f1_folder you have three input files, that contain the MS2 identified peptides (sing MASCOT) of three runs (three tecnical replicates ) from the CPTAC study 6.
you can download the relative raw files from https://goo.gl/ukbpCI, in order to run the next examples.
Matching between run module:
use : python moff_mbr.py -h
--inputF LOC_IN specify the folder of the input MS2 peptide files [REQUIRED]
--sample SAMPLE specify which replicate files are used fot mbr [regular expr. are valid]
--ext EXT specify the exstension of the input file (txt as default value)
--log_file_name LOG_LABEL a label name for the log file (moFF_mbr.log as default log file name)
--filt_width W_FILT iwidth value of the filter (k * mean(Dist_Malahobis) , k = 2 as default)
--out_filt OUT_FLAG filter outlier in each rt time allignment (active as default)
--weight_comb W_COMB weights for model combination combination : 0 for no weight (default) 1 weighted devised by model errors.
python moff_mbr.py --inputF f1_folder/
It runs the mbr modules and save the output files in a subfolder called 'mbr_output' inside the folder given in input.
The mbr module will take all the .txt files in your input folder as replicates. (to select specific files or different extension see below))
In the f1_folder/mbr_output you will find the same number of the input files, but they will have a new field called 'matched' that specifies which peptides are matched (1) or the not (0)
The rt field of the matched peptide contains the predicted rt retentioins time.
if your input files inside your working fodler have another exstension like (.list, etc) you can use :
use : python --inputF f1_folder/ --ext list ( Do not specify '.list' but only 'list')
if you need to select specific input files from your working folder ( choose ) , you can use an regular expression as:
use : python --inputF f1_folder/ --sample *_6A (you can also use --ext option if you need)
the mbr will output a log file (moFF_mbr.log as default log file name) with all the details and it is saved inside the --inputF given in inout
Apex module:
use python moff.py -h
--input NAME specify the input file with the of MS2 peptides
--tol TOLL specify the tollerance parameter in ppm
--rt_w RT_WINDOW specify rt window for xic (minute). Default value is 3 min
--rt_p RT_P_WINDOW specify the time windows for the peak ( minute). Default value is 0.1
--rt_p_match RT_P_WINDOW_MATCH specify the time windows for the matched peptide peak ( minute). Default value is 0.4
--raw_repo RAW specify the raw file repository
--output_folder LOC_OUT specify the folder output
python moff.mbr --input f1_folder/20080311_CPTAC6_07_6A005.txt --raw_rep f1_folder/ --tol 1O
it run the apex module on the input file , extraxing the apex intesity from the respective raw file in folder --raw_repo.
In the output files, moFF just add the following fields to your origin input file:
"intensity" intensity, taking the highest peak in the XIC
"rt_peak" rt of the highest peak
"lwhm" left width half maximun of the signal in seconds
"rwhm" right width half maximun of the signal in seconds
"SNR" signal-to-noise
"log_L_R" log ratio of lwhm over rwhm (peak shape )
"log_int" log 2 of the intesity
It generates a .log file (with same name of input file ) that contains detailesd information for each peak retrieved.
This module determines automaticaly if the input file contains matched peptides or not.
REMARK : the raw file names MUST be the same of the input file otherwise the script give you an error !
python moff.mbr --input f1_folder/20080311_CPTAC6_07_6A005.txt --raw_rep f1_folder/ --tol 1O --output_folder output_moff
It will put the results in the folder output_moff
Run the entire workflow (Mbr + Apex ) :
use python moff_all.py
--inputF LOC_IN specify the folder of the input MS2 peptide list files
--sample SAMPLE specify witch replicated use for mbr reg_exp are valid
--ext EXT specify the file extentention of the input like
--log_file_name LOG_LABEL a label name to use for the log file
--filt_width W_FILT width value of the filter k * mean(Dist_Malahobis)
--out_filt OUT_FLAG filter outlier in each rt time allignment
--weight_comb W_COMB weights for model combination combination : 0 for no weight 1 weighted devised by trein err of the model.
--tol TOLL specify the tollerance parameter in ppm
--rt_w RT_WINDOW specify rt window for xic (minute). Default value is 3 min
--rt_p RT_P_WINDOW specify the time windows for the peak ( minute). Default value is 0.1
--rt_p_match RT_P_WINDOW_MATCH specify the time windows for the matched peptide peak ( minute). Default value is 0.4
--raw_repo RAW specify the raw file repository
--output_folder LOC_OUT specify the folder output
python moff_all.py --inputF f1_folder/ --raw_repo f1_folder/ --output_folder output_moff
The option are the same of the two modules, the the output mbr files are stores in the folder f1_folder/mbr_output and the result of the apex module are stored in output_moff
Also the log files are stored in the respective folders