Computational studies on horseshoe shape pocket of human orexin receptor type 2 and boat conformation of suvorexant by molecular dynamics simulations

Interaction between the Human OX2 Orexin Receptor and Suvorexant and Some of Its Analogues: SAPT (DFT) Interaction Energy Decomposition Analysis

In this study, the interaction energy (E_int) of suvorexant (as an orexin receptor antagonist) and some of its analogues with the important residues of the human OX2 orexin receptor, determined by molecular docking, is calculated using the symmetry-adapted perturbation theory-density functional theory (SAPT (DFT)) method. Also, the important residues with the dominant interaction with each ligand are determined based on the obtained SAPT (DFT) interaction energies. To analyze the interaction of the receptor with each ligand, the decomposition of E_int to its constituent components including electrostatic (E_ele), exchange (E_ex), induction (E_ind), and exchange-induction (E_ex-ind), dispersion (E_disp), and exchange-dispersion (E_ex-disp) is performed. The change of interaction energy components with the replacement of the benzoxazole part of suvorexant by pyrimidine containing different functional groups, thieno pyrimidine, and furo pyrimidine is also investigated, separately. It is found that the change in E_int, due to these replacements, is controlled more by the variation of the electrostatic interaction energy component of E_int than by the other interaction energy components. A linear correlation (R² = 0.91) is found for the variation of E_int versus experimental ligand-binding affinities. Also, the existence of the linear correlation for the variation of the interaction energy components with experimental ligand-binding affinities is investigated. The variation of the electrostatic component versus experimental ligand-binding affinities shows a more linear correlation compared to the other interaction energy components.

Insights From Molecular Dynamics Simulations of a Number of G-Protein Coupled Receptor Targets for the Treatment of Pain and Opioid Use Disorders

Effective treatments for pain management remain elusive due to the dangerous side-effects of current gold-standard opioid analgesics, including the respiratory depression that has led to skyrocketing death rates from opioid overdoses over the past decade. In an attempt to address the horrific opioid crisis worldwide, the National Institute on Drug Abuse has recently proposed boosting research on specific pharmacological mechanisms mediated by a number of G protein-coupled receptors (GPCRs). This research is expected to expedite the discovery of medications for opioid overdose and opioid use disorders, leading toward a safer and more effective treatment of pain. Here, we review mechanistic insights from recent all-atom molecular dynamics simulations of a specific subset of GPCRs for which high-resolution experimental structures are available, including opioid, cannabinoid, orexin, metabotropic glutamate, and dopamine receptor subtypes.

Molecular Dynamics Exploration of Selectivity of Dual Inhibitors 5M7, 65X, and 65Z toward Fatty Acid Binding Proteins 4 and 5

Designing highly selective inhibitors of fatty acid binding proteins 4 and 5 (FABP4 and FABP5) is of importance for treatment of some diseases related with inflammation, metabolism, and tumor growth. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were performed to probe binding selectivity of three inhibitors (5M7, 65X, and 65Z) to FABP4/FABP5 with K_i values of 0.022/0.50 μM, 0.011/0.086 μM, and 0.016/0.12 μM, respectively. The results not only suggest that all inhibitors associate more tightly with FABP4 than FABP5, but also prove that the main forces driving the selective bindings of inhibitors to FABP4 and FABP5 stem from the difference in the van der Waals interactions and polar interactions of inhibitors with two proteins. Meanwhile, a residue-based free energy decomposition method was applied to reveal molecular basis that inhibitors selectively interact with individual residues of two different proteins. The calculated results show that the binding difference of inhibitors to the residues (Phe16, Phe19), (Ala33, Gly36), (Phe57, Leu60), (Ala75, Ala78), (Arg126, Arg129), and (Tyr128, Tyr131) in (FABP4, FABP5) drive the selectivity of inhibitors toward FABP4 and FABP5. This study will provide great help for further design of effective drugs to protect against a series of metabolic diseases, arteriosclerosis, and inflammation.