CHARMM Force-Field Parameters for Morphine, Heroin, and Oliceridine, and Conformational Dynamics of Opioid Drugs

Opioid drug binding to specialized G protein-coupled receptors (GPCRs) can lead to analgesia upon activation via downstream G(1) protein signaling and to severe side effects via activation of the beta-arrestin signaling pathway. Knowledge of how different opioid drugs interact with receptors is essential, as it can inform and guide the design of safer therapeutics. We performed quantum and classical mechanical computations to explore the potential energy landscape of four opioid drugs: morphine and its derivatives heroin and fentanyl and for the unrelated oliceridine. From potential energy profiles for bond twists and from interactions between opioids and water, we derived a set of force-field parameters that allow a good description of structural properties and intermolecular interactions of the opioids. Potential of mean force profiles computed from molecular dynamics simulations indicate that fentanyl and oliceridine have complex energy landscapes with relatively small energy penalties, suggesting that interactions with the receptor could select different binding poses of the drugs.