Molecular dynamics simulation of mammalian 15S-lipoxygenase with AMBER force field

A molecular dynamics simulation study of mononuclear iron 15S-lipoxygenase (15S-LOX) from rabbit reticulocytes was performed to investigate its structure and dynamics; newly developed AMBER force field parameters were employed for the first coordination sphere of the catalytic iron (II). The results obtained from this study demonstrate that the structural features of the catalytic iron coordination site are in good agreement with available data obtained from experiments. The motional flexibility of the N-terminal β-barrel domain is greater than the C-terminal catalytic domain; flexibility was assessed in terms of B-factors and secondary structure calculations. The significant features obtained for the relative motional flexibility of these two domains of 15S-LOX in solution as well as the isolated C-terminal domain were analyzed in terms of radius of gyration and maximum diameter, which correlated well with the structural flexibility of 15-lipoxygenase-1 in solution as probed by small-angle X-ray scattering. The motional flexibility indicates interdomain motion between the N-terminal β-barrel and the C-terminal catalytic domain; this was further verified by the evaluation of central bending in the solvated LOX molecule, which identified an unstructured stretch of amino acids as the interdomain linker. The average bending angle confirmed significant central bending between these two domains, which was linked to the high degree of motional freedom of the N-terminal β-barrel domain in aqueous solutions. This can be considered to have biological relevance for membrane binding as well as for regulating the catalytic domain.