# HG changeset patch # User galaxyp # Date 1588786551 14400 # Node ID 2133b0be850afa377d728af99c08ea7b67b08422 # Parent 7f432d87c82cb2d8a8c5ec0a0075bdcfdfb0640f "planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/maxquant commit 4e77eeda8a112fb50af00325a5164b986c16fc5c" diff -r 7f432d87c82c -r 2133b0be850a maxquant.xml --- a/maxquant.xml Wed Apr 15 11:17:42 2020 -0400 +++ b/maxquant.xml Wed May 06 13:35:51 2020 -0400 @@ -1,4 +1,4 @@ - + @@ -81,15 +81,26 @@ calcPeakProperties: ${search_opts.calc_peak_properties} #set writeMzTab = "mzTab" in $output_opts.output writeMzTab: $writeMzTab - #if $lfq_opts.do_ibaq.ibaq == 'True': - ibaq: ${lfq_opts.do_ibaq.ibaq} - ibaqLogFit: ${lfq_opts.do_ibaq.ibaqLogFit} + #if $protein_quant.lfq_opts.do_ibaq.ibaq == 'True': + ibaq: ${protein_quant.lfq_opts.do_ibaq.ibaq} + ibaqLogFit: ${protein_quant.lfq_opts.do_ibaq.ibaqLogFit} #end if - separateLfq: ${lfq_opts.separateLfq} - lfqStabilizeLargeRatios: ${lfq_opts.lfqStabilizeLargeRatios} - lfqRequireMsms: ${lfq_opts.lfqRequireMsms} - advancedSiteIntensities: ${lfq_opts.advancedSiteIntensities} + separateLfq: ${protein_quant.lfq_opts.separateLfq} + lfqStabilizeLargeRatios: ${protein_quant.lfq_opts.lfqStabilizeLargeRatios} + lfqRequireMsms: ${protein_quant.lfq_opts.lfqRequireMsms} + advancedSiteIntensities: ${protein_quant.lfq_opts.advancedSiteIntensities} matchBetweenRuns: ${search_opts.match_between_runs} + includeContaminants: ${search_opts.incl_contaminants} + quantMode: ${protein_quant.peptides_for_quantification} + restrictProteinQuantification: ${protein_quant.only_unmod_prot.unmod_prot} + #if $protein_quant.only_unmod_prot.unmod_prot == 'True': + #if $protein_quant.only_unmod_prot.mods_used_prot_quant: + restrictMods: [${protein_quant.only_unmod_prot.mods_used_prot_quant}] + #else: + restrictMods: [] + #end if + useCounterparts: $protein_quant.only_unmod_prot.discard_unmod_cpart_peptides + #end if paramGroups: #for $pg in $paramGroups: #set names = [re.sub('@SUBSTITUTION_RX@', '_', str($n.element_identifier)) for $n in $pg.files] @@ -111,6 +122,7 @@ #else: enzymes: [] #end if + enzymeMode: ${pg.digestion_mode} #if $pg.quant_method.select_quant_method == 'silac': labelMods: #if $pg.quant_method.light_labels: @@ -139,9 +151,10 @@ #if $pg.quant_method.select_quant_method == 'reporter_ion_ms2': lcmsRunType: 'Reporter ion MS2' reporterMassTolerance: 0.003 - reporterPif: 0 + reporterPif: ${pg.quant_method.reporter_pif} reporterFraction: 0 reporterBasePeakRatio: 0 + filterPif: ${pg.quant_method.filter_pif} isobaricLabels: #if $pg.quant_method.iso_labels.labeling == 'custom': #for $l in $pg.quant_method.iso_labels.iso_label: @@ -216,7 +229,7 @@ + help="Specify one or more FASTA databases."/> @@ -239,49 +252,133 @@
+ instructions see the help text."/> + label="maximum peptide mass [Da]" value="4600" + help="Peptides that are heavier than this mass will be discarded in the Andromeda search."/> + label="minimum unique peptides" value="0" + help="The minimum number of unique peptides a protein group should have to be considered as identified and reported in the final table."/> + label="Calculate peak properties" + truevalue="True" falsevalue="False" + help="If checked, several quantities characterizing peaks and isotopes patterns are calculated. This may lead to a substantial increase in computation time."/> + label="Match between runs" + truevalue="True" falsevalue="False" + help="Identifications are transferred to non-sequenced or non-identified MS features in other LC-MS runs."/> + +
- + +
+ + + + + + + + + + + + + + + + + + + + + + + +
+ + + + + + + + + + + + + + + +
+ +
+ + + label="missed cleavages" value="2" + help="The number of missed cleavages that are maximally tolerated in the in-silico digestion of the protien sequences."/> - + multiple="true" + help="Select zero or more fixed modifications. They will always be attached to any occurence of the respective amino acid."> + + multiple="true" + help="Select zero or more variable modifications. Do not specify label modifications here, neither ms1 level labels, like SILAC, nor isobaric labels."> - + + multiple="true" + help="Select zero or more enzymes. The enzymes used for generating the in silico peptides for the Andromeda search."> - + + + + + + + + + + @@ -292,33 +389,34 @@ - - - - - - - - - + + + + + + + + + - - - - + + + + @@ -334,6 +432,7 @@ + @@ -345,62 +444,37 @@ - - + + - - + + - + + + -
- - - - - - - - - - - - - - - -
+ -
+ label="Only create config files? (i.e. for testing)" + truevalue="True" falsevalue=""/> @@ -430,6 +504,7 @@ + @@ -464,6 +539,14 @@ + +
+ + + + + +
@@ -485,6 +568,8 @@ + + @@ -497,6 +582,8 @@ + + diff -r 7f432d87c82c -r 2133b0be850a mqparam.py --- a/mqparam.py Wed Apr 15 11:17:42 2020 -0400 +++ b/mqparam.py Wed May 06 13:35:51 2020 -0400 @@ -308,6 +308,23 @@ .format(key)) node.text = str(value) + def set_list_param(self, key, values): + """Set a list parameter. + Args: + key: (string) XML tag of the parameter + values: the lit of values of the parameter XML node + Returns: + None + """ + node = self._root.find(key) + if node is None: + raise ValueError('Element {} not found in parameter file' + .format(key)) + node.clear() + node.tag = key + for e in values: + et_add_child(node, name='string', text=e) + def _from_yaml(self, conf): """Read a yaml config file. Args: @@ -317,6 +334,7 @@ """ with open(conf) as f: conf_dict = yaml.safe_load(f.read()) + paramGroups = conf_dict.pop('paramGroups') self.add_infiles([pg.pop('files') for pg in paramGroups]) for i, pg in enumerate(paramGroups): @@ -339,7 +357,10 @@ else: raise Exception('No fasta files provided.') for key in conf_dict: - self.set_simple_param(key, conf_dict[key]) + if key in ['restrictMods']: + self.set_list_param(key, conf_dict[key]) + else: + self.set_simple_param(key, conf_dict[key]) def write(self, mqpar_out): """Write pretty formatted xml parameter file. diff -r 7f432d87c82c -r 2133b0be850a test-data/01/combined/txt/parameters.txt --- a/test-data/01/combined/txt/parameters.txt Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/01/combined/txt/parameters.txt Wed May 06 13:35:51 2020 -0400 @@ -18,9 +18,9 @@ Min. razor peptides 1 Min. peptides 1 Use only unmodified peptides and True -Modifications included in protein quantification Oxidation (M);Acetyl (Protein N-term) +Modifications included in protein quantification Peptides used for protein quantification Razor -Discard unmodified counterpart peptides True +Discard unmodified counterpart peptides False Label min. ratio count 2 Use delta score False iBAQ False diff -r 7f432d87c82c -r 2133b0be850a test-data/01/config.yml --- a/test-data/01/config.yml Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/01/config.yml Wed May 06 13:35:51 2020 -0400 @@ -13,10 +13,16 @@ lfqRequireMsms: True advancedSiteIntensities: True matchBetweenRuns: False + includeContaminants: True + quantMode: 1 + restrictProteinQuantification: True + restrictMods: [] + useCounterparts: True paramGroups: - files: ['BSA_min_23.mzXML'] maxMissedCleavages: 1 fixedModifications: [Carbamidomethyl (C)] variableModifications: [Oxidation (M)] enzymes: [Trypsin/P] + enzymeMode: 0 \ No newline at end of file diff -r 7f432d87c82c -r 2133b0be850a test-data/01/mqpar.xml --- a/test-data/01/mqpar.xml Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/01/mqpar.xml Wed May 06 13:35:51 2020 -0400 @@ -77,16 +77,13 @@ 1 1 0 - False + True True False 2 True - - Oxidation (M) - Acetyl (Protein N-term) - + 0 0 0 diff -r 7f432d87c82c -r 2133b0be850a test-data/02/config.yml --- a/test-data/02/config.yml Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/02/config.yml Wed May 06 13:35:51 2020 -0400 @@ -13,17 +13,24 @@ lfqRequireMsms: True advancedSiteIntensities: True matchBetweenRuns: False + includeContaminants: False + quantMode: 0 + restrictProteinQuantification: True + restrictMods: [Pro5,Pro6] + useCounterparts: True paramGroups: - files: ['BSA_min_23.mzXML'] maxMissedCleavages: 2 fixedModifications: [] variableModifications: [Oxidation (M)] enzymes: [Trypsin/P] + enzymeMode: 0 lcmsRunType: 'Reporter ion MS2' reporterMassTolerance: 0.003 - reporterPif: 0 + reporterPif: 0.6 reporterFraction: 0 reporterBasePeakRatio: 0 + filterPif: False isobaricLabels: - ['TMT2plex-Lys126','TMT2plex-Nter126',0.0,0.0,6.7,3.0,True] - ['TMT2plex-Lys127','TMT2plex-Nter127',0.0,0.0,0.0,0.0,True] @@ -32,11 +39,13 @@ fixedModifications: [] variableModifications: [] enzymes: [] + enzymeMode: 0 lcmsRunType: 'Reporter ion MS2' reporterMassTolerance: 0.003 - reporterPif: 0 + reporterPif: 0.8 reporterFraction: 0 reporterBasePeakRatio: 0 + filterPif: True isobaricLabels: - [iTRAQ4plex-Lys114,iTRAQ4plex-Nter114,0,0,0,0,False] - [iTRAQ4plex-Lys115,iTRAQ4plex-Nter115,0,0,0,0,False] diff -r 7f432d87c82c -r 2133b0be850a test-data/02/mqpar.xml --- a/test-data/02/mqpar.xml Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/02/mqpar.xml Wed May 06 13:35:51 2020 -0400 @@ -27,7 +27,7 @@ True revert all - True + False 4600 True True @@ -77,15 +77,15 @@ 1 1 0 - False + True True False 2 True - Oxidation (M) - Acetyl (Protein N-term) + Pro5 + Pro6 0 0 @@ -93,7 +93,7 @@ 0 15 0 - 1 + 0 200 @@ -267,7 +267,7 @@ 0 0 0.003 - 0 + 0.6 False 0 0 @@ -442,8 +442,8 @@ 0 0 0.003 - 0 - False + 0.8 + True 0 0 0 diff -r 7f432d87c82c -r 2133b0be850a test-data/03/config.yml --- a/test-data/03/config.yml Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/03/config.yml Wed May 06 13:35:51 2020 -0400 @@ -15,12 +15,18 @@ lfqRequireMsms: True advancedSiteIntensities: True matchBetweenRuns: False + includeContaminants: True + quantMode: 1 + restrictProteinQuantification: True + restrictMods: [] + useCounterparts: True paramGroups: - files: ['BSA_min_22.mzxml'] maxMissedCleavages: 1 fixedModifications: [Carbamidomethyl (C)] variableModifications: [Oxidation (M)] enzymes: [Trypsin/P] + enzymeMode: 0 lfqMode: 1 lfqSkipNorm: True lfqMinEdgesPerNode: 3 diff -r 7f432d87c82c -r 2133b0be850a test-data/03/mqpar.xml --- a/test-data/03/mqpar.xml Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/03/mqpar.xml Wed May 06 13:35:51 2020 -0400 @@ -77,16 +77,13 @@ 1 1 1 - False + True True False 2 True - - Oxidation (M) - Acetyl (Protein N-term) - + 0 0 0 diff -r 7f432d87c82c -r 2133b0be850a test-data/04/config.yml --- a/test-data/04/config.yml Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/04/config.yml Wed May 06 13:35:51 2020 -0400 @@ -13,12 +13,18 @@ lfqRequireMsms: True advancedSiteIntensities: True matchBetweenRuns: False + includeContaminants: True + quantMode: 1 + restrictProteinQuantification: True + restrictMods: [] + useCounterparts: True paramGroups: - files: ['BSA_min_23.mzXML', 'BSA_min_22.mzxml'] maxMissedCleavages: 2 fixedModifications: [Carbamidomethyl (C)] variableModifications: [Oxidation (M)] enzymes: [Trypsin/P] + enzymeMode: 0 labelMods: - [Arg6,Lys4] - [] diff -r 7f432d87c82c -r 2133b0be850a test-data/04/mqpar.xml --- a/test-data/04/mqpar.xml Wed Apr 15 11:17:42 2020 -0400 +++ b/test-data/04/mqpar.xml Wed May 06 13:35:51 2020 -0400 @@ -77,16 +77,13 @@ 1 1 0 - False + True True False 2 True - - Oxidation (M) - Acetyl (Protein N-term) - + 0 0 0