1
|
Stampelou M, Ladds G, Kolocouris A. Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A 3 Receptor. J Phys Chem B 2024; 128:914-936. [PMID: 38236582 DOI: 10.1021/acs.jpcb.3c05986] [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/19/2024]
Abstract
A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.
Collapse
Affiliation(s)
- Margarita Stampelou
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Graham Ladds
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Antonios Kolocouris
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| |
Collapse
|
2
|
Guidetti L, Zappia A, Scalvini L, Ferrari FR, Giorgio C, Castelli R, Galvani F, Vacondio F, Rivara S, Mor M, Urbinati C, Rusnati M, Tognolini M, Lodola A. Molecular Determinants of EphA2 and EphB2 Antagonism Enable the Design of Ligands with Improved Selectivity. J Chem Inf Model 2023; 63:6900-6911. [PMID: 37910792 PMCID: PMC10647059 DOI: 10.1021/acs.jcim.3c01064] [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] [Received: 07/13/2023] [Revised: 10/13/2023] [Accepted: 10/18/2023] [Indexed: 11/03/2023]
Abstract
With the aim of identifying novel antagonists selective for the EphA receptor family, a combined experimental and computational approach was taken to investigate the molecular basis of the recognition between a prototypical Eph-ephrin antagonist (UniPR1447) and two representative receptors of the EphA and EphB subfamilies, namely, EphA2 and EphB2 receptors. The conformational free-energy surface (FES) of the binding state of UniPR1447 within the ligand binding domain of EphA2 and EphB2, reconstructed from molecular dynamics (MD) simulations performed on the microsecond time scale, was exploited to drive the design and synthesis of a novel antagonist selective for EphA2 over the EphB2 receptor. The availability of compounds with this pharmacological profile will help discriminate the importance of these two receptors in the insurgence and progression of cancer.
Collapse
Affiliation(s)
- Lorenzo Guidetti
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Alfonso Zappia
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Laura Scalvini
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Francesca Romana Ferrari
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Carmine Giorgio
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Riccardo Castelli
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Francesca Galvani
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Federica Vacondio
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Silvia Rivara
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Marco Mor
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
- Microbiome
Research Hub, Università degli Studi
di Parma, Parco Area
delle scienze 11/A, I- 43124 Parma, Italy
| | - Chiara Urbinati
- Dipartimento
di Medicina Molecolare Traslazionale, Università
degli Studi di Brescia, Brescia 25121, Italy
| | - Marco Rusnati
- Dipartimento
di Medicina Molecolare Traslazionale, Università
degli Studi di Brescia, Brescia 25121, Italy
| | - Massimiliano Tognolini
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| | - Alessio Lodola
- Dipartimento
di Scienze degli Alimenti e del Farmaco, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I- 43124 Parma, Italy
| |
Collapse
|
3
|
Buigues P, Gehrke S, Badaoui M, Dudas B, Mandana G, Qi T, Bottegoni G, Rosta E. Investigating the Unbinding of Muscarinic Antagonists from the Muscarinic 3 Receptor. J Chem Theory Comput 2023; 19:5260-5272. [PMID: 37458730 PMCID: PMC10413856 DOI: 10.1021/acs.jctc.3c00023] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Indexed: 08/09/2023]
Abstract
Patient symptom relief is often heavily influenced by the residence time of the inhibitor-target complex. For the human muscarinic receptor 3 (hMR3), tiotropium is a long-acting bronchodilator used in conditions such as asthma or chronic obstructive pulmonary disease (COPD). The mechanistic insights into this inhibitor remain unclear; specifically, the elucidation of the main factors determining the unbinding rates could help develop the next generation of antimuscarinic agents. Using our novel unbinding algorithm, we were able to investigate ligand dissociation from hMR3. The unbinding paths of tiotropium and two of its analogues, N-methylscopolamin and homatropine methylbromide, show a consistent qualitative mechanism and allow us to identify the structural bottleneck of the process. Furthermore, our machine learning-based analysis identified key roles of the ECL2/TM5 junction involved in the transition state. Additionally, our results point to relevant changes at the intracellular end of the TM6 helix leading to the ICL3 kinase domain, highlighting the closest residue L482. This residue is located right between two main protein binding sites involved in signal transduction for hMR3's activation and regulation. We also highlight key pharmacophores of tiotropium that play determining roles in the unbinding kinetics and could aid toward drug design and lead optimization.
Collapse
Affiliation(s)
- Pedro
J. Buigues
- Department
of Physics and Astronomy, University College
London, London WC1E 6BT, United
Kingdom
| | - Sascha Gehrke
- Department
of Physics and Astronomy, University College
London, London WC1E 6BT, United
Kingdom
| | - Magd Badaoui
- Department
of Physics and Astronomy, University College
London, London WC1E 6BT, United
Kingdom
| | - Balint Dudas
- Department
of Physics and Astronomy, University College
London, London WC1E 6BT, United
Kingdom
| | - Gaurav Mandana
- Department
of Physics and Astronomy, University College
London, London WC1E 6BT, United
Kingdom
| | - Tianyun Qi
- Department
of Physics and Astronomy, University College
London, London WC1E 6BT, United
Kingdom
| | - Giovanni Bottegoni
- Dipartimento
di Scienze Biomolecolari (DISB), University
of Urbino, Urbino Piazza Rinascimento, 6, Urbino 61029, Italy
- Institute
of Clinical Sciences, University of Birmingham, Edgbaston, B15 2TT Birmingham, United Kingdom
| | - Edina Rosta
- Department
of Physics and Astronomy, University College
London, London WC1E 6BT, United
Kingdom
| |
Collapse
|