Molecular dynamics (MD) simulation of our homology model of CBTS with GGPP were embedded in a TIP3P water box (water density 0.997 g/ml, water probe radius 1.4 Å) and neutralized with a 0.9 mM NaCl concentration, and simulated using periodic boundary conditions at
T = 298 K using the YASARA2-force field, an 1 fs integration time, a
NPT ensemble and a cutoff of 8 Å for non-bonded interactions. Long-range coulomb electrostatics were treated using the Particle-Mesh-Ewald approach (PME) with a grid spacing of < 1 Å. The size of the simulation box was 125.00 × 156.83 × 123.79 Å
3. Charges and force field parameters for GGPP and cations A-C were obtained from the AutoSmiles force field parameter assignment is
YASARA Structure [23] (
link). QM-optimized structure of cation A from Hong
et al. [24] (
link) was docked in CBTS using the AutoDockVina program implemented in
YASARA Structure using a rectangular 14 × 14 × 14 Å
3 simulation cell around residue R316, V340, N347, S451, A487, D499, A764 and Y836. The docking simulations were performed treating all atoms as rigid, and binding modes were characterized based on a cluster analysis. Cations B and C were relaxed within the active site by energy minimization and a 1 ns MD-simulation using the YASARA2 force field [23] (
link). MD simulations were performed using
YASARA Structure (version 15.8.31, YASARA Biosciences, [23] (
link)a).
Schrepfer P., Ugur I., Klumpe S., Loll B., Kaila V.R, & Brück T. (2020). Exploring the catalytic cascade of cembranoid biosynthesis by combination of genetic engineering and molecular simulations. Computational and Structural Biotechnology Journal, 18, 1819-1829.