1
|
Pawnikar S, Magenheimer BS, Joshi K, Nevarez-Munoz E, Haldane A, Maser RL, Miao Y. Activation of polycystin-1 signaling by binding of stalk-derived peptide agonists. eLife 2024; 13:RP95992. [PMID: 39373641 PMCID: PMC11458180 DOI: 10.7554/elife.95992] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/08/2024] Open
Abstract
Polycystin-1 (PC1) is the protein product of the PKD1 gene whose mutation causes autosomal dominant Polycystic Kidney Disease (ADPKD). PC1 is an atypical G protein-coupled receptor (GPCR) with an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal and membrane-embedded C-terminal (CTF) fragments. Recently, activation of PC1 CTF signaling was shown to be regulated by a stalk tethered agonist (TA), resembling the mechanism observed for adhesion GPCRs. Here, synthetic peptides of the first 9- (p9), 17- (p17), and 21-residues (p21) of the PC1 stalk TA were shown to re-activate signaling by a stalkless CTF mutant in human cell culture assays. Novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) simulations elucidated binding conformations of p9, p17, and p21 and revealed multiple specific binding regions to the stalkless CTF. Peptide agonists binding to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of stalk TA-mediated PC1 CTF activation. Additional sequence coevolution analyses showed the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. These insights into the structural dynamic mechanism of PC1 activation by TA peptide agonists provide an in-depth understanding that will facilitate the development of therapeutics targeting PC1 for ADPKD treatment.
Collapse
Affiliation(s)
- Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of KansasLawrenceUnited States
| | - Brenda S Magenheimer
- Clinical Laboratory Sciences, University of Kansas Medical CenterKansas CityUnited States
- The Jared Grantham Kidney Institute, University of Kansas Medical CenterKansas CityUnited States
| | - Keya Joshi
- Department of Pharmacology and Computational Medicine Program, University of North CarolinaChapel HillUnited States
| | - Ericka Nevarez-Munoz
- Clinical Laboratory Sciences, University of Kansas Medical CenterKansas CityUnited States
| | - Allan Haldane
- Department of Physics, and Center for Biophysics and Computational Biology, Temple UniversityPhiladelphiaUnited States
| | - Robin L Maser
- Clinical Laboratory Sciences, University of Kansas Medical CenterKansas CityUnited States
- The Jared Grantham Kidney Institute, University of Kansas Medical CenterKansas CityUnited States
- Department of Biochemistry and Molecular Biology, University of Kansas Medical CenterKansas CityUnited States
| | - Yinglong Miao
- Department of Pharmacology and Computational Medicine Program, University of North CarolinaChapel HillUnited States
| |
Collapse
|
2
|
Liu Y, Wang K, Cao F, Gao N, Li W. Interactions between Inhibitors and 5-Lipoxygenase: Insights from Gaussian Accelerated Molecular Dynamics and Markov State Models. Int J Mol Sci 2024; 25:8295. [PMID: 39125865 PMCID: PMC11311652 DOI: 10.3390/ijms25158295] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2024] [Revised: 07/27/2024] [Accepted: 07/28/2024] [Indexed: 08/12/2024] Open
Abstract
Inflammation is a protective stress response triggered by external stimuli, with 5-lipoxygenase (5LOX) playing a pivotal role as a potent mediator of the leukotriene (Lts) inflammatory pathway. Nordihydroguaiaretic acid (NDGA) functions as a natural orthosteric inhibitor of 5LOX, while 3-acetyl-11-keto-β-boswellic acid (AKBA) acts as a natural allosteric inhibitor targeting 5LOX. However, the precise mechanisms of inhibition have remained unclear. In this study, Gaussian accelerated molecular dynamics (GaMD) simulation was employed to elucidate the inhibitory mechanisms of NDGA and AKBA on 5LOX. It was found that the orthosteric inhibitor NDGA was tightly bound in the protein's active pocket, occupying the active site and inhibiting the catalytic activity of the 5LOX enzyme through competitive inhibition. The binding of the allosteric inhibitor AKBA induced significant changes at the distal active site, leading to a conformational shift of residues 168-173 from a loop to an α-helix and significant negative correlated motions between residues 285-290 and 375-400, reducing the distance between these segments. In the simulation, the volume of the active cavity in the stable conformation of the protein was reduced, hindering the substrate's entry into the active cavity and, thereby, inhibiting protein activity through allosteric effects. Ultimately, Markov state models (MSM) were used to identify and classify the metastable states of proteins, revealing the transition times between different conformational states. In summary, this study provides theoretical insights into the inhibition mechanisms of 5LOX by AKBA and NDGA, offering new perspectives for the development of novel inhibitors specifically targeting 5LOX, with potential implications for anti-inflammatory drug development.
Collapse
Affiliation(s)
- Yuyang Liu
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (Y.L.); (K.W.)
| | - Kaiyu Wang
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (Y.L.); (K.W.)
| | - Fuyan Cao
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China;
| | - Nan Gao
- Key Laboratory of Polyoxometalate and Reticular Material Chemistry of Ministry of Education, Faculty of Chemistry, Northeast Normal University, Changchun 130024, China
| | - Wannan Li
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (Y.L.); (K.W.)
| |
Collapse
|
3
|
Pradhan B, Pavan M, Fisher CL, Salmaso V, Wan TC, Keyes RF, Rollison N, Suresh RR, Kumar TS, Gao ZG, Smith BC, Auchampach JA, Jacobson KA. Lipid Trolling to Optimize A 3 Adenosine Receptor-Positive Allosteric Modulators (PAMs). J Med Chem 2024; 67:12221-12247. [PMID: 38959401 DOI: 10.1021/acs.jmedchem.4c00944] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2024]
Abstract
A3 adenosine receptor (A3AR) positive allosteric modulators (PAMs) (2,4-disubstituted-1H-imidazo[4,5-c]quinolin-4-amines) allosterically increase the Emax of A3AR agonists, but not potency, due to concurrent orthosteric antagonism. Following mutagenesis/homology modeling of the proposed lipid-exposed allosteric binding site on the cytosolic side, we functionalized the scaffold, including heteroatom substitutions and exocyclic phenylamine extensions, to increase allosteric binding. Strategically appended linear alkyl-alkynyl chains with terminal amino/guanidino groups improved allosteric effects at both human and mouse A3ARs. The chain length, functionality, and attachment position were varied to modulate A3AR PAM activity. For example, 26 (MRS8247, p-alkyne-linked 8 methylenes) and homologues increased agonist Cl-IB-MECA's Emax and potency ([35S]GTPγS binding). The putative mechanism involves a flexible, terminally cationic chain penetrating the lipid environment for stable electrostatic anchoring to cytosolic phospholipid head groups, suggesting "lipid trolling", supported by molecular dynamic simulation of the active-state model. Thus, we have improved A3AR PAM activity through rational design based on an extrahelical, lipidic binding site.
Collapse
Affiliation(s)
- Balaram Pradhan
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, United States
| | - Matteo Pavan
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, United States
| | - Courtney L Fisher
- Department of Pharmacology & Toxicology and the Cardiovascular Center, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, Wisconsin 53226, United States
| | - Veronica Salmaso
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, United States
- Department of Pharmaceutical and Pharmacological Sciences, University of Padua, 35131 Padua, Italy
| | - Tina C Wan
- Department of Pharmacology & Toxicology and the Cardiovascular Center, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, Wisconsin 53226, United States
| | - Robert F Keyes
- Department of Biochemistry and the Program in Chemical Biology, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, Wisconsin 53226, United States
| | - Noah Rollison
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, United States
| | - R Rama Suresh
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, United States
| | - T Santhosh Kumar
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, United States
| | - Zhan-Guo Gao
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, United States
| | - Brian C Smith
- Department of Biochemistry and the Program in Chemical Biology, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, Wisconsin 53226, United States
| | - John A Auchampach
- Department of Pharmacology & Toxicology and the Cardiovascular Center, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, Wisconsin 53226, United States
| | - Kenneth A Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 9000 Rockville Pike, Bethesda, Maryland 20892, United States
| |
Collapse
|
4
|
Pawnikar S, Magenheimer BS, Joshi K, Munoz EN, Haldane A, Maser RL, Miao Y. Activation of Polycystin-1 Signaling by Binding of Stalk-derived Peptide Agonists. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.06.574465. [PMID: 38260358 PMCID: PMC10802338 DOI: 10.1101/2024.01.06.574465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
Polycystin-1 (PC1) is the membrane protein product of the PKD1 gene whose mutation is responsible for 85% of the cases of autosomal dominant polycystic kidney disease (ADPKD). ADPKD is primarily characterized by the formation of renal cysts and potential kidney failure. PC1 is an atypical G protein-coupled receptor (GPCR) consisting of 11 transmembrane helices and an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal (NTF) and membrane-embedded C-terminal (CTF) fragments. Recently, signaling activation of the PC1 CTF was shown to be regulated by a stalk tethered agonist (TA), a distinct mechanism observed in the adhesion GPCR family. A novel allosteric activation pathway was elucidated for the PC1 CTF through a combination of Gaussian accelerated molecular dynamics (GaMD), mutagenesis and cellular signaling experiments. Here, we show that synthetic, soluble peptides with 7 to 21 residues derived from the stalk TA, in particular, peptides including the first 9 residues (p9), 17 residues (p17) and 21 residues (p21) exhibited the ability to re-activate signaling by a stalkless PC1 CTF mutant in cellular assays. To reveal molecular mechanisms of stalk peptide-mediated signaling activation, we have applied a novel Peptide GaMD (Pep-GaMD) algorithm to elucidate binding conformations of selected stalk peptide agonists p9, p17 and p21 to the stalkless PC1 CTF. The simulations revealed multiple specific binding regions of the stalk peptide agonists to the PC1 protein including an "intermediate" bound yet inactive state. Our Pep-GaMD simulation findings were consistent with the cellular assay experimental data. Binding of peptide agonists to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of the PC1 CTF signaling activation mechanism. Using sequence covariation analysis of PC1 homologs, we further showed that the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. Therefore, structural dynamic insights into the mechanisms of PC1 activation by stalk-derived peptide agonists have enabled an in-depth understanding of PC1 signaling. They will form a foundation for development of PC1 as a therapeutic target for the treatment of ADPKD.
Collapse
Affiliation(s)
- Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Brenda S. Magenheimer
- Clinical Laboratory Sciences, University of Kansas Medical Center, Kansas City, KS 66160
- The Jared Grantham Kidney Institute, University of Kansas Medical Center, Kansas City, KS 66160
| | - Keya Joshi
- Department of Pharmacology and Computational Medicine Program, University of North Carolina – Chapel Hill, Chapel Hill, NC 27599
| | - Ericka Nevarez Munoz
- Clinical Laboratory Sciences, University of Kansas Medical Center, Kansas City, KS 66160
| | - Allan Haldane
- Dept of Physics, and Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA 19122
| | - Robin L. Maser
- Departments of Biochemistry and Molecular Biology, University of Kansas Medical Center, Kansas City, KS 66160
- Clinical Laboratory Sciences, University of Kansas Medical Center, Kansas City, KS 66160
- The Jared Grantham Kidney Institute, University of Kansas Medical Center, Kansas City, KS 66160
| | - Yinglong Miao
- Department of Pharmacology and Computational Medicine Program, University of North Carolina – Chapel Hill, Chapel Hill, NC 27599
| |
Collapse
|
5
|
Akhter S, Tang Z, Wang J, Haboro M, Holmstrom ED, Wang J, Miao Y. Mechanism of Ligand Binding to Theophylline RNA Aptamer. J Chem Inf Model 2024; 64:1017-1029. [PMID: 38226603 PMCID: PMC11058067 DOI: 10.1021/acs.jcim.3c01454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2024]
Abstract
Studying RNA-ligand interactions and quantifying their binding thermodynamics and kinetics are of particular relevance in the field of drug discovery. Here, we combined biochemical binding assays and accelerated molecular simulations to investigate ligand binding and dissociation in RNA using the theophylline-binding RNA as a model system. All-atom simulations using a Ligand Gaussian accelerated Molecular Dynamics method (LiGaMD) have captured repetitive binding and dissociation of theophylline and caffeine to RNA. Theophylline's binding free energy and kinetic rate constants align with our experimental data, while caffeine's binding affinity is over 10,000 times weaker, and its kinetics could not be determined. LiGaMD simulations allowed us to identify distinct low-energy conformations and multiple ligand binding pathways to RNA. Simulations revealed a "conformational selection" mechanism for ligand binding to the flexible RNA aptamer, which provides important mechanistic insights into ligand binding to the theophylline-binding model. Our findings suggest that compound docking using a structural ensemble of representative RNA conformations would be necessary for structure-based drug design of flexible RNA.
Collapse
Affiliation(s)
- Sana Akhter
- Computational Biology Program and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Zhichao Tang
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
| | - Jinan Wang
- Computational Biology Program and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Mercy Haboro
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
| | - Erik D Holmstrom
- Department of Molecular Biosciences and Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Jingxin Wang
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Computational Biology Program and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
6
|
Madhu MK, Shewani K, Murarka RK. Biased Signaling in Mutated Variants of β 2-Adrenergic Receptor: Insights from Molecular Dynamics Simulations. J Chem Inf Model 2024; 64:449-469. [PMID: 38194225 DOI: 10.1021/acs.jcim.3c01481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
The molecular basis of receptor bias in G protein-coupled receptors (GPCRs) caused by mutations that preferentially activate specific intracellular transducers over others remains poorly understood. Two experimentally identified biased variants of β2-adrenergic receptors (β2AR), a prototypical GPCR, are a triple mutant (T68F, Y132A, and Y219A) and a single mutant (Y219A); the former bias the receptor toward the β-arrestin pathway by disfavoring G protein engagement, while the latter induces G protein signaling explicitly due to selection against GPCR kinases (GRKs) that phosphorylate the receptor as a prerequisite of β-arrestin binding. Though rigorous characterizations have revealed functional implications of these mutations, the atomistic origin of the observed transducer selectivity is not clear. In this study, we investigated the allosteric mechanism of receptor bias in β2AR using microseconds of all-atom Gaussian accelerated molecular dynamics (GaMD) simulations. Our observations reveal distinct rearrangements in transmembrane helices, intracellular loop 3, and critical residues R1313.50 and Y3267.53 in the conserved motifs D(E)RY and NPxxY for the mutant receptors, leading to their specific transducer interactions. Moreover, partial dissociation of G protein from the receptor core is observed in the simulations of the triple mutant in contrast to the single mutant and wild-type receptor. The reorganization of allosteric communications from the extracellular agonist BI-167107 to the intracellular receptor-transducer interfaces drives the conformational rearrangements responsible for receptor bias in the single and triple mutants. The molecular insights into receptor bias of β2AR presented here could improve the understanding of biased signaling in GPCRs, potentially opening new avenues for designing novel therapeutics with fewer side-effects and superior efficacy.
Collapse
Affiliation(s)
- Midhun K Madhu
- Department of Biological Sciences, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| | - Kunal Shewani
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| | - Rajesh K Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| |
Collapse
|
7
|
Koirala K, Joshi K, Adediwura V, Wang J, Do H, Miao Y. Accelerating Molecular Dynamics Simulations for Drug Discovery. Methods Mol Biol 2024; 2714:187-202. [PMID: 37676600 DOI: 10.1007/978-1-0716-3441-7_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
Accurate prediction of ligand binding thermodynamics and kinetics is crucial in drug design. However, it remains challenging for conventional molecular dynamics (MD) simulations due to sampling issues. Gaussian accelerated MD (GaMD) is an enhanced sampling method that adds a harmonic boost to overcome energy barriers, which has demonstrated significant benefits in exploring protein-ligand interactions. Especially, the ligand GaMD (LiGaMD) applies a selective boost potential to the ligand nonbonded potential energy, significantly improving sampling for ligand binding and dissociation. Furthermore, a selective boost potential is applied to the potential of both ligand and protein residues around binding pocket in LiGaMD2 to further increase the sampling of protein-ligand interaction. LiGaMD and LiGaMD2 simulations could capture repetitive ligand binding and unbinding events within microsecond simulations, allowing to simultaneously characterize ligand binding thermodynamics and kinetics, which is expected to greatly facilitate drug design. In this chapter, we provide a brief review of the status of LiGaMD in drug discovery and outline its usage.
Collapse
Affiliation(s)
- Kushal Koirala
- Computational Biology Program and Department of Molecular Biosciences, The University of Kansas, Lawrence, KS, USA
| | - Keya Joshi
- Computational Biology Program and Department of Molecular Biosciences, The University of Kansas, Lawrence, KS, USA
| | - Victor Adediwura
- Computational Biology Program and Department of Molecular Biosciences, The University of Kansas, Lawrence, KS, USA
| | - Jinan Wang
- Computational Biology Program and Department of Molecular Biosciences, The University of Kansas, Lawrence, KS, USA
| | - Hung Do
- Computational Biology Program and Department of Molecular Biosciences, The University of Kansas, Lawrence, KS, USA
| | - Yinglong Miao
- Computational Biology Program and Department of Molecular Biosciences, The University of Kansas, Lawrence, KS, USA.
| |
Collapse
|
8
|
Adediwura VA, Miao Y. Mechanistic Insights into Peptide Binding and Deactivation of an Adhesion G Protein-Coupled Receptor. Molecules 2023; 29:164. [PMID: 38202747 PMCID: PMC10780249 DOI: 10.3390/molecules29010164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2023] [Revised: 12/20/2023] [Accepted: 12/23/2023] [Indexed: 01/12/2024] Open
Abstract
Adhesion G protein-coupled receptors (ADGRGs) play critical roles in the reproductive, neurological, cardiovascular, and endocrine systems. In particular, ADGRG2 plays a significant role in Ewing sarcoma cell proliferation, parathyroid cell function, and male fertility. In 2022, a cryo-EM structure was reported for the active ADGRG2 bound by an optimized peptide agonist IP15 and the Gs protein. The IP15 peptide agonist was also modified to antagonists 4PH-E and 4PH-D with mutations of the 4PH residue to Glu and Asp, respectively. However, experimental structures of inactive antagonist-bound ADGRs remain to be resolved, and the activation mechanism of ADGRs such as ADGRG2 is poorly understood. Here, we applied Gaussian accelerated molecular dynamics (GaMD) simulations to probe conformational dynamics of the agonist- and antagonist-bound ADGRG2. By performing GaMD simulations, we were able to identify important low-energy conformations of ADGRG2 in the active, intermediate, and inactive states, as well as explore the binding conformations of each peptide. Moreover, our simulations revealed critical peptide-receptor residue interactions during the deactivation of ADGRG2. In conclusion, through GaMD simulations, we uncovered mechanistic insights into peptide (agonist and antagonist) binding and deactivation of the ADGRG2. These findings will potentially facilitate rational design of new peptide modulators of ADGRG2 and other ADGRs.
Collapse
Affiliation(s)
| | - Yinglong Miao
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA;
| |
Collapse
|
9
|
Perfilova VN, Muzyko EA, Taran AS, Shevchenko AA, Naumenko LV. Problems and prospects for finding new pharmacological agents among adenosine receptor agonists, antagonists, or their allosteric modulators for the treatment of cardiovascular diseases. BIOMEDITSINSKAIA KHIMIIA 2023; 69:353-370. [PMID: 38153051 DOI: 10.18097/pbmc20236906353] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2023]
Abstract
A1-adenosine receptors (A1AR) are widely distributed in the human body and mediate many different effects. They are abundantly present in the cardiovascular system, where they control angiogenesis, vascular tone, heart rate, and conduction. This makes the cardiovascular system A1AR an attractive target for the treatment of cardiovascular diseases (CVD). The review summarizes the literature data on the structure and functioning of A1AR, and analyzes their involvement in the formation of myocardial hypertrophy, ischemia-reperfusion damage, various types of heart rhythm disorders, chronic heart failure, and arterial hypertension. Special attention is paid to the role of some allosteric regulators of A1AR as potential agents for the CVD treatment.
Collapse
Affiliation(s)
- V N Perfilova
- Volgograd State Medical University, Volgograd, Russia; Volgograd Medical Research Center, Volgograd, Russia
| | - E A Muzyko
- Volgograd State Medical University, Volgograd, Russia
| | - A S Taran
- Volgograd State Medical University, Volgograd, Russia
| | | | - L V Naumenko
- Volgograd State Medical University, Volgograd, Russia
| |
Collapse
|
10
|
Do HN, Wang J, Miao Y. Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors. JACS AU 2023; 3:3165-3180. [PMID: 38034960 PMCID: PMC10685416 DOI: 10.1021/jacsau.3c00503] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/26/2023] [Revised: 10/10/2023] [Accepted: 10/18/2023] [Indexed: 12/02/2023]
Abstract
G-protein-coupled receptors (GPCRs) make up the largest superfamily of human membrane proteins and represent primary targets of ∼1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon the binding of positive and negative allosteric modulators (PAMs and NAMs). The mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy prOfiling Workflow (GLOW). GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence and absence of the modulator. DL and free energy calculations revealed significantly reduced dynamic fluctuations and conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G-protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes. Therefore, GPCR allostery exhibits a dynamic "conformational selection" mechanism. In the absence of available modulator-bound structures as for most current GPCRs, it is critical to use a structural ensemble of representative GPCR conformations rather than a single structure for compound docking ("ensemble docking"), which will potentially improve structure-based design of novel allosteric drugs of GPCRs.
Collapse
Affiliation(s)
| | - Jinan Wang
- Computational Biology Program
and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | | |
Collapse
|
11
|
Maria-Solano MA, Choi S. Dynamic allosteric networks drive adenosine A 1 receptor activation and G-protein coupling. eLife 2023; 12:RP90773. [PMID: 37656635 PMCID: PMC10473838 DOI: 10.7554/elife.90773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/03/2023] Open
Abstract
G-protein coupled receptors (GPCRs) present specific activation pathways and signaling among receptor subtypes. Hence, an extensive knowledge of the structural dynamics of the receptor is critical for the development of therapeutics. Here, we target the adenosine A1 receptor (A1R), for which a negligible number of drugs have been approved. We combine molecular dynamics simulations, enhanced sampling techniques, network theory and pocket detection to decipher the activation pathway of A1R, decode the allosteric networks and identify transient pockets. The A1R activation pathway reveal hidden intermediate and pre-active states together with the inactive and fully-active states observed experimentally. The protein energy networks computed throughout these conformational states successfully unravel the extra and intracellular allosteric centers and the communication pathways that couples them. We observe that the allosteric networks are dynamic, being increased along activation and fine-tuned in presence of the trimeric G-proteins. Overlap of transient pockets and energy networks uncover how the allosteric coupling between pockets and distinct functional regions of the receptor is altered along activation. By an in-depth analysis of the bridge between activation pathway, energy networks and transient pockets, we provide a further understanding of A1R. This information can be useful to ease the design of allosteric modulators for A1R.
Collapse
Affiliation(s)
- Miguel A Maria-Solano
- Global AI Drug Discovery Center, College of Pharmacy and Graduate School of Pharmaceutical Science, Ewha Womans UniversitySeoulRepublic of Korea
| | - Sun Choi
- Global AI Drug Discovery Center, College of Pharmacy and Graduate School of Pharmaceutical Science, Ewha Womans UniversitySeoulRepublic of Korea
| |
Collapse
|
12
|
Nguyen ATN, Tran QL, Baltos JA, McNeill SM, Nguyen DTN, May LT. Small molecule allosteric modulation of the adenosine A 1 receptor. Front Endocrinol (Lausanne) 2023; 14:1184360. [PMID: 37435481 PMCID: PMC10331460 DOI: 10.3389/fendo.2023.1184360] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/11/2023] [Accepted: 05/23/2023] [Indexed: 07/13/2023] Open
Abstract
G protein-coupled receptors (GPCRs) represent the target for approximately a third of FDA-approved small molecule drugs. The adenosine A1 receptor (A1R), one of four adenosine GPCR subtypes, has important (patho)physiological roles in humans. A1R has well-established roles in the regulation of the cardiovascular and nervous systems, where it has been identified as a potential therapeutic target for a number of conditions, including cardiac ischemia-reperfusion injury, cognition, epilepsy, and neuropathic pain. A1R small molecule drugs, typically orthosteric ligands, have undergone clinical trials. To date, none have progressed into the clinic, predominantly due to dose-limiting unwanted effects. The development of A1R allosteric modulators that target a topographically distinct binding site represent a promising approach to overcome current limitations. Pharmacological parameters of allosteric ligands, including affinity, efficacy and cooperativity, can be optimized to regulate A1R activity with high subtype, spatial and temporal selectivity. This review aims to offer insights into the A1R as a potential therapeutic target and highlight recent advances in the structural understanding of A1R allosteric modulation.
Collapse
Affiliation(s)
- Anh T. N. Nguyen
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, Australia
| | - Quan L. Tran
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, Australia
| | - Jo-Anne Baltos
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, Australia
| | - Samantha M. McNeill
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, Australia
| | - Diep T. N. Nguyen
- Department of Information Technology, Faculty of Engineering and Technology, Vietnam National University, Hanoi, Vietnam
| | - Lauren T. May
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, Australia
| |
Collapse
|
13
|
Do H, Wang J, Miao Y. Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors. RESEARCH SQUARE 2023:rs.3.rs-2543463. [PMID: 36865316 PMCID: PMC9980202 DOI: 10.21203/rs.3.rs-2543463/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/22/2023]
Abstract
G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~ 1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.
Collapse
|
14
|
Do HN, Wang J, Miao Y. Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.15.524128. [PMID: 36711515 PMCID: PMC9882226 DOI: 10.1101/2023.01.15.524128] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.
Collapse
Affiliation(s)
- Hung N. Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047
| |
Collapse
|
15
|
Xiaoli A, Yuzhen N, Qiong Y, Yang L, Yao X, Bing Z. Investigating the Dynamic Binding Behavior of PMX53 Cooperating with Allosteric Antagonist NDT9513727 to C5a Anaphylatoxin Chemotactic Receptor 1 through Gaussian Accelerated Molecular Dynamics and Free-Energy Perturbation Simulations. ACS Chem Neurosci 2022; 13:3502-3511. [PMID: 36428153 DOI: 10.1021/acschemneuro.2c00556] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
Abstract
C5a anaphylatoxin chemotactic receptor 1 (C5aR1) is an important target in anti-inflammatory therapeutics. The cyclic peptide antagonist PMX53 binds to the orthosteric site located in the extracellular vestibule of C5aR1, and the non-peptide antagonist NDT9513727 binds to the allosteric site formed by the middle region of TM3 (trans-membrane helix), TM4, and TM5. We catch a sight of the variational binding mode of PMX53 during the Gaussian accelerated molecular dynamic (GaMD) simulations. In the binary complex of C5aR1 and PMX53, the PMX53 takes a dynamic binding mechanism during the simulation. Namely, the side chain of Arg6 of PMX53 extends to TM6-TM7 (pose 1) or swings to TM5 (pose 2), forming a salt bridge with Glu199. Meanwhile, in the ternary complex of C5aR1 with PMX53 and NDT9513727, the side chain of Arg6 of PMX53 swings to TM5 (pose 2) from extending to TM6-TM7 (pose 1) at the beginning of the GaMD simulation. In subsequent simulation, PMX53 stabilizes in the pose 2 binding mode by forming a stable salt bridge with Glu199. The free-energy perturbation (FEP) calculations demonstrate that pose 1 (ΔGbinding = -10.94 kcal/mol) is more stable in the binary complex and pose 2 (ΔGbinding = -7.91 kcal/mol) is unstable because of highly dynamic TM5. NDT9513727 interacts directly with TM4 and TM5 and stabilizes the hydrophobic stack between the extracellular sides of the two helices. Therefore, pose 2 (ΔGbinding = -16.27 kcal/mol) is notably stable than pose 1 (ΔGbinding = -9.78 kcal/mol) in the ternary complex. The identification of a novel binding mode of PMX53 and the detailed structural information of PMX53 interacting with a receptor obtained by GaMD simulations will be helpful in designing potent antagonists of C5aR1.
Collapse
Affiliation(s)
- An Xiaoli
- Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China
| | - Niu Yuzhen
- Shandong Laboratory of Yantai Advanced Materials and Green Manufacturing, Yantai 264006, China.,Yantai Zhongke Research Institute of Advanced Materials and Green Chemical Engineering, Yantai 264006, China
| | - Yang Qiong
- Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China
| | - Lei Yang
- Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China
| | - Xiaojun Yao
- Dr. Neher's Biophysics Laboratory for Innovative Drug Discovery, State Key Laboratory of Quality Research in Chinese Medicine, Macau Institute for Applied Research in Medicine and Health, Macau University of Science and Technology, Taipa, Macau 999078, China
| | - Zhitong Bing
- Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China.,Advanced Energy Science and Technology Guangdong Laboratory, Huizhou 516000, China
| |
Collapse
|
16
|
The full activation mechanism of the adenosine A 1 receptor revealed by GaMD and Su-GaMD simulations. Proc Natl Acad Sci U S A 2022; 119:e2203702119. [PMID: 36215480 PMCID: PMC9586258 DOI: 10.1073/pnas.2203702119] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The full activation process of G protein-coupled receptor (GPCR) plays an important role in cellular signal transduction. However, it remains challenging to simulate the whole process in which the GPCR is recognized and activated by a ligand and then couples to the G protein on a reasonable simulation timescale. Here, we developed a molecular dynamics (MD) approach named supervised (Su) Gaussian accelerated MD (GaMD) by incorporating a tabu-like supervision algorithm into a standard GaMD simulation. By using this Su-GaMD method, from the active and inactive structure of adenosine A1 receptor (A1R), we successfully revealed the full activation mechanism of A1R, including adenosine (Ado)-A1R recognition, preactivation of A1R, and A1R-G protein recognition, in hundreds of nanoseconds of simulations. The binding of Ado to the extracellular side of A1R initiates conformational changes and the preactivation of A1R. In turn, the binding of Gi2 to the intracellular side of A1R causes a decrease in the volume of the extracellular orthosteric site and stabilizes the binding of Ado to A1R. Su-GaMD could be a useful tool to reconstruct or even predict ligand-protein and protein-protein recognition pathways on a short timescale. The intermediate states revealed in this study could provide more detailed complementary structural characterizations to facilitate the drug design of A1R in the future.
Collapse
|
17
|
IJzerman AP, Jacobson KA, Müller CE, Cronstein BN, Cunha RA. International Union of Basic and Clinical Pharmacology. CXII: Adenosine Receptors: A Further Update. Pharmacol Rev 2022; 74:340-372. [PMID: 35302044 PMCID: PMC8973513 DOI: 10.1124/pharmrev.121.000445] [Citation(s) in RCA: 74] [Impact Index Per Article: 37.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Our previous International Union of Basic and Clinical Pharmacology report on the nomenclature and classification of adenosine receptors (2011) contained a number of emerging developments with respect to this G protein-coupled receptor subfamily, including protein structure, protein oligomerization, protein diversity, and allosteric modulation by small molecules. Since then, a wealth of new data and results has been added, allowing us to explore novel concepts such as target binding kinetics and biased signaling of adenosine receptors, to examine a multitude of receptor structures and novel ligands, to gauge new pharmacology, and to evaluate clinical trials with adenosine receptor ligands. This review should therefore be considered a further update of our previous reports from 2001 and 2011. SIGNIFICANCE STATEMENT: Adenosine receptors (ARs) are of continuing interest for future treatment of chronic and acute disease conditions, including inflammatory diseases, neurodegenerative afflictions, and cancer. The design of AR agonists ("biased" or not) and antagonists is largely structure based now, thanks to the tremendous progress in AR structural biology. The A2A- and A2BAR appear to modulate the immune response in tumor biology. Many clinical trials for this indication are ongoing, whereas an A2AAR antagonist (istradefylline) has been approved as an anti-Parkinson agent.
Collapse
Affiliation(s)
- Adriaan P IJzerman
- Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands (A.P.IJ.); National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases, Molecular Recognition Section, Bethesda, Maryland (K.A.J.); Universität Bonn, Bonn, Germany (C.E.M.); New York University School of Medicine, New York, New York (B.N.C.); and Center for Neurosciences and Cell Biology and Faculty of Medicine, University of Coimbra, Coimbra, Portugal (R.A.C.)
| | - Kenneth A Jacobson
- Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands (A.P.IJ.); National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases, Molecular Recognition Section, Bethesda, Maryland (K.A.J.); Universität Bonn, Bonn, Germany (C.E.M.); New York University School of Medicine, New York, New York (B.N.C.); and Center for Neurosciences and Cell Biology and Faculty of Medicine, University of Coimbra, Coimbra, Portugal (R.A.C.)
| | - Christa E Müller
- Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands (A.P.IJ.); National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases, Molecular Recognition Section, Bethesda, Maryland (K.A.J.); Universität Bonn, Bonn, Germany (C.E.M.); New York University School of Medicine, New York, New York (B.N.C.); and Center for Neurosciences and Cell Biology and Faculty of Medicine, University of Coimbra, Coimbra, Portugal (R.A.C.)
| | - Bruce N Cronstein
- Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands (A.P.IJ.); National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases, Molecular Recognition Section, Bethesda, Maryland (K.A.J.); Universität Bonn, Bonn, Germany (C.E.M.); New York University School of Medicine, New York, New York (B.N.C.); and Center for Neurosciences and Cell Biology and Faculty of Medicine, University of Coimbra, Coimbra, Portugal (R.A.C.)
| | - Rodrigo A Cunha
- Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands (A.P.IJ.); National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases, Molecular Recognition Section, Bethesda, Maryland (K.A.J.); Universität Bonn, Bonn, Germany (C.E.M.); New York University School of Medicine, New York, New York (B.N.C.); and Center for Neurosciences and Cell Biology and Faculty of Medicine, University of Coimbra, Coimbra, Portugal (R.A.C.)
| |
Collapse
|
18
|
Wang J, Bhattarai A, Do HN, Akhter S, Miao Y. Molecular Simulations and Drug Discovery of Adenosine Receptors. Molecules 2022; 27:2054. [PMID: 35408454 PMCID: PMC9000248 DOI: 10.3390/molecules27072054] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Revised: 03/18/2022] [Accepted: 03/20/2022] [Indexed: 02/02/2023] Open
Abstract
G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.
Collapse
Affiliation(s)
| | | | | | | | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA; (J.W.); (A.B.); (H.N.D.); (S.A.)
| |
Collapse
|
19
|
Pawnikar S, Miao Y. Mechanism of Peptide Agonist Binding in CXCR4 Chemokine Receptor. Front Mol Biosci 2022; 9:821055. [PMID: 35359589 PMCID: PMC8963245 DOI: 10.3389/fmolb.2022.821055] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 02/14/2022] [Indexed: 01/07/2023] Open
Abstract
Chemokine receptors are key G-protein-coupled receptors (GPCRs) that control cell migration in immune system responses, development of cardiovascular and central nervous systems, and numerous diseases. In particular, the CXCR4 chemokine receptor promotes metastasis, tumor growth and angiogenesis in cancers. CXCR4 is also used as one of the two co-receptors for T-tropic HIV-1 entry into host cells. Therefore, CXCR4 serves as an important therapeutic target for treating cancers and HIV infection. Apart from the CXCL12 endogenous peptide agonist, previous studies suggested that the first 17 amino acids of CXCL12 are sufficient to activate CXCR4. Two 17-residue peptides with positions 1-4 mutated to RSVM and ASLW functioned as super and partial agonists of CXCR4, respectively. However, the mechanism of peptide agonist binding in CXCR4 remains unclear. Here, we have investigated this mechanism through all-atom simulations using a novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method. The Pep-GaMD simulations have allowed us to explore representative binding conformations of each peptide and identify critical low-energy states of CXCR4 activated by the super versus partial peptide agonists. Our simulations have provided important mechanistic insights into peptide agonist binding in CXCR4, which are expected to facilitate rational design of new peptide modulators of CXCR4 and other chemokine receptors.
Collapse
|
20
|
Egyed A, Kiss DJ, Keserű GM. The Impact of the Secondary Binding Pocket on the Pharmacology of Class A GPCRs. Front Pharmacol 2022; 13:847788. [PMID: 35355719 PMCID: PMC8959758 DOI: 10.3389/fphar.2022.847788] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2022] [Accepted: 02/01/2022] [Indexed: 12/19/2022] Open
Abstract
G-protein coupled receptors (GPCRs) are considered important therapeutic targets due to their pathophysiological significance and pharmacological relevance. Class A receptors represent the largest group of GPCRs that gives the highest number of validated drug targets. Endogenous ligands bind to the orthosteric binding pocket (OBP) embedded in the intrahelical space of the receptor. During the last 10 years, however, it has been turned out that in many receptors there is secondary binding pocket (SBP) located in the extracellular vestibule that is much less conserved. In some cases, it serves as a stable allosteric site harbouring allosteric ligands that modulate the pharmacology of orthosteric binders. In other cases it is used by bitopic compounds occupying both the OBP and SBP. In these terms, SBP binding moieties might influence the pharmacology of the bitopic ligands. Together with others, our research group showed that SBP binders contribute significantly to the affinity, selectivity, functional activity, functional selectivity and binding kinetics of bitopic ligands. Based on these observations we developed a structure-based protocol for designing bitopic compounds with desired pharmacological profile.
Collapse
Affiliation(s)
| | | | - György M. Keserű
- Medicinal Chemistry Research Group, Research Centre for Natural Sciences, Budapest, Hungary
| |
Collapse
|
21
|
Do HN, Wang J, Bhattarai A, Miao Y. GLOW: A Workflow Integrating Gaussian-Accelerated Molecular Dynamics and Deep Learning for Free Energy Profiling. J Chem Theory Comput 2022; 18:1423-1436. [PMID: 35200019 PMCID: PMC9773012 DOI: 10.1021/acs.jctc.1c01055] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We introduce a Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and free energy profiling workflow (GLOW) to predict molecular determinants and map free energy landscapes of biomolecules. All-atom GaMD-enhanced sampling simulations are first performed on biomolecules of interest. Structural contact maps are then calculated from GaMD simulation frames and transformed into images for building DL models using a convolutional neural network. Important structural contacts are further determined from DL models of attention maps of the structural contact gradients, which allow us to identify the system reaction coordinates. Finally, free energy profiles are calculated for the selected reaction coordinates through energetic reweighting of the GaMD simulations. We have also successfully demonstrated GLOW for the characterization of activation and allosteric modulation of a G protein-coupled receptor, using the adenosine A1 receptor (A1AR) as a model system. GLOW findings are highly consistent with previous experimental and computational studies of the A1AR, while also providing further mechanistic insights into the receptor function. In summary, GLOW provides a systematic approach to mapping free energy landscapes of biomolecules. The GLOW workflow and its user manual can be downloaded at http://miaolab.org/GLOW.
Collapse
Affiliation(s)
- Hung N. Do
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Jinan Wang
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Apurba Bhattarai
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Yinglong Miao
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047,Corresponding author:
| |
Collapse
|
22
|
Madhu MK, Debroy A, Murarka RK. Molecular Insights into Phosphorylation-Induced Allosteric Conformational Changes in a β 2-Adrenergic Receptor. J Phys Chem B 2022; 126:1917-1932. [PMID: 35196859 DOI: 10.1021/acs.jpcb.1c08610] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The large conformational flexibility of G protein-coupled receptors (GPCRs) has been a puzzle in structural and pharmacological studies for the past few decades. Apart from structural rearrangements induced by ligands, enzymatic phosphorylations by GPCR kinases (GRKs) at the carboxy-terminal tail (C-tail) of a GPCR also make conformational alterations to the transmembrane helices and facilitates the binding of one of its transducer proteins named β-arrestin. The phosphorylation-induced conformational transition of the receptor that causes specific binding to β-arrestin but prevents the association of other transducers such as G proteins lacks atomistic understanding and is elusive to experimental studies. Using microseconds of all-atom conventional and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigate the allosteric mechanism of phosphorylation induced-conformational changes in β2-adrenergic receptor, a well-characterized GPCR model system. Free energy profiles reveal that the phosphorylated receptor samples a new conformational state in addition to the canonical active state corroborating with recent nuclear magnetic resonance experimental findings. The new state has a smaller intracellular cavity that is likely to accommodate β-arrestin better than G protein. Using contact map and inter-residue interaction energy calculations, we found the phosphorylated C-tail adheres to the cytosolic surface of the transmembrane domain of the receptor. Transfer entropy calculations show that the C-tail residues drive the correlated motions of TM residues, and the allosteric signal is relayed via several residues at the cytosolic surface. Our results also illustrate how the redistribution of inter-residue nonbonding interaction couples with the allosteric communication from the phosphorylated C-tail to the transmembrane. Atomistic insight into phosphorylation-induced β-arrestin specific conformation is therapeutically important to design drugs with higher efficacy and fewer side effects. Our results, therefore, open novel opportunities to fine-tune β-arrestin bias in GPCR signaling.
Collapse
Affiliation(s)
- Midhun K Madhu
- Department of Biological Sciences, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| | - Annesha Debroy
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| | - Rajesh K Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| |
Collapse
|
23
|
Pawnikar S, Bhattarai A, Wang J, Miao Y. Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives. Adv Appl Bioinform Chem 2022; 15:1-19. [PMID: 35023931 PMCID: PMC8747661 DOI: 10.2147/aabc.s247950] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2020] [Accepted: 02/20/2021] [Indexed: 12/12/2022] Open
Abstract
Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.
Collapse
Affiliation(s)
- Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| |
Collapse
|
24
|
Challenges and frontiers of computational modelling of biomolecular recognition. QRB DISCOVERY 2022. [DOI: 10.1017/qrd.2022.11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
Abstract
Biomolecular recognition including binding of small molecules, peptides and proteins to their target receptors plays a key role in cellular function and has been targeted for therapeutic drug design. However, the high flexibility of biomolecules and slow binding and dissociation processes have presented challenges for computational modelling. Here, we review the challenges and computational approaches developed to characterise biomolecular binding, including molecular docking, molecular dynamics simulations (especially enhanced sampling) and machine learning. Further improvements are still needed in order to accurately and efficiently characterise binding structures, mechanisms, thermodynamics and kinetics of biomolecules in the future.
Collapse
|
25
|
Rehman AU, Lu S, Khan AA, Khurshid B, Rasheed S, Wadood A, Zhang J. Hidden allosteric sites and De-Novo drug design. Expert Opin Drug Discov 2021; 17:283-295. [PMID: 34933653 DOI: 10.1080/17460441.2022.2017876] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
INTRODUCTION Hidden allosteric sites are not visible in apo-crystal structures, but they may be visible in holo-structures when a certain ligand binds and maintains the ligand intended conformation. Several computational and experimental techniques have been used to investigate these hidden sites but identifying them remains a challenge. AREAS COVERED This review provides a summary of the many theoretical approaches for predicting hidden allosteric sites in disease-related proteins. Furthermore, promising cases have been thoroughly examined to reveal the hidden allosteric site and its modulator. EXPERT OPINION In the recent past, with the development in scientific techniques and bioinformatics tools, the number of drug targets for complex human diseases has significantly increased but unfortunately most of these targets are undruggable due to several reasons. Alternative strategies such as finding cryptic (hidden) allosteric sites are an attractive approach for exploitation of the discovery of new targets. These hidden sites are difficult to recognize compared to allosteric sites, mainly due to a lack of visibility in the crystal structure. In our opinion, after many years of development, MD simulations are finally becoming successful for obtaining a detailed molecular description of drug-target interaction.
Collapse
Affiliation(s)
- Ashfaq Ur Rehman
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China
| | - Shaoyong Lu
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China
| | - Abdul Aziz Khan
- Bio-X Institutes, Key Laboratory for the Genetics of Development and Neuropsychiatric Disorders (Ministry of Education), Shanghai Key Laboratory of Psychotic Disorders, and Brain Science and Technology Research Center, Institute of Psychology and Behavioral Sciences, Shanghai Jiao Tong University, Shanghai, China
| | - Beenish Khurshid
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Salman Rasheed
- National Center for Bioinformatics, Quaid-e-Azam University, Islamabad, Pakistan
| | - Abdul Wadood
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Jian Zhang
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China.,School of Pharmaceutical Sciences, Zhengzhou University, Zhengzhou, China
| |
Collapse
|
26
|
Wang J, Lan L, Wu X, Xu L, Miao Y. Mechanism of RNA recognition by a Musashi RNA-binding protein. Curr Res Struct Biol 2021; 4:10-20. [PMID: 34988468 PMCID: PMC8695263 DOI: 10.1016/j.crstbi.2021.12.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Revised: 10/31/2021] [Accepted: 12/07/2021] [Indexed: 12/21/2022] Open
Abstract
The Musashi RNA-binding proteins (RBPs) regulate translation of target mRNAs and maintenance of cell stemness and tumorigenesis. Musashi-1 (MSI1), long considered as an intestinal and neural stem cell marker, has been more recently found to be over expressed in many cancers. It has served as an important drug target for treating acute myeloid leukemia and solid tumors such as ovarian, colorectal and bladder cancer. One of the reported binding targets of MSI1 is Numb, a negative regulator of the Notch signaling. However, the dynamic mechanism of Numb RNA binding to MSI1 remains unknown, largely hindering effective drug design targeting this critical interaction. Here, we have performed extensive all-atom microsecond-timescale simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method, which successfully captured multiple times of spontaneous and highly accurate binding of the Numb RNA from bulk solvent to the MSI1 protein target site. GaMD simulations revealed that Numb RNA binding to MSI1 involved largely induced fit in both the RNA and protein. The simulations also identified important low-energy intermediate conformational states during RNA binding, in which Numb interacted mainly with the β2-β3 loop and C terminus of MSI1. The mechanistic understanding of RNA binding obtained from our GaMD simulations is expected to facilitate rational structure-based drug design targeting MSI1 and other RBPs.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology, University of Kansas, Lawrence, KS, 66047, USA
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Lan Lan
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Xiaoqing Wu
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Liang Xu
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
- Department of Radiation Oncology, The University of Kansas Cancer Center, Kansas City, KS, 66160, USA
| | - Yinglong Miao
- Center for Computational Biology, University of Kansas, Lawrence, KS, 66047, USA
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| |
Collapse
|
27
|
Saini A, Patel R, Gaba S, Singh G, Gupta GD, Monga V. Adenosine receptor antagonists: Recent advances and therapeutic perspective. Eur J Med Chem 2021; 227:113907. [PMID: 34695776 DOI: 10.1016/j.ejmech.2021.113907] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Revised: 10/05/2021] [Accepted: 10/05/2021] [Indexed: 12/14/2022]
Abstract
Adenosine is an endogenous purine-based nucleoside expressed nearly in all body tissues. It regulates various body functions by activating four G-protein coupled receptors, A1, A2A, A2B, and A3. These receptors are widely acknowledged as drug targets for treating different neurological, metabolic, and inflammatory diseases. Although numerous adenosine receptor inhibitors have been developed worldwide, achieving target selectivity is still a big hurdle in drug development. However, the identification of specific radioligands-based affinity assay, fluorescent ligands, and MS-based ligand assay have contributed to the development of selective and potent adenosine ligands. In recent years various small heterocyclic-based molecules have shown some promising results. Istradefylline has been approved for treating Parkinson's in Japan, while preladenant, tozadenant, CVT-6883, MRS-1523, and many more are under different phases of clinical development. The present review is focused on the quest to develop potent and selective adenosine inhibitors from 2013 to early 2021 by various research groups. The review also highlights their biological activity, selectivity, structure-activity relationship, molecular docking, and mechanistic studies. A special emphsesis on drug designing strategies has been also given the manuscript. The comprehensive compilation of research work carried out in the field will provide inevitable scope for designing and developing novel adenosine inhibitors with improved selectivity and efficacy.
Collapse
Affiliation(s)
- Anjali Saini
- Department of Pharmaceutical Chemistry, ISF College of Pharmacy, GT Road, Ghal Kalan, Moga, 142001, Punjab, India
| | - Rajiv Patel
- Department of Pharmaceutical Chemistry, ISF College of Pharmacy, GT Road, Ghal Kalan, Moga, 142001, Punjab, India
| | - Sobhi Gaba
- Department of Pharmaceutical Chemistry, ISF College of Pharmacy, GT Road, Ghal Kalan, Moga, 142001, Punjab, India
| | - Gurpreet Singh
- Department of Pharmaceutical Chemistry, ISF College of Pharmacy, GT Road, Ghal Kalan, Moga, 142001, Punjab, India.
| | - G D Gupta
- Department of Pharmaceutics, ISF College of Pharmacy, GT Road, Ghal Kalan, Moga, 142001, Punjab, India
| | - Vikramdeep Monga
- Department of Pharmaceutical Chemistry, ISF College of Pharmacy, GT Road, Ghal Kalan, Moga, 142001, Punjab, India.
| |
Collapse
|
28
|
Wang J, Arantes PR, Bhattarai A, Hsu RV, Pawnikar S, Huang YMM, Palermo G, Miao Y. Gaussian accelerated molecular dynamics (GaMD): principles and applications. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2021; 11:e1521. [PMID: 34899998 PMCID: PMC8658739 DOI: 10.1002/wcms.1521] [Citation(s) in RCA: 95] [Impact Index Per Article: 31.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2020] [Accepted: 01/28/2021] [Indexed: 12/20/2022]
Abstract
Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, KS, 66047, United States
| | - Pablo R Arantes
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr, Lawrence, KS, 66047, United States
| | - Rohaine V Hsu
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, KS, 66047, United States
| | - Yu-Ming M Huang
- Department of Physics & Astronomy, Wayne State University, 666 W Hancock St, Detroit, MI 48207, USA
| | - Giulia Palermo
- Department of Bioengineering and Department of Chemistry, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, Kansas 66047, United States
| |
Collapse
|
29
|
Draper-Joyce CJ, Bhola R, Wang J, Bhattarai A, Nguyen ATN, Cowie-Kent I, O'Sullivan K, Chia LY, Venugopal H, Valant C, Thal DM, Wootten D, Panel N, Carlsson J, Christie MJ, White PJ, Scammells P, May LT, Sexton PM, Danev R, Miao Y, Glukhova A, Imlach WL, Christopoulos A. Positive allosteric mechanisms of adenosine A 1 receptor-mediated analgesia. Nature 2021; 597:571-576. [PMID: 34497422 PMCID: PMC8711093 DOI: 10.1038/s41586-021-03897-2] [Citation(s) in RCA: 77] [Impact Index Per Article: 25.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Accepted: 08/11/2021] [Indexed: 02/08/2023]
Abstract
The adenosine A1 receptor (A1R) is a promising therapeutic target for non-opioid analgesic agents to treat neuropathic pain1,2. However, development of analgesic orthosteric A1R agonists has failed because of a lack of sufficient on-target selectivity as well as off-tissue adverse effects3. Here we show that [2-amino-4-(3,5-bis(trifluoromethyl)phenyl)thiophen-3-yl)(4-chlorophenyl)methanone] (MIPS521), a positive allosteric modulator of the A1R, exhibits analgesic efficacy in rats in vivo through modulation of the increased levels of endogenous adenosine that occur in the spinal cord of rats with neuropathic pain. We also report the structure of the A1R co-bound to adenosine, MIPS521 and a Gi2 heterotrimer, revealing an extrahelical lipid-detergent-facing allosteric binding pocket that involves transmembrane helixes 1, 6 and 7. Molecular dynamics simulations and ligand kinetic binding experiments support a mechanism whereby MIPS521 stabilizes the adenosine-receptor-G protein complex. This study provides proof of concept for structure-based allosteric drug design of non-opioid analgesic agents that are specific to disease contexts.
Collapse
Affiliation(s)
- Christopher J Draper-Joyce
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
- The Florey Institute of Neuroscience and Mental Health, University of Melbourne, Parkville, Victoria, Australia
| | - Rebecca Bhola
- Department of Physiology, Monash Biomedicine Discovery Institute, Monash University, Clayton, Victoria, Australia
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, USA
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, USA
| | - Anh T N Nguyen
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - India Cowie-Kent
- Department of Physiology, Monash Biomedicine Discovery Institute, Monash University, Clayton, Victoria, Australia
| | - Kelly O'Sullivan
- Department of Physiology, Monash Biomedicine Discovery Institute, Monash University, Clayton, Victoria, Australia
| | - Ling Yeong Chia
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Hariprasad Venugopal
- Department of Biochemistry and Molecular Biology, Monash University, Clayton, Victoria, Australia
| | - Celine Valant
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - David M Thal
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Denise Wootten
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Nicolas Panel
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Jens Carlsson
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Macdonald J Christie
- Discipline of Pharmacology, Faculty of Medicine and Health, University of Sydney, Camperdown, New South Wales, Australia
| | - Paul J White
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Peter Scammells
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Lauren T May
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Patrick M Sexton
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Radostin Danev
- Graduate School of Medicine, University of Tokyo, Tokyo, Japan
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, USA
| | - Alisa Glukhova
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia.
- Structural Biology Division, Walter and Eliza Hall Institute of Medical Research, Parkville, Victoria, Australia.
- Department of Biochemistry and Pharmacology, University of Melbourne, Parkville, Victoria, Australia.
| | - Wendy L Imlach
- Department of Physiology, Monash Biomedicine Discovery Institute, Monash University, Clayton, Victoria, Australia.
| | - Arthur Christopoulos
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia.
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia.
| |
Collapse
|
30
|
An X, Bai Q, Bing Z, Liu H, Yao X. Insights into the molecular mechanism of positive cooperativity between partial agonist MK-8666 and full allosteric agonist AP8 of hGPR40 by Gaussian accelerated molecular dynamics (GaMD) simulations. Comput Struct Biotechnol J 2021; 19:3978-3989. [PMID: 34377364 PMCID: PMC8313488 DOI: 10.1016/j.csbj.2021.07.008] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Revised: 07/07/2021] [Accepted: 07/08/2021] [Indexed: 10/29/2022] Open
Abstract
Activation of human free fatty acid receptor 1 (FFAR1, also called hGPR40) enhances insulin secretion in a glucose-dependent manner. Hence, the development of selective agonist targeting hGPR40 has been proposed as a therapeutic strategy of type 2 diabetes mellitus. Some agonists targeting hGPR40 were reported. The radioligand-binding studies and the crystal structures reveal that there are multiple sites on GPR40, and there exists positive binding cooperativity between the partial agonist MK-8666 and full allosteric agonist (AgoPAM) AP8. In this work, we carried out long-time Gaussian accelerated molecular dynamics (GaMD) simulations on hGPR40 to shed light on the mechanism of the cooperativity between the two agonists at different sites. Our results reveal that the induced-fit conformational coupling is bidirectional between the two sites. The movements and rotations of TM3, TM4, TM5 and TM6 due to their inherent flexibility are crucial in coupling the conformational changes of the two agonists binding sites. These helices adopt similar conformational states upon alternative ligand or both ligands binding. The Leu1384.57, Leu1865.42 and Leu1905.46 play roles in coordinating the rearrangements of residues in the two pockets, which makes the movements of residues in the two sites like gear movements. These results provide detailed information at the atomic level about the conformational coupling between different sites of GPR40, and also provide the structural information for further design of new agonists of GPR40. In addition, these results suggest that it is necessary by considering the effect of other site bound in structure-based ligands discovery.
Collapse
Affiliation(s)
- Xiaoli An
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry, Lanzhou University, Lanzhou 730000, China
| | - Qifeng Bai
- School of Basic Medical Science, Lanzhou University, Lanzhou, China
| | - Zhitong Bing
- Institute of Modern Physics of Chinese Academy of Sciences, Gansu Province, Lanzhou, China
| | - Huanxiang Liu
- School of Pharmacy, Lanzhou University, Lanzhou 730000, China
| | - Xiaojun Yao
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry, Lanzhou University, Lanzhou 730000, China.,State Key Laboratory of Quality Research in Chinese Medicine, Macau Institute for Applied Research in Medicine and Health, Macau University of Science and Technology, Taipa, Macau, China
| |
Collapse
|
31
|
Abstract
Extracellular nucleosides and nucleotides have widespread functions in responding to physiological stress. The "purinome" encompasses 4 G-protein-coupled receptors (GPCRs) for adenosine, 8 GPCRs activated by nucleotides, 7 adenosine 5'-triphosphate-gated P2X ion channels, as well as the associated enzymes and transporters that regulate native agonist levels. Purinergic signaling modulators, such as receptor agonists and antagonists, have potential for treating chronic pain. Adenosine and its analogues potently suppress nociception in preclinical models by activating A1 and/or A3 adenosine receptors (ARs), but safely harnessing this pathway to clinically treat pain has not been achieved. Both A2AAR agonists and antagonists are efficacious in pain models. Highly selective A3AR agonists offer a novel approach to treat chronic pain. We have explored the structure activity relationship of nucleoside derivatives at this subtype using a computational structure-based approach. Novel A3AR agonists for pain control containing a bicyclic ring system (bicyclo [3.1.0] hexane) in place of ribose were designed and screened using an in vivo phenotypic model, which reflected both pharmacokinetic and pharmacodynamic parameters. High specificity (>10,000-fold selective for A3AR) was achieved with the aid of receptor homology models based on related GPCR structures. These A3AR agonists are well tolerated in vivo and highly efficacious in models of chronic neuropathic pain. Furthermore, signaling molecules acting at P2X3, P2X4, P2X7, and P2Y12Rs play critical roles in maladaptive pain neuroplasticity, and their antagonists reduce chronic or inflammatory pain, and, therefore, purine receptor modulation is a promising approach for future pain therapeutics. Structurally novel antagonists for these nucleotide receptors were discovered recently.
Collapse
|
32
|
Do HN, Akhter S, Miao Y. Pathways and Mechanism of Caffeine Binding to Human Adenosine A 2A Receptor. Front Mol Biosci 2021; 8:673170. [PMID: 33987207 PMCID: PMC8111288 DOI: 10.3389/fmolb.2021.673170] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Accepted: 03/24/2021] [Indexed: 11/13/2022] Open
Abstract
Caffeine (CFF) is a common antagonist to the four subtypes of adenosine G-protein-coupled receptors (GPCRs), which are critical drug targets for treating heart failure, cancer, and neurological diseases. However, the pathways and mechanism of CFF binding to the target receptors remain unclear. In this study, we have performed all-atom-enhanced sampling simulations using a robust Gaussian-accelerated molecular dynamics (GaMD) method to elucidate the binding mechanism of CFF to human adenosine A2A receptor (A2AAR). Multiple 500–1,000 ns GaMD simulations captured both binding and dissociation of CFF in the A2AAR. The GaMD-predicted binding poses of CFF were highly consistent with the x-ray crystal conformations with a characteristic hydrogen bond formed between CFF and residue N6.55 in the receptor. In addition, a low-energy intermediate binding conformation was revealed for CFF at the receptor extracellular mouth between ECL2 and TM1. While the ligand-binding pathways of the A2AAR were found similar to those of other class A GPCRs identified from previous studies, the ECL2 with high sequence divergence serves as an attractive target site for designing allosteric modulators as selective drugs of the A2AAR.
Collapse
Affiliation(s)
- Hung N Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| | - Sana Akhter
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| |
Collapse
|
33
|
Deganutti G, Barkan K, Ladds G, Reynolds CA. Multisite Model of Allosterism for the Adenosine A1 Receptor. J Chem Inf Model 2021; 61:2001-2015. [PMID: 33779168 DOI: 10.1021/acs.jcim.0c01331] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Despite being a target for about one-third of approved drugs, G protein-coupled receptors (GPCRs) still represent a tremendous reservoir for therapeutic strategies against diseases. For example, several cardiovascular and central nervous system conditions could benefit from clinical agents that activate the adenosine 1 receptor (A1R); however, the pursuit of A1R agonists for clinical use is usually impeded by both on- and off-target side effects. One of the possible strategies to overcome this issue is the development of positive allosteric modulators (PAMs) capable of selectively enhancing the effect of a specific receptor subtype and triggering functional selectivity (a phenomenon also referred to as bias). Intriguingly, besides enforcing the effect of agonists upon binding to an allosteric site, most of the A1R PAMs display intrinsic partial agonism and orthosteric competition with antagonists. To rationalize this behavior, we simulated the binding of the prototypical PAMs PD81723 and VCP171, the full-agonist NECA, the antagonist 13B, and the bitopic agonist VCP746. We propose that a single PAM can bind several A1R sites rather than a unique allosteric pocket, reconciling the structure-activity relationship and the mutagenesis results.
Collapse
Affiliation(s)
- Giuseppe Deganutti
- Centre for Sport, Exercise and Life Sciences, Faculty of Health and Life Sciences, Coventry University, Alison Gingell Building, Coventry CV1 5FB, U.K
| | - Kerry Barkan
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Graham Ladds
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Christopher A Reynolds
- Centre for Sport, Exercise and Life Sciences, Faculty of Health and Life Sciences, Coventry University, Alison Gingell Building, Coventry CV1 5FB, U.K
| |
Collapse
|
34
|
Deganutti G, Barkan K, Preti B, Leuenberger M, Wall M, Frenguelli BG, Lochner M, Ladds G, Reynolds CA. Deciphering the Agonist Binding Mechanism to the Adenosine A1 Receptor. ACS Pharmacol Transl Sci 2021; 4:314-326. [PMID: 33615181 PMCID: PMC7887845 DOI: 10.1021/acsptsci.0c00195] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2020] [Indexed: 12/21/2022]
Abstract
Despite being among the most characterized G protein-coupled receptors (GPCRs), adenosine receptors (ARs) have always been a difficult target in drug design. To date, no agonist other than the natural effector and the diagnostic regadenoson has been approved for human use. Recently, the structure of the adenosine A1 receptor (A1R) was determined in the active, Gi protein complexed state; this has important repercussions for structure-based drug design. Here, we employed supervised molecular dynamics simulations and mutagenesis experiments to extend the structural knowledge of the binding of selective agonists to A1R. Our results identify new residues involved in the association and dissociation pathway, they suggest the binding mode of N6-cyclopentyladenosine (CPA) related ligands, and they highlight the dramatic effect that chemical modifications can have on the overall binding mechanism, paving the way for the rational development of a structure-kinetics relationship of A1R agonists.
Collapse
Affiliation(s)
- Giuseppe Deganutti
- Centre
for Sport, Exercise and Life Sciences, Faculty of Health and Life
Sciences, Coventry University, Alison Gingell Building, Coventry CV1 5FB, U.K.
| | - Kerry Barkan
- Department
of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge, CB2 1PD, U.K.
| | - Barbara Preti
- Institute
of Biochemistry and Molecular Medicine, University of Bern, 3012 Bern, Switzerland
| | - Michele Leuenberger
- Institute
of Biochemistry and Molecular Medicine, University of Bern, 3012 Bern, Switzerland
| | - Mark Wall
- School
of Life Sciences, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Bruno G. Frenguelli
- School
of Life Sciences, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Martin Lochner
- Institute
of Biochemistry and Molecular Medicine, University of Bern, 3012 Bern, Switzerland
| | - Graham Ladds
- Department
of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge, CB2 1PD, U.K.
| | - Christopher A. Reynolds
- Centre
for Sport, Exercise and Life Sciences, Faculty of Health and Life
Sciences, Coventry University, Alison Gingell Building, Coventry CV1 5FB, U.K.
| |
Collapse
|
35
|
Bhattarai A, Pawnikar S, Miao Y. Mechanism of ligand recognition by human ACE2 receptor. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2020. [PMID: 33140043 PMCID: PMC7605550 DOI: 10.1101/2020.10.30.362749] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Angiotensin converting enzyme 2 (ACE2) plays a key role in renin-angiotensin system regulation and amino acid homeostasis. Human ACE2 acts as the receptor for severe acute respiratory syndrome coronaviruses SARS-CoV and SARS-CoV-2. ACE2 is also widely expressed in epithelial cells of lungs, heart, kidney and pancreas. It is considered an important drug target for treating SARS-CoV-2, as well as pulmonary diseases, heart failure, hypertension, renal diseases and diabetes. Despite the critical importance, the mechanism of ligand binding to the human ACE2 receptor remains unknown. Here, we address this challenge through all-atom simulations using a novel ligand Gaussian accelerated molecular dynamics (LiGaMD) method. Microsecond LiGaMD simulations have successfully captured both binding and unbinding of the MLN-4760 inhibitor in the ACE2 receptor. In the ligand unbound state, the ACE2 receptor samples distinct Open, Partially Open and Closed conformations. Ligand binding biases the receptor conformational ensemble towards the Closed state. The LiGaMD simulations thus suggest a conformational selection mechanism for ligand recognition by the ACE2 receptor. Our simulation findings are expected to facilitate rational drug design of ACE2 against coronaviruses and other related human diseases.
Collapse
Affiliation(s)
- Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| | - Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| |
Collapse
|
36
|
Miao Y, Bhattarai A, Wang J. Ligand Gaussian Accelerated Molecular Dynamics (LiGaMD): Characterization of Ligand Binding Thermodynamics and Kinetics. J Chem Theory Comput 2020; 16:5526-5547. [PMID: 32692556 DOI: 10.1021/acs.jctc.0c00395] [Citation(s) in RCA: 105] [Impact Index Per Article: 26.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Calculations of ligand binding free energies and kinetic rates are important for drug design. However, such tasks have proven challenging in computational chemistry and biophysics. To address this challenge, we have developed a new computational method, ligand Gaussian accelerated molecular dynamics (LiGaMD), which selectively boosts the ligand nonbonded interaction potential energy based on the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling technique. Another boost potential could be applied to the remaining potential energy of the entire system in a dual-boost algorithm (LiGaMD_Dual) to facilitate ligand binding. LiGaMD has been demonstrated on host-guest and protein-ligand binding model systems. Repetitive guest binding and unbinding in the β-cyclodextrin host were observed in hundreds-of-nanosecond LiGaMD_Dual simulations. The calculated guest binding free energies agreed excellently with experimental data with <1.0 kcal/mol errors. Compared with converged microsecond-time scale conventional molecular dynamics simulations, the sampling errors of LiGaMD_Dual simulations were also <1.0 kcal/mol. Accelerations of ligand kinetic rate constants in LiGaMD simulations were properly estimated using Kramers' rate theory. Furthermore, LiGaMD allowed us to capture repetitive dissociation and binding of the benzamidine inhibitor in trypsin within 1 μs simulations. The calculated ligand binding free energy and kinetic rate constants compared well with the experimental data. In summary, LiGaMD provides a powerful enhanced sampling approach for characterizing ligand binding thermodynamics and kinetics simultaneously, which is expected to facilitate computer-aided drug design.
Collapse
Affiliation(s)
- Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
37
|
Allosteric interactions in the parathyroid hormone GPCR-arrestin complex formation. Nat Chem Biol 2020; 16:1096-1104. [PMID: 32632293 PMCID: PMC7502484 DOI: 10.1038/s41589-020-0567-0] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Accepted: 05/13/2020] [Indexed: 11/09/2022]
Abstract
Peptide ligands of class B G-protein-coupled receptors act via a two-step binding process, but the essential mechanisms that link their extracellular binding to intracellular receptor-arrestin interactions are not fully understood. Using NMR, crosslinking coupled to mass spectrometry, signaling experiments and computational approaches on the parathyroid hormone (PTH) type 1 receptor (PTHR), we show that initial binding of the PTH C-terminal part constrains the conformation of the flexible PTH N-terminal signaling epitope before a second binding event occurs. A 'hot-spot' PTH residue, His9, that inserts into the PTHR transmembrane domain at this second step allosterically engages receptor-arrestin coupling. A conformational change in PTHR intracellular loop 3 permits favorable interactions with β-arrestin's finger loop. These results unveil structural determinants for PTHR-arrestin complex formation and reveal that the two-step binding mechanism proceeds via cooperative fluctuations between ligand and receptor, which extend to other class B G-protein-coupled receptors.
Collapse
|
38
|
Bitencourt-Ferreira G, de Azevedo WF. Molecular Dynamics Simulations with NAMD2. Methods Mol Biol 2020; 2053:109-124. [PMID: 31452102 DOI: 10.1007/978-1-4939-9752-7_8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
X-ray diffraction crystallography is the primary technique to determine the three-dimensional structures of biomolecules. Although a robust method, X-ray crystallography is not able to access the dynamical behavior of macromolecules. To do so, we have to carry out molecular dynamics simulations taking as an initial system the three-dimensional structure obtained from experimental techniques or generated using homology modeling. In this chapter, we describe in detail a tutorial to carry out molecular dynamics simulations using the program NAMD2. We chose as a molecular system to simulate the structure of human cyclin-dependent kinase 2.
Collapse
Affiliation(s)
- Gabriela Bitencourt-Ferreira
- Escola de Ciências da Saúde, Pontifícia Universidade Católica do Rio Grande do Sul-PUCRS, Porto Alegre, RS, Brazil
| | - Walter Filgueira de Azevedo
- Escola de Ciências da Saúde, Pontifícia Universidade Católica do Rio Grande do Sul-PUCRS, Porto Alegre, RS, Brazil.
| |
Collapse
|
39
|
Pathway and mechanism of drug binding to chemokine receptors revealed by accelerated molecular simulations. Future Med Chem 2020; 12:1213-1225. [PMID: 32515227 DOI: 10.4155/fmc-2020-0044] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Background: Chemokine GPCRs play key roles in biology and medicine. Particularly, CXCR4 promotes cancer metastasis and facilitate HIV entry into host cells. Plerixafor (PLX) is a CXCR4 drug, but the pathway and binding site of PLX in CXCR4 remain unknown. Results & methodology: We have performed molecular docking and all-atom simulations using Gaussian accelerated molecular dynamics (GaMD), which are consistent with previous mutation experiments, suggesting that PLX binds to the orthosteric site of CXCR4 as an antagonist. The GaMD simulations further revealed an intermediate allosteric binding site at the extracellular mouth of CXCR4. Conclusion: The newly identified allosteric site can be targeted for novel drug design targeting CXCR4 and other chemokine receptors.
Collapse
|
40
|
Salmaso V, Jacobson KA. In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors. Biomolecules 2020; 10:E812. [PMID: 32466404 PMCID: PMC7356333 DOI: 10.3390/biom10060812] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 05/21/2020] [Accepted: 05/22/2020] [Indexed: 12/14/2022] Open
Abstract
Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.
Collapse
Affiliation(s)
| | - Kenneth A. Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA;
| |
Collapse
|
41
|
Bhattarai A, Wang J, Miao Y. Retrospective ensemble docking of allosteric modulators in an adenosine G-protein-coupled receptor. Biochim Biophys Acta Gen Subj 2020; 1864:129615. [PMID: 32298791 DOI: 10.1016/j.bbagen.2020.129615] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2020] [Revised: 02/26/2020] [Accepted: 04/08/2020] [Indexed: 12/18/2022]
Abstract
BACKGROUND Ensemble docking has proven useful in drug discovery and development. It increases the hit rate by incorporating receptor flexibility into molecular docking as demonstrated on important drug targets including G-protein-coupled receptors (GPCRs). Adenosine A1 receptor (A1AR) is a key GPCR that has been targeted for treating cardiac ischemia-reperfusion injuries, neuropathic pain and renal diseases. Development of allosteric modulators, compounds binding to distinct and less conserved GPCR target sites compared with agonists and antagonists, has attracted increasing interest for designing selective drugs of the A1AR. Despite significant advances, more effective approaches are needed to discover potent and selective allosteric modulators of the A1AR. METHODS Ensemble docking that integrates Gaussian accelerated molecular dynamic (GaMD) simulations and molecular docking using Autodock has been implemented for retrospective docking of known positive allosteric modulators (PAMs) in the A1AR. RESULTS Ensemble docking outperforms docking of the receptor cryo-EM structure. The calculated docking enrichment factors (EFs) and the area under the receiver operating characteristic curves (AUC) are significantly increased. CONCLUSIONS Receptor ensembles generated from GaMD simulations are able to increase the success rate of discovering PAMs of A1AR. It is important to account for receptor flexibility through GaMD simulations and flexible docking. GENERAL SIGNIFICANCE Ensemble docking is a promising approach for drug discovery targeting flexible receptors.
Collapse
Affiliation(s)
- Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA.
| |
Collapse
|
42
|
Deganutti G, Moro S, Reynolds CA. A Supervised Molecular Dynamics Approach to Unbiased Ligand–Protein Unbinding. J Chem Inf Model 2020; 60:1804-1817. [DOI: 10.1021/acs.jcim.9b01094] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Giuseppe Deganutti
- School of Life Sciences, University of Essex, Wivenhoe Park, Colchester, CO4 3SQ, United Kingdom
| | - Stefano Moro
- Molecular Modeling Section, Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Via Marzolo 5, 35131, Padova, Italy
| | - Christopher A. Reynolds
- School of Life Sciences, University of Essex, Wivenhoe Park, Colchester, CO4 3SQ, United Kingdom
| |
Collapse
|
43
|
Bhattarai A, Wang J, Miao Y. G-Protein-Coupled Receptor-Membrane Interactions Depend on the Receptor Activation State. J Comput Chem 2020; 41:460-471. [PMID: 31602675 PMCID: PMC7026935 DOI: 10.1002/jcc.26082] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2019] [Revised: 07/25/2019] [Accepted: 09/03/2019] [Indexed: 12/20/2022]
Abstract
G-protein-coupled receptors (GPCRs) are the largest family of human membrane proteins and serve as primary targets of approximately one-third of currently marketed drugs. In particular, adenosine A1 receptor (A1 AR) is an important therapeutic target for treating cardiac ischemia-reperfusion injuries, neuropathic pain, and renal diseases. As a prototypical GPCR, the A1 AR is located within a phospholipid membrane bilayer and transmits cellular signals by changing between different conformational states. It is important to elucidate the lipid-protein interactions in order to understand the functional mechanism of GPCRs. Here, all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method were performed on both the inactive (antagonist bound) and active (agonist and G-protein bound) A1 AR, which was embedded in a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) lipid bilayer. In the GaMD simulations, the membrane lipids played a key role in stabilizing different conformational states of the A1 AR. Our simulations further identified important regions of the receptor that interacted distinctly with the lipids in highly correlated manner. Activation of the A1 AR led to differential dynamics in the upper and lower leaflets of the lipid bilayer. In summary, GaMD enhanced simulations have revealed strongly coupled dynamics of the GPCR and lipids that depend on the receptor activation state. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| |
Collapse
|
44
|
Jones LH, Kelly JW. Structure-based design and analysis of SuFEx chemical probes. RSC Med Chem 2020; 11:10-17. [PMID: 33479601 DOI: 10.1039/c9md00542k] [Citation(s) in RCA: 74] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2019] [Accepted: 12/16/2019] [Indexed: 12/11/2022] Open
Abstract
The discerning reactivity of sulfur(vi)-fluoride exchange (SuFEx) chemistry has enabled the context-specific labeling of protein binding sites by chemical probes that incorporate these versatile warheads. Emerging information from protein-probe structures and proteomic mapping experiments is helping advance our understanding of the protein microenvironment that dictates the reactivity of targetable amino acid residues. This review explores these new findings that should influence the future rational design of SuFEx probes for a multitude of applications in chemical biology and drug discovery.
Collapse
Affiliation(s)
- Lyn H Jones
- Center for Protein Degradation , Dana-Farber Cancer Institute , 360 Longwood Avenue , Boston , MA 02215 , USA .
| | - Jeffery W Kelly
- Departments of Chemistry and Molecular Medicine , The Scripps Research Institute , La Jolla , CA , USA
| |
Collapse
|
45
|
Balogh G, Komáromi I, Bereczky Z. The mechanism of high affinity pentasaccharide binding to antithrombin, insights from Gaussian accelerated molecular dynamics simulations. J Biomol Struct Dyn 2019; 38:4718-4732. [PMID: 31686597 DOI: 10.1080/07391102.2019.1688194] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from 'static' structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a 'non-activated' AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. 'Generalized correlation' matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Gábor Balogh
- Division of Clinical Laboratory Science, Department of Laboratory Medicine, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.,Kálmán Laki Doctoral School, University of Debrecen, Debrecen, Hungary
| | - István Komáromi
- Division of Clinical Laboratory Science, Department of Laboratory Medicine, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Zsuzsanna Bereczky
- Division of Clinical Laboratory Science, Department of Laboratory Medicine, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| |
Collapse
|
46
|
Li Y, Wang M, Gao N, Li D, Lin J. The effect of dimerization on the activation and conformational dynamics of adenosine A 1 receptor. Phys Chem Chem Phys 2019; 21:22763-22773. [PMID: 31595279 DOI: 10.1039/c9cp04060a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The adenosine A1 receptor (A1R) is one of four adenosine receptors in humans, which are involved in the function of the cardiovascular, respiratory and central nervous systems. Experimental results indicate that A1R can form a homodimer and that the protomer-protomer interaction in the A1R dimer is related to certain pharmacological characteristics of A1R activation. In this work, we performed docking, metadynamics simulation, conventional molecular dynamics simulations, Gaussian-accelerated molecular dynamics simulations, potential of mean force calculations, dynamic cross-correlation motions analysis and community network analysis to study the binding mode of 5'-N-ethylcarboxamidoadenosine (NECA) to A1R and the effect of dimerization on the activation of A1R. Our results show that NECA binds to A1R in a similar mode to adenosine in the A1R crystal structure and NECA in the A2AR crystal structure. The A1R homodimer can be activated by one or two agonists with NECA occupying its orthosteric pockets in one (which we call the NECA-A1R system) or both protomers (which we call the dNECA-A1R system). In the NECA-A1R system, activation is predicated in the protomer without NECA bound. In the dNECA-A1R system, only one protomer achieves the active state. These findings suggest an asymmetrical activation mechanism of the homodimer and a negative cooperativity between the two protomers. We envision that our results may further facilitate the drug development of A1R.
Collapse
Affiliation(s)
- Yang Li
- State Key Laboratory of Medicinal Chemical Biology, College of Pharmacy and Tianjin Key Laboratory of Molecular Drug Research, Nankai University, Haihe Education Park, 38 Tongyan Road, Tianjin 300353, People's Republic of China.
| | | | | | | | | |
Collapse
|
47
|
Wang J, Miao Y. Mechanistic Insights into Specific G Protein Interactions with Adenosine Receptors. J Phys Chem B 2019; 123:6462-6473. [PMID: 31283874 DOI: 10.1021/acs.jpcb.9b04867] [Citation(s) in RCA: 75] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Coupling between G-protein-coupled receptors (GPCRs) and the G proteins is a key step in cellular signaling. Despite extensive experimental and computational studies, the mechanism of specific GPCR-G protein coupling remains poorly understood. This has greatly hindered effective drug design of GPCRs that are primary targets of ∼1/3 of currently marketed drugs. Here, we have employed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method to decipher the mechanism of the GPCR-G protein interactions. Adenosine receptors (ARs) were used as model systems based on very recently determined cryo-EM structures of the A1AR and A2AAR coupled with the Gi and Gs proteins, respectively. Changing the Gi protein to the Gs led to increased fluctuations in the A1AR and agonist adenosine (ADO), while agonist 5'-N-ethylcarboxamidoadenosine (NECA) binding in the A2AAR could be still stabilized upon changing the Gs protein to the Gi. Free energy calculations identified one stable low-energy conformation for each of the A1AR-Gi and A2AAR-Gs complexes as in the cryo-EM structures, similarly for the A2AAR-Gi complex. In contrast, the ADO agonist and Gs protein sampled multiple conformations in the A1AR-Gs system. GaMD simulations thus indicated that the A1AR preferred to couple with the Gi protein to the Gs, while the A2AAR could couple with both the Gs and Gi proteins, being highly consistent with experimental findings of the ARs. More importantly, detailed analysis of the atomic simulations showed that the specific AR-G protein coupling resulted from remarkably complementary residue interactions at the protein interface, involving mainly the receptor transmembrane 6 helix and the Gα α5 helix and α4-β6 loop. In summary, the GaMD simulations have provided unprecedented insights into the dynamic mechanism of specific GPCR-G protein interactions at an atomistic level.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences , University of Kansas , Lawrence , Kansas 66047 , United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences , University of Kansas , Lawrence , Kansas 66047 , United States
| |
Collapse
|
48
|
Mahmod Al-Qattan MN, Mordi MN. Molecular Basis of Modulating Adenosine Receptors Activities. Curr Pharm Des 2019; 25:817-831. [DOI: 10.2174/1381612825666190304122624] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2019] [Accepted: 02/26/2019] [Indexed: 01/04/2023]
Abstract
Modulating cellular processes through extracellular chemical stimuli is medicinally an attractive approach to control disease conditions. GPCRs are the most important group of transmembranal receptors that produce different patterns of activations using intracellular mediators (such as G-proteins and Beta-arrestins). Adenosine receptors (ARs) belong to GPCR class and are divided into A1AR, A2AAR, A2BAR and A3AR. ARs control different physiological activities thus considered valuable target to control neural, heart, inflammatory and other metabolic disorders. Targeting ARs using small molecules essentially works by binding orthosteric and/or allosteric sites of the receptors. Although targeting orthosteric site is considered typical to modulate receptor activity, allosteric sites provide better subtype selectivity, saturable modulation of activity and variable activation patterns. Each receptor exists in dynamical equilibrium between conformational ensembles. The equilibrium is affected by receptor interaction with other molecules. Changing the population of conformational ensembles of the receptor is the method by which orthosteric, allosteric and other cellular components control receptor signaling. Herein, the interactions of ARs with orthosteric, allosteric ligands as well as intracellular mediators are described. A quinary interaction model for the receptor is proposed and energy wells for major conformational ensembles are retrieved.
Collapse
Affiliation(s)
| | - Mohd Nizam Mordi
- Centre For Drug Research, Universiti Sains Malaysia, 11800 Gelugor, Penang, Malaysia
| |
Collapse
|
49
|
Lamim Ribeiro JM, Filizola M. Allostery in G protein-coupled receptors investigated by molecular dynamics simulations. Curr Opin Struct Biol 2019; 55:121-128. [PMID: 31096158 DOI: 10.1016/j.sbi.2019.03.016] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Revised: 03/07/2019] [Accepted: 03/11/2019] [Indexed: 01/14/2023]
Abstract
G-protein-coupled receptors (GPCRs) are allosteric signaling machines that trigger distinct functional responses depending on the particular conformational state they adopt upon binding. This so-called GPCR functional selectivity is prompted by ligands of different efficacy binding at orthosteric or allosteric sites on the receptor, as well as by interactions with intracellular protein partners or other receptor types. Molecular dynamics (MD) simulations can provide important mechanistic, thermodynamic, and kinetic insights into these interactions at a level of molecular detail that is necessary to rightly inform modern drug discovery. Here, we review the most recent MD contributions to understanding GPCR allostery, with an emphasis on their strengths and limitations.
Collapse
Affiliation(s)
- João Marcelo Lamim Ribeiro
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, One Gustave L. Levy Place, Box 1677, New York, NY, 10029, USA
| | - Marta Filizola
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, One Gustave L. Levy Place, Box 1677, New York, NY, 10029, USA.
| |
Collapse
|
50
|
Jacobson KA, Tosh DK, Jain S, Gao ZG. Historical and Current Adenosine Receptor Agonists in Preclinical and Clinical Development. Front Cell Neurosci 2019; 13:124. [PMID: 30983976 PMCID: PMC6447611 DOI: 10.3389/fncel.2019.00124] [Citation(s) in RCA: 136] [Impact Index Per Article: 27.2] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2019] [Accepted: 03/13/2019] [Indexed: 12/22/2022] Open
Abstract
Adenosine receptors (ARs) function in the body’s response to conditions of pathology and stress associated with a functional imbalance, such as in the supply and demand of energy/oxygen/nutrients. Extracellular adenosine concentrations vary widely to raise or lower the basal activation of four subtypes of ARs. Endogenous adenosine can correct an energy imbalance during hypoxia and other stress, for example, by slowing the heart rate by A1AR activation or increasing the blood supply to heart muscle by the A2AAR. Moreover, exogenous AR agonists, antagonists, or allosteric modulators can be applied for therapeutic benefit, and medicinal chemists working toward that goal have reported thousands of such agents. Thus, numerous clinical trials have ensued, using promising agents to modulate adenosinergic signaling, most of which have not succeeded. Currently, short-acting, parenteral agonists, adenosine and Regadenoson, are the only AR agonists approved for human use. However, new concepts and compounds are currently being developed and applied toward preclinical and clinical evaluation, and initial results are encouraging. This review focuses on key compounds as AR agonists and positive allosteric modulators (PAMs) for disease treatment or diagnosis. AR agonists for treating inflammation, pain, cancer, non-alcoholic steatohepatitis, angina, sickle cell disease, ischemic conditions and diabetes have been under development. Multiple clinical trials with two A3AR agonists are ongoing.
Collapse
Affiliation(s)
- Kenneth A Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, United States
| | - Dilip K Tosh
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, United States
| | - Shanu Jain
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, United States
| | - Zhan-Guo Gao
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, United States
| |
Collapse
|