The development of an Amber-compatible organosilane force field for drug-like small molecules

As members of the group IVA elements, silicon and carbon have long been thought of as isosteres of each other in drug design. However, the lack of silicon parameters in current main stream force fields hinders the computational study of this important element in drug discovery. Thus, in this study, we attempted to supplement the parameters of organosilanes in the General Amber Force Field (GAFF2). The parameters have been designed following the principles of GAFF2 to make it compatible with the Amber force field family. The accuracy of the parameters was discussed by comparing the pair interaction energy, the liquid properties, and the structures and alchemical binding free energy differences for a set of protein–ligand complexes.