Calculate protein interface ΔΔG values for a collection of proteins complexes given a list of mutations of interest for each complex. FoldX's AnalyseComplex command is run on each WT complex, then the required mutated models are predicted with BuildModel and analysed again with AnalyseComplex to determine the difference in binding strength.
- FoldX 5
- Python 3
- Snakemake
- Numpy
- Pandas
- Clone this repo
- Place your complex model PDB files as
data/complex/{complex id}/model.pdb
- Add an
individual_list
file to each complex directory, as specified by FoldX's BuildModel - Run Snakemake
A number of simple modifications can be easily made:
- Generate all possible interface mutants based on the output of the WT interface analysis (see my Mutfunc: SARS-CoV-2 pipeline)
- Increase RAM requirements for individual complexes (see
get_mut_complex_ram
inSnakefile
) - Analyse interfaces other than between chain A and B (see
get_complex_tsv_files
inSnakefile
)