Poses (avg) | Affinity (wild) | Affinity (mutant) | ΔAffinity | ΔAffinity expected | Status |
---|---|---|---|---|---|
1 | -6.60 | -6.67 | -0.07 | ΔAffinity < 0 | ✓ |
1-3 | -6.56 | -6.61 | -0.050000000000001 | ΔAffinity < 0 | ✓ |
1-10 | -6.15 | -6.16 | -0.0099999999999998 | ΔAffinity < 0 | ✓ |
Poses (avg) | Affinity (wild) | Affinity (mutant) | ΔAffinity | ΔAffinity expected | Status |
---|---|---|---|---|---|
1 | -6.60 | -5.80 | 0.8 | ΔAffinity > 0 | ✓ |
1-3 | -6.56 | -5.60 | 0.96 | ΔAffinity > 0 | ✓ |
1-10 | -5.39 | -5.43 | -0.04 | ΔAffinity > 0 | x |
Three-dimensional structures were modeled by homology using MODELLER. Sequences were collected in UniProt. We searched for the E.C. number of beta-glucosidases (3.2.1.21) and collected sequences of family GH1 (~3000). For each sequence, we selected the structure of highest identity using BLAST against the PDB database to be used as a template. Sequences with identity lower than 25% were discarded. 100 models were constructed for each sequence. The best model was selected using the lowest value of the DOPE score. Mutants were modeled using MODELLER's mutation script. The mutations were based on beneficial mutations described in the literature.
PDB models were minimized using AMBER tools. We performed 1000 steps of minimization in the vacuum (parameters: imin=1, maxcyc=1000, ncyc=750, cut=10.0, ntb=0, igb=0
).
We performed multiple alignment using Clustal Omega (default parameters). We used Python scripts to detect conservation.
We performed docking of glucose (product) and cellobiose (substrate) for wild and mutant beta-glucosidases using AutoDock Vina (parameters: exhaustiveness = 20; num_modes = 10
). The box center was defined using the median coordinate of the last atom of the two catalytic glutamates. The experiment was performed in triplicate. We collected ten poses for each docking and analyzed the average of affinity energies predicted by the software for (i) pose 1, (ii) poses 1-3, and (iii) poses 1-10.