Boron containing compounds (BCCs) aroused increasing interest in the scientific community due to their wide application as drugs in various fields. In order to design new compounds hopefully endowed with pharmacological activity and also investigate their conformational behavior, the support of computational studies is crucial. Nevertheless, the suitable molecular mechanics parameterization and the force fields needed to perform these simulations are not completely available for this class of molecules. In this paper, Amber force field parameters for phenyl-, benzyl-, benzylamino-, and methylamino-boronates, a group of boron-containing compounds involved in different branches of the medicinal chemistry, were created. The robustness of the obtained data was confirmed through molecular dynamics simulations on ligand/β-lactamases covalent complexes. The ligand torsional angles, populated over the trajectory frames, were confirmed by values found in the ligand geometries, located through optimizations at the DFT/B3LYP/6-31g(d) level, using water as a solvent. In summary, this study successfully provided a library of parameters, opening the possibility to perform molecular dynamics simulations of this class of boron-containing compounds.
Keywords: MD simulations; amber; boron; covalent ligand; force field; paramfit; β-lactamases.