Optimization of CHARMM force field parameters for ryanodine receptor inhibitory drug dantrolene using FFTK and FFParam

Ryanodine receptors (RyRs) are large intracellular ligand-gated calcium release ion channels. Mutations in human RyR1 in combination with a volatile anesthetic or muscle relaxant are known to cause leaky RyRs resulting in malignant hyperthermia (MH). This has long been primarily treated with the RyR inhibitory drug dantrolene. Alternatives to dantrolene as a RyR inhibitor may be found through computer-aided drug design. Additionally, molecular dynamics (MD) studies of dantrolene interacting with RyRs may reveal its full mechanism of action. The availability of accurate force field parameters is important for the success of both.

In this study, force field parameters for dantrolene were obtained from the CHARMM General Force Field (CGenFF) program and optimized using the force field toolkit (FFTK) and FFParam programs. The obtained parameters were then validated by a comparison between calculated and experimental IR spectra and normal mode analysis, among other techniques.