I am working with a small molecule ligand that contains selenium. I am having trouble parameterizing it using CGenFF/ParamChem. I can't get it to successfully run through any of the parameterization options in CHARMM-GUI either. Does anyone have *any* experience with/advice regarding parameterizing and running MD of selenium containing small molecules? Any hints/suggestions would be very much appreciated!
With a quick web search I found
Dataset of AMBER force field parameters of drugs, natural products and steroids for simulations using GROMACS.
Jennifer Loschwitz, Anna JÃ¤ckering, Monika Keutmann, Maryam Olagunju, Olujide O.Olubiyi, Birgit Strodel
Data in Brief,Volume 35, April 2021, 106948
They used GAFF to derive models using S instead of Se, defined some modified Se parameters, and changed S to Se after all the model building, which may be an approach to use with CGenFF and CHARMM-GUI. The Se parameters in that paper may or may not be transferrable, and should be validated before use with CHARMM force fields.
Thank you very much for sharing this nice idea from the paper. This is perhaps what I will do as well. Seems like Se parameters are not available in many of the popular forcefields.