Package: pssm-promoter-tool
PSSM Promoter ToolThe tool applies CORPSE (Codon Restrained Promoter Silencing) method and inverted CORPSE (iCORPSE) to the provided gene sequence.-35 and -10 promoters, along with the additional non-canonical sequence motifs, are predicted based on the Salis Lab Promoter Calculator (https://github.com/hsalis/SalisLabCode/tree/master/Promoter_Calculator).
Position-specific scoring matrix (PSSM) is applied to all the synonymous codon variants of the promoters associated with the lowest and highest transcription rates in order to maximally decrease (CORPSE) or increase the transcription rate (iCORPSE).
The output CSV file/files contain synonymous codon promoters and sequence motifs for the minimal and maximal transcriptional rates, along with the non-canonical sequence motifs for forward and reverse strands.For a web (Google Colab) version, please navigate tohttps://colab.research.google.com/drive/171iBNCrA1hMS-LpX_qaFani34HiueTZO?usp=sharing.INSTALLATION:Install required libraries using pip:pip install pssm-promoter-toolDownload and unpack the archive with files using "Code"->"Download ZIP" buttons in the right corner athttps://github.com/ellinium/pssm_promoter_tool.
Or usegit clone https://github.com/ellinium/pssm-promoter-tool[! it's necessary to have a Git account and be authorised in the system to run the command !]USAGE:The tool requires a text or fasta file with a nucleotide sequence of a gene to process.
From the folder with the downloaded files, run:python pssm_promoter_calculator.py <file_name>where 'file_name' is a path to the file with a gene sequence (TXT or FASTA format).By default, the command processes 40 PSSM promoter combinations with the highest PSSM scores for transcription rate maximisation or 40 combinations with the lowest PSSM scores for transcription rate minimisation.
To run all the promoter combinations, use an additional argument 'all':python pssm_promoter_calculator.py <file_name> allDepending on the result, up to four output CSV files can be generated:PSSMPromoterCalculator_MIN_FWD_results.csv - contains promoters to minimise transcription rate (forward strand)PSSMPromoterCalculator_MIN_REV_results.csv - contains promoters to minimise transcription rate (reverse strand)PSSMPromoterCalculator_MAX_FWD_results.csv - contains promoters to maximise transcription rate (forward strand)PSSMPromoterCalculator_MAX_REV_results.csv - contains promoters to maximise transcription rate (reverse strand).The output file fields in the CSV files contain data from Salis' Promoter calculator and additional fields:Type - 'Original Promoter' - for original sequence promoters; 'Modified Promoter' - for synonymous promoters.TSS -  transcriptional start Site position (nt)Tx_rate - transcription unitiation rate (au)Tx_rate_FoldChange - the fold change between the original transcription rate and re-calculated with PSSM primersframe - a reading frame for the promotershex35 - -35 hexamer sequence (an upstream 6-nucleotide site called the −35 motif)PSSM_hex35 - position-specific scoring matrix value for the -35 motifhex35_sequence - contains 9nt sequence that includes -hex35 promoter if it's in the 2nd or 3rd framehex35_aa - an amino acid sequence for the -35 motifhex10 - -10 hexamer sequence (a downstream 6-nucleotide site called the −10 motif)PSSM_hex10 - position-specific scoring matrix value for the -10 motifhex10_sequence - contains 9nt sequence that includes -hex35 promoter if it's in the 2nd or 3rd framehex10_aa - an amino acid sequence for the -10 motifUP - an upstream sequence (UP). A 20-nucleotide region that appears upstream of the −35 motif.spacer - a spacer sequence that separates the −10 and −35 motifsdisc - discriminator sequence. A typically 6-nucleotide region in between the −10 motif and TSS.ITR - the first 20 transcribed nucleotides called the initial transcribed region (ITR)new_gene_sequence - contains a gene sequence (nt) with substituted (Type = 'Modified Promoter') or  original promoters (Type = 'Original Promoter')promoter_sequence - contains -35 motif, spacer and - 10 motifdG_total - total Gibbs free energy change (kcal/mol)dG_10 - Gibbs binding free Energy for -10 hexamer (kcal/mol)dG_35 - Gibbs binding free Energy for -35 hexamer (kcal/mol)dG_disc - Gibbs free energy penalty for non-optimal discriminator element (kcal/mol)dG_ITR - Gibbs free energy change for R-loop formation at the initial transcribed region (kcal/mol)dG_ext10 - Gibbs binding free energy for the extended -10 hexamer (kcal/mol)dG_spacer - Gibbs free energy penalty for non-optimal spacing (kcal/mol)dG_UP - Gibbs binding free energy for the upstream element (kcal/mol)dG_bind - binding Gibbs free energyUP_position - a position of the UP elementhex35_position - a position of the -35 motifspacer_position - a position of the spacerhex10_position - a position of the -10 motifdisc_position - a position of the discriminatorReferences:Logel DY, Trofimova E, Jaschke PR. Codon-Restrained Method for Both Eliminating and Creating Intragenic Bacterial Promoters. ACS Synth Biol. 2022 Jan 19;acssynbio.1c00359. Available fromhttps://pubs.acs.org/doi/10.1021/acssynbio.1c00359. doi: 10.1021/acssynbio.1c00359LaFleur TL, Hossain A, Salis HM. Automated model-predictive design of synthetic promoters to control transcriptional profiles in bacteria. Nat Commun. 2022 Sep 2;13(1):5159. Available fromhttps://www.nature.com/articles/s41467-022-32829-5. doi: 10.1038/s41467-022-32829-5