Squalene–hopene cyclase catalyzes the cyclization of squalene
to hopanoids. A previous study has identified a network of tunnels
in the protein, where water molecules have been indicated to move.
Blocking these tunnels by site-directed mutagenesis was found to change
the activation entropy of the catalytic reaction from positive to
negative with a concomitant lowering of the activation enthalpy. As
a consequence, some variants are faster and others are slower than
the wild type (wt) in vitro under optimal reaction conditions for
the wt. In this study, molecular dynamics (MD) simulations have been
performed for the wt and the variants to investigate how the mutations
affect the protein structure and the water flow in the enzyme, hypothetically
influencing the activation parameters. Interestingly, the tunnel-obstructing
variants are associated with an increased flow of water in the active
site, particularly close to the catalytic residue Asp376. MD simulations
with the substrate present in the active site indicate that the distance
for the rate-determining proton transfer between Asp376 and the substrate
is longer in the tunnel-obstructing protein variants than in the wt.
On the basis of the previous experimental results and the current
MD results, we propose that the tunnel-obstructing variants, at least
partly, could operate by a different catalytic mechanism, where the
proton transfer may have contributions from a Grotthuss-like mechanism.