z-logo
open-access-imgOpen Access
Combined docking and machine learning identify key molecular determinants of ligand pharmacological activity on β2 adrenoceptor
Author(s) -
JiménezRosés Mireia,
Morgan Bradley Angus,
Jimenez Sigstad Maria,
Tran Thuy Duong Zoe,
Srivastava Rohini,
Bunsuz Asuman,
BorregaRomán Leire,
Hompluem Pattarin,
Cullum Sean A.,
Harwood Clare R.,
Koers Eline J.,
Sykes David A.,
Styles Iain B.,
Veprintsev Dmitry B.
Publication year - 2022
Publication title -
pharmacology research and perspectives
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.975
H-Index - 27
ISSN - 2052-1707
DOI - 10.1002/prp2.994
Subject(s) - g protein coupled receptor , docking (animal) , drug discovery , agonism , computational biology , agonist , antagonism , chemistry , ligand (biochemistry) , antagonist , receptor , functional selectivity , pharmacology , stereochemistry , biology , medicine , biochemistry , nursing , politics , political science , law
G protein‐coupled receptors (GPCRs) are valuable therapeutic targets for many diseases. A central question of GPCR drug discovery is to understand what determines the agonism or antagonism of ligands that bind them. Ligands exert their action via the interactions in the ligand binding pocket. We hypothesized that there is a common set of receptor interactions made by ligands of diverse structures that mediate their action and that among a large dataset of different ligands, the functionally important interactions will be over‐represented. We computationally docked ~2700 known β2AR ligands to multiple β2AR structures, generating ca 75 000 docking poses and predicted all atomic interactions between the receptor and the ligand. We used machine learning (ML) techniques to identify specific interactions that correlate with the agonist or antagonist activity of these ligands. We demonstrate with the application of ML methods that it is possible to identify the key interactions associated with agonism or antagonism of ligands. The most representative interactions for agonist ligands involve K97 2.68×67 , F194 ECL2 , S203 5.42×43 , S204 5.43×44 , S207 5.46×641 , H296 6.58×58 , and K305 7.32×31 . Meanwhile, the antagonist ligands made interactions with W286 6.48×48 and Y316 7.43×42 , both residues considered to be important in GPCR activation. The interpretation of ML analysis in human understandable form allowed us to construct an exquisitely detailed structure‐activity relationship that identifies small changes to the ligands that invert their pharmacological activity and thus helps to guide the drug discovery process. This approach can be readily applied to any drug target.

The content you want is available to Zendy users.

Already have an account? Click here to sign in.
Having issues? You can contact us here