The Performance of Different Water Models on the Structure and Function of Cytochrome P450 Enzymes

Abstract

Modelling approaches and modern simulations to investigate the biomolecular structure and function rely on various methods, one among which is the choice of water model. Water molecules play a crucial role in all sorts of chemistry. Cytochrome P450 (CYP450), in particular, water molecules are crucial for the formation of active oxidants which perform the oxidation and metabolism of several substrates. Computational chemistry tools such as MD simulations and QM/MM calculations, nowadays, have become complementary tools to study the structure and functions of CYP450 enzymes, and therefore, accurate modeling of water molecules is crucial. In the present study, we have highlighted the behavior of the three most widely used water models-TIP3P, SPC/E, and OPC for three different CYP450 enzymes-CYP450BM3, CYP450OleT, and CYP450BSβ during MD simulations and QM/MM calculations. We studied the various properties such as RMSD, RMSF, H-bond, water occupancy in the first solvation shell, and Hydrogen Atom Transfer (HAT) using QM/MM calculations and compared them for all the three water models.Our study shows that, the stability of the enzymes structure is well maintained in all the three water models. However, OPC water model performs well for the polar active sites, i.e., in CYP450OleT, CYP450BSβ while the TIP3P water model is superior for the hydrophobic site such as CYP450BM3 .

Molecular Dynamics Simulations Reveal the Conformational Transition of GH33 Sialidases

“…A hydrogen bond (gmx_hbond) was defined if the acceptor–donor distance was less than 0.35 nm and the acceptor–hydrogen–donor angle was less than 30° [ 46 ]. The RMSD (gmx_rms) and RMSF (gmx_rmsf) were determined to examine the structural stability and flexibility after eliminating the overall translational and rotational movements by superimposing the C α atoms of each snapshot structure onto the initial structure using least-squares fitting [ 47 , 48 ]. The interaction energy (gmx_energy) between amino acid residues and substrate was calculated based on the CHARMM36 all-atom force field [ 49 ].…”

Multiscale QM/MM modelling of catalytic systems with ChemShell

“…2) are among the most studied systems using QM/MM simulations with ChemShell, and continue to be a very active topic of research. [36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51] In recent work Kalita et al used a combination of molecular dynamics simulations and QM/MM to explore oxygen binding and the protonation steps leading to formation of compound I in two archetypal P450 systems. 52 The team used these techniques also to understand the behaviour of a bioengineered CYP450 BME variant that catalyses a reaction (C-H amination) that natural P450s do not.…”

Theoretical study on hydrogen transfer in the dissociation of dimethyl disulfide radical cations

“…Because water molecules do not directly participate in these three HTs and they influence free energy barriers through electrostatic interaction, and even H-bonding interaction. 33,34 These may also influence the detailed process in the reaction. For example, H7 transfers from the C4 to the S2 atom and the S2-C4 bond forms after TS4, while H7 still transfers from the C4 to the S2 atom but the C1-S3 bond forms after TS4-w, as the H-bond formed with a water molecule is conducive to stabilize charge and thus stabilize a molecule.…”