DFT calculations were used to obtain parameters compatible with the CHARMM36 force field for iron-sulfur clusters (Fe-S) of the type [Fe4S4]2+ that are coordinated to dissimilatory adenosine-5'-phosphosulfate reductase (APSrAB). Classical molecular dynamics (MD) simulations were performed on two APSrAB systems to validate the parameters and verify the stability of the studied systems. The time analysis of the parameters inserted into the force field was in reasonable agreement with the experimental X-ray diffraction data. The analysis of the time evolution of the studied systems indicated that these systems and, in particular, the clusters in their respective cavities had a good stability and were in agreement with what was observed in previous works. The parameters obtained provide the basis for the study of APSrAB as well as other systems that contain [Fe4S4]2+ through the CHARMM36 force field.
Keywords: DFT; Dissimilatory adenosine-5’-phosphosulfate reductase; force field parametrization; iron-sulfur clusters; molecular dynamics simulation.