Research Article

AUTO-MUTE 2.0: A Portable Framework with Enhanced Capabilities for Predicting Protein Functional Consequences upon Mutation

Table 1

Prediction tables obtained for a sample text file of mutants (mutants.txt) using stability_changes_ddG.pl. For the initial run, the flag -t 0 enables the automatic retrieval of PDB files, which are saved in the main AUTO-MUTE 2.0 folder. The flag -m 0 requests that predictions (sample_output_ddG_RF.txt) be made using our previously trained random forest (RF) supervised classification model. For each requested mutant, a predicted stability category and a confidence measure (%) for the prediction are provided in the table under column headings labeled Stability and Confid., respectively. The flag -t 1 is used with the second run since the PDB files have already been downloaded, and the flag -m 2 requests that predictions (sample_output_ddG_REPTree.txt) be made using our previously trained tree regression (REPTree) model. In this case, a predicted stability change value for each requested mutant is provided in the table under the column heading labeled ddG. A glossary of terms (Glossary.txt) defining all table column headings to the right of the predictions is provided in the main folder. Briefly, parameters Vol., sT, Loc., and Num. characterize the mutated position C-alpha vertex with respect to the tessellation of the protein structure, while SS refers to the secondary structure at the mutated position, and these parameters are also discussed in Section 2.2 of the Methods. Not shown here are similar predictions based on two additional models that are available with this stability_changes_ddG.pl program: support vector machine classification (-m 1) and support vector regression (-m 3). Details for all the programs are available from the README file in the main folder, as well as by typing “perl <program_name>” without any flags at the command-line and hitting enter.
(a)

mutants.txt

3phvAD25E257
3phvAT31A257
4lyzAA31V257
4lyzAD101R257
1benDL15Q257
1rn1CQ25K257
1g35BH69N257
1g35BL90M257
1pohAK49E257

(b)

C:∖AutoMute2>perl stability_changes_ddG.pl -t 0 -f mutants.txt -m 0 -o sample_output_ddG_RF.txt

sample_output_ddG_RF.txt

PDB_IDChainMutationStabilityConfid.Temp.pHVol.sTLoc.Num.SS

3phvAD25EDecreased0.6225714.60.16B0C
3phvAT31ADecreased0.6425714.10.14B0S
4lyzAA31VIncreased0.8625719.50.10B0H
4lyzAD101RIncreased0.8525717.70.18S6C
1benDL15QDecreased0.9025710.90.13U1H
1rn1CQ25KIncreased0.962578.70.14S6H
1g35BH69NDecreased0.6725715.20.18S4S
1g35BL90MDecreased0.8525718.90.13B0H
1pohAK49EIncreased1.0025710.10.12S4H

(c)

C:∖AutoMute2>perl stability_changes_ddG.pl -t 1 -f mutants.txt -m 2 -o sample_output_ddG_REPTree.txt

sample_output_ddG_REPTree.txt

PDB_IDChainMutationddGTemp.pHVol.sTLoc.Num.SS

3phvAD25E−1.9525714.60.16B0C
3phvAT31A−1.3825714.10.14B0S
4lyzAA31V1.2325719.50.10B0H
4lyzAD101R0.2425717.70.18S6C
1benDL15Q−1.6425710.90.13U1H
1rn1CQ25K0.612578.70.14S6H
1g35BH69N−0.6025715.20.18S4S
1g35BL90M−0.9025718.90.13B0H
1pohAK49E1.3725710.10.12S4H