Thermotoga neapolitana β-glucosidase Download PDB file

Protein: 1,4-beta-D-glucan glucohydrolase
Organism: Thermotoga neapolitana
Lenght: 425
UniprotID: P0C946 (Download sequence)
Template: 1OD0:A (90.59%) | [φ,ψ] plot



Acid/base catalytic
Conserved (50-79%)
Conserved (> 80%)
Extrapoled mutation
Catalytic pocket (0-6.5Å from ligand)
α-helix | 3-helix | π-helix
β-bridge | β-ladder | β-strand
Loop | bend | hydrogen bonded turn
Community 1
Community 2

Extrapoled mutations

| Cellobiose docking

Poses (avg) Affinity (wild) Affinity (mutant) ΔAffinity ΔAffinity expected Status
1 -6.07 -6.30 -0.23 ΔAffinity < 0 ✓
1-3 -5.80 -6.10 -0.3 ΔAffinity < 0 ✓
1-10 -5.01 -5.09 -0.08 ΔAffinity < 0 ✓

| Glucose docking

Poses (avg) Affinity (wild) Affinity (mutant) ΔAffinity ΔAffinity expected Status
1 -6.07 -5.70 0.37 ΔAffinity > 0 ✓
1-3 -5.80 -5.67 0.13 ΔAffinity > 0 ✓
1-10 -5.17 -5.39 -0.22 ΔAffinity > 0 x

| Cellobiose docking

Poses (avg) Affinity (wild) Affinity (mutant) ΔAffinity ΔAffinity expected Status
1 -6.07 -6.03 0.04 ΔAffinity < 0 x
1-3 -5.80 -5.82 -0.02 ΔAffinity < 0 ✓
1-10 -5.01 -4.97 0.04 ΔAffinity < 0 x

| Glucose docking

Poses (avg) Affinity (wild) Affinity (mutant) ΔAffinity ΔAffinity expected Status
1 -6.07 -5.50 0.57 ΔAffinity > 0 ✓
1-3 -5.80 -5.37 0.43 ΔAffinity > 0 ✓
1-10 -5.17 -5.07 0.1 ΔAffinity > 0 ✓



Methodology

Comparative modeling

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.

Minimization

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).

Multiple alignment

We performed multiple alignment using Clustal Omega (default parameters). We used Python scripts to detect conservation.

Docking

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.