Insights to the Binding of a Selective Adenosine A 3 Receptor Antagonist Using Molecular Dynamic Simulations, MM-PBSA and MM-GBSA Free Energy Calculations, and Mutagenesis

• Published in: Journal of chemical information and modeling 2019-12, Vol.59 (12), p.5183-5197
• Main Authors: Lagarias, Panagiotis, Barkan, Kerry, Tzortzini, Eva, Stampelou, Margarita, Vrontaki, Eleni, Ladds, Graham, Kolocouris, Antonios
• Format: Article
• Language: English
• ISSN: 1549-9596; 1549-960X
• DOI: 10.1021/acs.jcim.9b00751

Adenosine A receptor (A R) is a promising drug target cancer and for a number of other conditions like inflammatory diseases, including asthma and rheumatoid arthritis, glaucoma, chronic obstructive pulmonary disease, and ischemic injury. Currently, there is no experimentally determined structure of A R. We explored the binding profile of O4-{[3-(2,6-dichlorophenyl)-5-methylisoxazol-4-yl]carbonyl}-2-methyl-1,3-thiazole-4-carbohydroximamide (K18), which is a new specific and competitive antagonist at the orthosteric binding site of A R. MD simulations and MM-GBSA calculations of the WT A R in complex with K18 combined with in vitro mutagenic studies show that the most plausible binding conformation for the dichlorophenyl group of K18 is oriented toward trans-membrane helices (TM) 5, 6 and reveal important residues for binding. Further, MM-GBSA calculations distinguish mutations that reduce or maintain or increase antagonistic activity. Our studies show that selectivity of K18 toward A R is defined not only by direct interactions with residues within the orthosteric binding area but also by remote residues playing a significant role. Although V169 is considered to be a selectivity filter for A R binders, when it was mutated to glutamic acid, K18 maintained antagonistic potency, in agreement with our previous results obtained for agonists binding profile investigation. Mutation of the direct interacting residue L90 in the low region and the remote L264 in the middle/upper region to alanine increases antagonistic potency, suggesting an empty space in the orthosteric area available for increasing antagonist potency. These results approve the computational model for the description of K18 binding at A R, which we previously performed for agonists binding to A R, and the design of more effective antagonists based on K18.