1
|
Saporiti S, Bianchi D, Ben Mariem O, Rossi M, Guerrini U, Eberini I, Centola F. In silico evaluation of the role of Fab glycosylation in cetuximab antibody dynamics. Front Immunol 2024; 15:1429600. [PMID: 39185413 PMCID: PMC11342397 DOI: 10.3389/fimmu.2024.1429600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2024] [Accepted: 07/23/2024] [Indexed: 08/27/2024] Open
Abstract
Introduction N-glycosylation is a post-translational modification that is highly important for the development of monoclonal antibodies (mAbs), as it regulates their biological activity, particularly in terms of immune effector functions. While typically added at the Fc level, approximately 15-25% of circulating antibodies exhibit glycosylation in the Fab domains as well. To the best of our knowledge, cetuximab (Erbitux®) is the only therapeutic antibody presenting Fab glycosylation approved world-wide targeting the epidermal growth factor receptor for the treatment of metastatic-colorectal and head and neck cancers. Additionally, it can trigger antibody-dependent cell cytotoxicity (ADCC), a response that typically is influenced by N-glycosylation at Fc level. However, the role of Fab glycosylation in cetuximab remains poorly understood. Hence, this study aims to investigate the structural role of Fab glycosylation on the conformational behavior of cetuximab. Methods The study was performed in silico via accelerated molecular dynamics simulations. The commercial cetuximab was compared to its form without Fab glycosylation and structural descriptors were evaluated to establish conformational differences. Results The results clearly show a correlation between the Fab glycosylation and structural descriptors that may modulate the conformational freedom of the antibody, potentially affecting Fc effector functions, and suggesting a negative role of Fab glycosylation on the interaction with FcγRIIIa. Conclusion Fab glycosylation of cetuximab is the most critical challenge for biosimilar development, but the differences highlighted in this work with respect to its aglycosylated form can improve the knowledge and represent also a great opportunity to develop novel strategies of biotherapeutics.
Collapse
Affiliation(s)
- Simona Saporiti
- Analytical Excellence and Program Management, Merck Serono S.p.A., Rome, Italy
| | - Davide Bianchi
- Dipartimento di Scienze Farmacologiche e Biomolecolari, Università degli Studi di Milano, Milan, Italy
| | - Omar Ben Mariem
- Dipartimento di Scienze Farmacologiche e Biomolecolari, Università degli Studi di Milano, Milan, Italy
| | - Mara Rossi
- Analytical Excellence and Program Management, Merck Serono S.p.A., Rome, Italy
| | - Uliano Guerrini
- Dipartimento di Scienze Farmacologiche e Biomolecolari, Università degli Studi di Milano, Milan, Italy
| | - Ivano Eberini
- Dipartimento di Scienze Farmacologiche e Biomolecolari & Data Science Research Center (DSRC), Università degli Studi di Milano, Milan, Italy
| | - Fabio Centola
- Analytical Excellence and Program Management, Merck Serono S.p.A., Rome, Italy
| |
Collapse
|
2
|
Tu G, Gong Y, Yao X, Liu Q, Xue W, Zhang R. Pathways and mechanism of MRTX1133 binding to KRAS G12D elucidated by molecular dynamics simulations and Markov state models. Int J Biol Macromol 2024; 274:133374. [PMID: 38925182 DOI: 10.1016/j.ijbiomac.2024.133374] [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] [Received: 04/12/2024] [Revised: 06/19/2024] [Accepted: 06/21/2024] [Indexed: 06/28/2024]
Abstract
KRAS G12D is the most common oncogenic mutation identified in several types of cancer. Therefore, design of inhibitors targeting KRAS G12D represents a promising strategy for anticancer therapy. MRTX1133 is a highly potent inhibitor (approximate experiment Kd ≈ 0.0002 nM) of KRAS G12D and is currently in Phase 1/2 study, however, pathways of the compound binding to KRAS G12D has remained unknown, and the mechanism underlying the complicated dynamic process are challenging to capture experimentally, which hinder the structure-based anti-cancer drug design. Here, using MRTX1133 as a probe, unbiased molecular dynamics (MD) was used to simulate the process of MRTX1133 spontaneously binding to KRAS G12D. In six of 42 independent MD simulation (a total of 99 μs), MRTX1133 was observed to successfully associate with KRAS G12D. The kinetically metastable states refer to the potential pathways of MRTX1133 binding to KRAS G12D were revealed by Markov state models (MSM) analysis. Additionally, 8 key residues that are essential for MRTX1133 recognition and tight binding at the preferred low energy states were identified by MM/GBSA analysis. In sum, this study provides a new perspective on understanding the pathways and mechanism of MRTX1133 binding to KRAS G12D.
Collapse
Affiliation(s)
- Gao Tu
- Department of Pharmacy, The Second Affiliated Hospital, Army Medical University, 183 Xinqiao Road, Chongqing 400037, China; 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, 999078, Macau
| | - Yaguo Gong
- 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, 999078, Macau
| | - Xiaojun Yao
- Centre for Artificial Intelligence Driven Drug Discovery, Faculty of Applied Sciences, Macao Polytechnic University, 999078, Macau.
| | - Qing Liu
- Suzhou Institute for Advance Research, University of Science and Technology of China, Suzhou, China
| | - Weiwei Xue
- School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China.
| | - Rong Zhang
- Department of Pharmacy, The Second Affiliated Hospital, Army Medical University, 183 Xinqiao Road, Chongqing 400037, China.
| |
Collapse
|
3
|
Wang J, Miao Y. Ligand Gaussian Accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides. J Chem Theory Comput 2024; 20:5829-5841. [PMID: 39002136 DOI: 10.1021/acs.jctc.4c00502] [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: 07/15/2024]
Abstract
Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional molecular dynamics (cMD), due to limited simulation time scales. Based on our previously developed ligand Gaussian accelerated molecular dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3″, in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding, and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as the model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 μs simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 were in agreement with the available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligands and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina-Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina-Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
4
|
Gaiser BI, Danielsen M, Xu X, Røpke Jørgensen K, Fronik P, Märcher-Rørsted E, Wróbel TM, Liu X, Mosolff Mathiesen J, Sejer Pedersen D. Bitopic Ligands Support the Presence of a Metastable Binding Site at the β 2 Adrenergic Receptor. J Med Chem 2024; 67:11053-11068. [PMID: 38952152 DOI: 10.1021/acs.jmedchem.4c00578] [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: 07/03/2024]
Abstract
Metastable binding sites (MBS) have been observed in a multitude of molecular dynamics simulations and can be considered low affinity allosteric binding sites (ABS) that function as stepping stones as the ligand moves toward the orthosteric binding site (OBS). Herein, we show that MBS can be utilized as ABS in ligand design, resulting in ligands with improved binding kinetics. Four homobivalent bitopic ligands (1-4) were designed by molecular docking of (S)-alprenolol ((S)-ALP) in the cocrystal structure of the β2 adrenergic receptor (β2AR) bound to the antagonist ALP. Ligand 4 displayed a potency and affinity similar to (S)-ALP, but with a >4-fold increase in residence time. The proposed binding mode was confirmed by X-ray crystallography of ligand 4 in complex with the β2AR. This ligand design principle can find applications beyond the β2AR and G protein-coupled receptors (GPCRs) as a general approach for improving the pharmacological profile of orthosteric ligands by targeting the OBS and an MBS simultaneously.
Collapse
Affiliation(s)
- Birgit Isabel Gaiser
- Department of Drug Design and Pharmacology, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| | - Mia Danielsen
- Department of Drug Design and Pharmacology, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| | - Xinyu Xu
- State Key Laboratory of Membrane Biology, Tsinghua-Peking Center for Life Sciences, Beijing Frontier Research Center for Biological Structure, School of Pharmaceutical Sciences, Tsinghua University, Beijing 100084 ,China
| | - Kira Røpke Jørgensen
- Department of Drug Design and Pharmacology, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| | - Philipp Fronik
- Department of Drug Design and Pharmacology, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| | - Emil Märcher-Rørsted
- Department of Drug Design and Pharmacology, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| | - Tomasz M Wróbel
- Department of Drug Design and Pharmacology, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
- Department of Synthesis and Chemical Technology of Pharmaceutical Substances, Medical University of Lublin, Chodźki 4a, 20093 Lublin, Poland
| | - Xiangyu Liu
- State Key Laboratory of Membrane Biology, Tsinghua-Peking Center for Life Sciences, Beijing Frontier Research Center for Biological Structure, School of Pharmaceutical Sciences, Tsinghua University, Beijing 100084 ,China
| | - Jesper Mosolff Mathiesen
- Department of Drug Design and Pharmacology, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| | - Daniel Sejer Pedersen
- Department of Drug Design and Pharmacology, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| |
Collapse
|
5
|
Ghosh D, Biswas A, Radhakrishna M. Advanced computational approaches to understand protein aggregation. BIOPHYSICS REVIEWS 2024; 5:021302. [PMID: 38681860 PMCID: PMC11045254 DOI: 10.1063/5.0180691] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 03/18/2024] [Indexed: 05/01/2024]
Abstract
Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.
Collapse
Affiliation(s)
- Deepshikha Ghosh
- Department of Biological Sciences and Engineering, Indian Institute of Technology (IIT) Gandhinagar, Palaj, Gujarat 382355, India
| | - Anushka Biswas
- Department of Chemical Engineering, Indian Institute of Technology (IIT) Gandhinagar, Palaj, Gujarat 382355, India
| | | |
Collapse
|
6
|
Wang J, Miao Y. Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.06.592668. [PMID: 38766067 PMCID: PMC11100592 DOI: 10.1101/2024.05.06.592668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional Molecular Dynamics (cMD), due to limited simulation timescales. Based on our previously developed Ligand Gaussian accelerated Molecular Dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3", in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 microsecond simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 agreed with available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligand and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina – Chapel Hill, Chapel Hill, North Carolina, USA 27599
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina – Chapel Hill, Chapel Hill, North Carolina, USA 27599
| |
Collapse
|
7
|
Rahman MH, Hegazy L. Mechanism of antagonist ligand binding to REV-ERBα. Sci Rep 2024; 14:8401. [PMID: 38600172 PMCID: PMC11006950 DOI: 10.1038/s41598-024-58945-4] [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: 01/11/2024] [Accepted: 04/04/2024] [Indexed: 04/12/2024] Open
Abstract
REV-ERBα, a therapeutically promising nuclear hormone receptor, plays a crucial role in regulating various physiological processes such as the circadian clock, inflammation, and metabolism. However, the availability of chemical probes to investigate the pharmacology of this receptor is limited, with SR8278 being the only identified synthetic antagonist. Moreover, no X-ray crystal structures are currently available that demonstrate the binding of REV-ERBα to antagonist ligands. This lack of structural information impedes the development of targeted therapeutics. To address this issue, we employed Gaussian accelerated molecular dynamics (GaMD) simulations to investigate the binding pathway of SR8278 to REV-ERBα. For comparison, we also used GaMD to observe the ligand binding process of STL1267, for which an X-ray structure is available. GaMD simulations successfully captured the binding of both ligands to the receptor's orthosteric site and predicted the ligand binding pathway and important amino acid residues involved in the antagonist SR8278 binding. This study highlights the effectiveness of GaMD in investigating protein-ligand interactions, particularly in the context of drug recognition for nuclear hormone receptors.
Collapse
Affiliation(s)
- Mohammad Homaidur Rahman
- Center for Clinical Pharmacology, Washington University School of Medicine, University of Health Sciences and Pharmacy, St. Louis, MO, USA
- Department of Pharmaceutical and Administrative Sciences, University of Health Sciences and Pharmacy, St. Louis, MO, USA
| | - Lamees Hegazy
- Center for Clinical Pharmacology, Washington University School of Medicine, University of Health Sciences and Pharmacy, St. Louis, MO, USA.
- Department of Pharmaceutical and Administrative Sciences, University of Health Sciences and Pharmacy, St. Louis, MO, USA.
| |
Collapse
|
8
|
Chen J, Gou Q, Chen X, Song Y, Zhang F, Pu X. Exploring biased activation characteristics by molecular dynamics simulation and machine learning for the μ-opioid receptor. Phys Chem Chem Phys 2024; 26:10698-10710. [PMID: 38512140 DOI: 10.1039/d3cp05050e] [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: 03/22/2024]
Abstract
Biased ligands selectively activating specific downstream signaling pathways (termed as biased activation) exhibit significant therapeutic potential. However, the conformational characteristics revealed are very limited for the biased activation, which is not conducive to biased drug development. Motivated by the issue, we combine extensive accelerated molecular dynamics simulations and an interpretable deep learning model to probe the biased activation features for two complex systems constructed by the inactive μOR and two different biased agonists (G-protein-biased agonist TRV130 and β-arrestin-biased agonist endomorphin2). The results indicate that TRV130 binds deeper into the receptor core compared to endomorphin2, located between W2936.48 and D1142.50, and forms hydrogen bonding with D1142.50, while endomorphin2 binds above W2936.48. The G protein-biased agonist induces greater outward movements of the TM6 intracellular end, forming a typical active conformation, while the β-arrestin-biased agonist leads to a smaller extent of outward movements of TM6. Compared with TRV130, endomorphin2 causes more pronounced inward movements of the TM7 intracellular end and more complex conformational changes of H8 and ICL1. In addition, important residues determining the two different biased activation states were further identified by using an interpretable deep learning classification model, including some common biased activation residues across Class A GPCRs like some key residues on the TM2 extracellular end, ECL2, TM5 intracellular end, TM6 intracellular end, and TM7 intracellular end, and some specific important residues of ICL3 for μOR. The observations will provide valuable information for understanding the biased activation mechanism for GPCRs.
Collapse
Affiliation(s)
- Jianfang Chen
- College of Chemistry, Sichuan University, Chengdu 610064, China.
| | - Qiaoling Gou
- College of Chemistry, Sichuan University, Chengdu 610064, China.
| | - Xin Chen
- College of Chemistry, Sichuan University, Chengdu 610064, China.
| | - Yuanpeng Song
- College of Chemistry, Sichuan University, Chengdu 610064, China.
| | - Fuhui Zhang
- Graduate School, Sichuan University, Chengdu 610064, China
| | - Xuemei Pu
- College of Chemistry, Sichuan University, Chengdu 610064, China.
| |
Collapse
|
9
|
Montejo-López W, Sampieri-Cabrera R, Nicolás-Vázquez MI, Aceves-Hernández JM, Razo-Hernández RS. Analysing the effect caused by increasing the molecular volume in M1-AChR receptor agonists and antagonists: a structural and computational study. RSC Adv 2024; 14:8615-8640. [PMID: 38495977 PMCID: PMC10938299 DOI: 10.1039/d3ra07380g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Accepted: 03/04/2024] [Indexed: 03/19/2024] Open
Abstract
M1 muscarinic acetylcholine receptor (M1-AChR), a member of the G protein-coupled receptors (GPCR) family, plays a crucial role in learning and memory, making it an important drug target for Alzheimer's disease (AD) and schizophrenia. M1-AChR activation and deactivation have shown modifying effects in AD and PD preclinical models, respectively. However, understanding the pharmacology associated with M1-AChR activation or deactivation is complex, because of the low selectivity among muscarinic subtypes, hampering their therapeutic applications. In this regard, we constructed two quantitative structure-activity relationship (QSAR) models, one for M1-AChR agonists (total and partial), and the other for the antagonists. The binding mode of 59 structurally different compounds, including agonists and antagonists with experimental binding affinity values (pKi), were analyzed employing computational molecular docking over different structures of M1-AChR. Furthermore, we considered the interaction energy (Einter), the number of rotatable bonds (NRB), and lipophilicity (ilogP) for the construction of the QSAR model for agonists (R2 = 89.64, QLMO2 = 78, and Qext2 = 79.1). For the QSAR model of antagonists (R2 = 88.44, QLMO2 = 82, and Qext2 = 78.1) we considered the Einter, the fraction of sp3 carbons fCsp3, and lipophilicity (MlogP). Our results suggest that the ligand volume is a determinant to establish its biological activity (agonist or antagonist), causing changes in binding energy, and determining the affinity for M1-AChR.
Collapse
Affiliation(s)
- Wilber Montejo-López
- Departamento de Ciencias Químicas, Facultad de Estudios Superiores Cuautitlán Campo 1, Universidad Nacional Autónoma de México Avenida 1o de Mayo s/n, Colonia Santa María las Torres Cuautitlán Izcalli Estado de Mexico 54740 Mexico
| | - Raúl Sampieri-Cabrera
- Departamento de Fisiología, Facultad de Medicina, Universidad Nacional Autónoma de México, Centro de Ciencias de Complejidad, Universidad Nacional Autónoma de México Mexico
| | - María Inés Nicolás-Vázquez
- Departamento de Ciencias Químicas, Facultad de Estudios Superiores Cuautitlán Campo 1, Universidad Nacional Autónoma de México Avenida 1o de Mayo s/n, Colonia Santa María las Torres Cuautitlán Izcalli Estado de Mexico 54740 Mexico
| | - Juan Manuel Aceves-Hernández
- Unidad de Investigación Multidisciplinaria L14 (Alimentos, Micotoxinas, y Micotoxicosis), Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México Cuautitlán Izcalli Estado de Mexico 54714 Mexico
| | - Rodrigo Said Razo-Hernández
- Centro de Investigación en Dinámica Celular, Instituto de Investigación en Ciencias Básicas y Aplicadas, Universidad Autónoma del Estado de Morelos Av. Universidad 1001 Cuernavaca 62209 Mexico
| |
Collapse
|
10
|
Shenol A, Lückmann M, Trauelsen M, Lambrughi M, Tiberti M, Papaleo E, Frimurer TM, Schwartz TW. Molecular dynamics-based identification of binding pathways and two distinct high-affinity sites for succinate in succinate receptor 1/GPR91. Mol Cell 2024; 84:955-966.e4. [PMID: 38325379 DOI: 10.1016/j.molcel.2024.01.011] [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] [Received: 04/06/2023] [Revised: 11/30/2023] [Accepted: 01/16/2024] [Indexed: 02/09/2024]
Abstract
SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations.
Collapse
Affiliation(s)
- Aslihan Shenol
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Michael Lückmann
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Mette Trauelsen
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Matteo Lambrughi
- Cancer Structural Biology, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Matteo Tiberti
- Cancer Structural Biology, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Elena Papaleo
- Cancer Structural Biology, Danish Cancer Society Research Center, Copenhagen, Denmark; Cancer Systems Biology, Section for Bioinformatics, Department of Health and Technology, Technical University of Denmark, Lyngby, Denmark
| | - Thomas M Frimurer
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Thue W Schwartz
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
11
|
da Rocha MN, da Fonseca AM, Dantas ANM, Dos Santos HS, Marinho ES, Marinho GS. In Silico Study in MPO and Molecular Docking of the Synthetic Drynaran Analogues Against the Chronic Tinnitus: Modulation of the M1 Muscarinic Acetylcholine Receptor. Mol Biotechnol 2024; 66:254-269. [PMID: 37079267 DOI: 10.1007/s12033-023-00748-5] [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] [Received: 02/17/2023] [Accepted: 04/03/2023] [Indexed: 04/21/2023]
Abstract
Tinnitus is a syndrome that affects the human auditory system and is characterized by a perception of sounds in the absence of acoustic stimuli, or in total silence. Research indicates that muscarinic acetylcholine receptors (mAChRs), especially the M1 type, have a fundamental role in the alterations of auditory perceptions of tinnitus. Here, a series of computer-aided tools were used, from molecular surface analysis software to services available on the web for estimating pharmacokinetics and pharmacodynamics. The results infer that the low lipophilicity ligands, that is, the 1a-d alkyl furans, present the best pharmacokinetic profile, as compounds with an optimal alignment between permeability and clearance. However, only ligands 1a and 1b have properties that are safe for the central nervous system, the site of cholinergic modulation. These ligands showed similarity with compounds deposited in the European Molecular Biology Laboratory chemical (ChEMBL) database acting on the mAChRs M1 type, the target selected for the molecular docking test. The simulations suggest that the 1 g ligand can form the ligand-receptor complex with the best affinity energy order and that, together with the 1b ligand, they are competitive agonists in relation to the antagonist Tiotropium, in addition to acting in synergism with the drug Bromazepam in the treatment of chronic tinnitus.
Collapse
Affiliation(s)
- Matheus Nunes da Rocha
- Graduate Program in Natural Sciences, Center for Science and Technology, State University of Ceará, Fortaleza, CE, Brazil.
| | - Aluísio Marques da Fonseca
- Institute of Engineering and Sustainable Development, Academic Master in Sociobiodiversity and Sustainable Technologies, University of International Integration of Afro-Brazilian Lusofonia, Acarape, CE, Brazil
| | | | | | - Emmanuel Silva Marinho
- Graduate Program in Natural Sciences, Center for Science and Technology, State University of Ceará, Fortaleza, CE, Brazil
- Group of Theoretical Chemistry and Electrochemistry, State University of Ceará, Limoeiro Do Norte, CE, Brazil
| | - Gabrielle Silva Marinho
- Group of Theoretical Chemistry and Electrochemistry, State University of Ceará, Limoeiro Do Norte, CE, Brazil
| |
Collapse
|
12
|
AlRawashdeh S, Barakat KH. Applications of Molecular Dynamics Simulations in Drug Discovery. Methods Mol Biol 2024; 2714:127-141. [PMID: 37676596 DOI: 10.1007/978-1-0716-3441-7_7] [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
In the current drug development process, molecular dynamics (MD) simulations have proven to be very useful. This chapter provides an overview of the current applications of MD simulations in drug discovery, from detecting protein druggable sites and validating drug docking outcomes to exploring protein conformations and investigating the influence of mutations on its structure and functions. In addition, this chapter emphasizes various strategies to improve the conformational sampling efficiency in molecular dynamics simulations. With a growing computer power and developments in the production of force fields and MD techniques, the importance of MD simulations in helping the drug development process is projected to rise significantly in the future.
Collapse
Affiliation(s)
- Sara AlRawashdeh
- Faculty of Pharmacy and Pharmaceutical Sciences, University of Alberta, Edmonton, AB, Canada
| | - Khaled H Barakat
- Faculty of Pharmacy and Pharmaceutical Sciences, University of Alberta, Edmonton, AB, Canada.
| |
Collapse
|
13
|
Burger WAC, Pham V, Vuckovic Z, Powers AS, Mobbs JI, Laloudakis Y, Glukhova A, Wootten D, Tobin AB, Sexton PM, Paul SM, Felder CC, Danev R, Dror RO, Christopoulos A, Valant C, Thal DM. Xanomeline displays concomitant orthosteric and allosteric binding modes at the M 4 mAChR. Nat Commun 2023; 14:5440. [PMID: 37673901 PMCID: PMC10482975 DOI: 10.1038/s41467-023-41199-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Accepted: 08/26/2023] [Indexed: 09/08/2023] Open
Abstract
The M4 muscarinic acetylcholine receptor (M4 mAChR) has emerged as a drug target of high therapeutic interest due to its expression in regions of the brain involved in the regulation of psychosis, cognition, and addiction. The mAChR agonist, xanomeline, has provided significant improvement in the Positive and Negative Symptom Scale (PANSS) scores in a Phase II clinical trial for the treatment of patients suffering from schizophrenia. Here we report the active state cryo-EM structure of xanomeline bound to the human M4 mAChR in complex with the heterotrimeric Gi1 transducer protein. Unexpectedly, two molecules of xanomeline were found to concomitantly bind to the monomeric M4 mAChR, with one molecule bound in the orthosteric (acetylcholine-binding) site and a second molecule in an extracellular vestibular allosteric site. Molecular dynamic simulations supports the structural findings, and pharmacological validation confirmed that xanomeline acts as a dual orthosteric and allosteric ligand at the human M4 mAChR. These findings provide a basis for further understanding xanomeline's complex pharmacology and highlight the myriad of ways through which clinically relevant ligands can bind to and regulate GPCRs.
Collapse
Affiliation(s)
- Wessel A C Burger
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
- Australian Research Council Centre for Cryo-Electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | - Vi Pham
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | - Ziva Vuckovic
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | - Alexander S Powers
- Department of Chemistry, Stanford University, Stanford, CA, 94305, USA
- Departments of Computer Science, Structural Biology, and Molecular and Cellular Physiology, Stanford University, Stanford, CA, 94305, USA
| | - Jesse I Mobbs
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
- Australian Research Council Centre for Cryo-Electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | - Yianni Laloudakis
- Department of Chemistry, Stanford University, Stanford, CA, 94305, USA
| | - Alisa Glukhova
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
- Australian Research Council Centre for Cryo-Electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | - Denise Wootten
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
- Australian Research Council Centre for Cryo-Electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | - Andrew B Tobin
- The Advanced Research Centre (ARC), Centre for Translational Science, School of Biomolecular Sciences, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, G12 8QQ, UK
| | - Patrick M Sexton
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
- Australian Research Council Centre for Cryo-Electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | | | | | - Radostin Danev
- Graduate School of Medicine, University of Tokyo, N415, 7-3-1 Hongo, Bunkyo-ku, 113-0033, Tokyo, Japan
| | - Ron O Dror
- Departments of Computer Science, Structural Biology, and Molecular and Cellular Physiology, Stanford University, Stanford, CA, 94305, USA.
| | - Arthur Christopoulos
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia.
- Australian Research Council Centre for Cryo-Electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia.
- Neuromedicines Discovery Centre, Monash University, Parkville, VIC, 3052, Australia.
| | - Celine Valant
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia.
| | - David M Thal
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia.
- Australian Research Council Centre for Cryo-Electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia.
| |
Collapse
|
14
|
Galvani F, Pala D, Cuzzolin A, Scalvini L, Lodola A, Mor M, Rizzi A. Unbinding Kinetics of Muscarinic M3 Receptor Antagonists Explained by Metadynamics Simulations. J Chem Inf Model 2023; 63:2842-2856. [PMID: 37053454 PMCID: PMC10170513 DOI: 10.1021/acs.jcim.3c00042] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2023] [Indexed: 04/15/2023]
Abstract
The residence time (RT), the time for which a drug remains bound to its biological target, is a critical parameter for drug design. The prediction of this key kinetic property has been proven to be challenging and computationally demanding in the framework of atomistic simulations. In the present work, we setup and applied two distinct metadynamics protocols to estimate the RTs of muscarinic M3 receptor antagonists. In the first method, derived from the conformational flooding approach, the kinetics of unbinding is retrieved from a physics-based parameter known as the acceleration factor α (i.e., the running average over time of the potential deposited in the bound state). Such an approach is expected to recover the absolute RT value for a compound of interest. In the second method, known as the tMETA-D approach, a qualitative estimation of the RT is given by the time of simulation required to drive the ligand from the binding site to the solvent bulk. This approach has been developed to reproduce the change of experimental RTs for compounds targeting the same target. Our analysis shows that both computational protocols are able to rank compounds in agreement with their experimental RTs. Quantitative structure-kinetics relationship (SKR) models can be identified and employed to predict the impact of a chemical modification on the experimental RT once a calibration study has been performed.
Collapse
Affiliation(s)
- 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
| | - Daniele Pala
- Chemistry
Research and Drug Design Department, Chiesi
Farmaceutici S.p.A., Largo F. Belloli 11/A, 43122 Parma, Italy
| | - Alberto Cuzzolin
- Chemistry
Research and Drug Design Department, Chiesi
Farmaceutici S.p.A., Largo F. Belloli 11/A, 43122 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
| | - 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
| | - 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, University of Parma, Parco Area delle Scienze 11/A, I-43124 Parma, Italy
| | - Andrea Rizzi
- Chemistry
Research and Drug Design Department, Chiesi
Farmaceutici S.p.A., Largo F. Belloli 11/A, 43122 Parma, Italy
| |
Collapse
|
15
|
Wang J, Do HN, Koirala K, Miao Y. Predicting Biomolecular Binding Kinetics: A Review. J Chem Theory Comput 2023; 19:2135-2148. [PMID: 36989090 DOI: 10.1021/acs.jctc.2c01085] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
Biomolecular binding kinetics including the association (kon) and dissociation (koff) rates are critical parameters for therapeutic design of small-molecule drugs, peptides, and antibodies. Notably, the drug molecule residence time or dissociation rate has been shown to correlate with their efficacies better than binding affinities. A wide range of modeling approaches including quantitative structure-kinetic relationship models, Molecular Dynamics simulations, enhanced sampling, and Machine Learning has been developed to explore biomolecular binding and dissociation mechanisms and predict binding kinetic rates. Here, we review recent advances in computational modeling of biomolecular binding kinetics, with an outlook for future improvements.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Hung N Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Kushal Koirala
- 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
|
16
|
Ligand Gaussian Accelerated Molecular Dynamics 2 (LiGaMD2): Improved Calculations of Ligand Binding Thermodynamics and Kinetics with Closed Protein Pocket. J Chem Theory Comput 2023; 19:733-745. [PMID: 36706316 DOI: 10.1021/acs.jctc.2c01194] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Ligand binding thermodynamics and kinetics are critical parameters for drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics from molecular simulations due to limited simulation timescales. Protein dynamics, especially in the ligand binding pocket, often plays an important role in ligand binding. Based on our previously developed Ligand Gaussian accelerated molecular dynamics (LiGaMD), here we present LiGaMD2 in which a selective boost potential was applied to both the ligand and protein residues in the binding pocket to improve sampling of ligand binding and dissociation. To validate the performance of LiGaMD2, the T4 lysozyme (T4L) mutants with open and closed pockets bound by different ligands were chosen as model systems. LiGaMD2 could efficiently capture repetitive ligand dissociation and binding within microsecond simulations of all T4L systems. The obtained ligand binding kinetic rates and free energies agreed well with available experimental values and previous modeling results. Therefore, LiGaMD2 provides an improved approach to sample opening of closed protein pockets for ligand dissociation and binding, thereby allowing for efficient calculations of ligand binding thermodynamics and kinetics.
Collapse
|
17
|
Wu F, Jin S, Jiang Y, Jin X, Tang B, Niu Z, Liu X, Zhang Q, Zeng X, Li SZ. Pre-Training of Equivariant Graph Matching Networks with Conformation Flexibility for Drug Binding. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2022; 9:e2203796. [PMID: 36202759 PMCID: PMC9685463 DOI: 10.1002/advs.202203796] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Revised: 09/07/2022] [Indexed: 05/16/2023]
Abstract
The latest biological findings observe that the motionless "lock-and-key" theory is not generally applicable and that changes in atomic sites and binding pose can provide important information for understanding drug binding. However, the computational expenditure limits the growth of protein trajectory-related studies, thus hindering the possibility of supervised learning. A spatial-temporal pre-training method based on the modified equivariant graph matching networks, dubbed ProtMD which has two specially designed self-supervised learning tasks: atom-level prompt-based denoising generative task and conformation-level snapshot ordering task to seize the flexibility information inside molecular dynamics (MD) trajectories with very fine temporal resolutions is presented. The ProtMD can grant the encoder network the capacity to capture the time-dependent geometric mobility of conformations along MD trajectories. Two downstream tasks are chosen to verify the effectiveness of ProtMD through linear detection and task-specific fine-tuning. A huge improvement from current state-of-the-art methods, with a decrease of 4.3% in root mean square error for the binding affinity problem and an average increase of 13.8% in the area under receiver operating characteristic curve and the area under the precision-recall curve for the ligand efficacy problem is observed. The results demonstrate a strong correlation between the magnitude of conformation's motion in the 3D space and the strength with which the ligand binds with its receptor.
Collapse
Affiliation(s)
- Fang Wu
- School of EngineeringWestlake UniversityHangzhou310024China
- MindRank AI Ltd.Hangzhou310000China
| | - Shuting Jin
- MindRank AI Ltd.Hangzhou310000China
- School of InformaticsXiamen UniversityXiamen361005China
| | | | | | | | | | - Xiangrong Liu
- School of InformaticsXiamen UniversityXiamen361005China
| | - Qiang Zhang
- ZJU‐Hangzhou Global Scientific and Technological Innovation CenterHangzhou311200China
- College of Computer Science and TechnologyZhejiang UniversityHangzhou310013China
| | - Xiangxiang Zeng
- School of Information Science and EngineeringHunan UniversityHunan410082China
| | - Stan Z. Li
- School of EngineeringWestlake UniversityHangzhou310024China
| |
Collapse
|
18
|
Khan SH, Braet SM, Koehler SJ, Elacqua E, Anand GS, Okafor CD. Ligand-induced shifts in conformational ensembles that describe transcriptional activation. eLife 2022; 11:e80140. [PMID: 36222302 PMCID: PMC9555869 DOI: 10.7554/elife.80140] [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: 05/10/2022] [Accepted: 09/14/2022] [Indexed: 11/15/2022] Open
Abstract
Nuclear receptors function as ligand-regulated transcription factors whose ability to regulate diverse physiological processes is closely linked with conformational changes induced upon ligand binding. Understanding how conformational populations of nuclear receptors are shifted by various ligands could illuminate strategies for the design of synthetic modulators to regulate specific transcriptional programs. Here, we investigate ligand-induced conformational changes using a reconstructed, ancestral nuclear receptor. By making substitutions at a key position, we engineer receptor variants with altered ligand specificities. We combine cellular and biophysical experiments to characterize transcriptional activity, as well as elucidate mechanisms underlying altered transcription in receptor variants. We then use atomistic molecular dynamics (MD) simulations with enhanced sampling to generate ensembles of wildtype and engineered receptors in combination with multiple ligands, followed by conformational analysis and correlation of MD-based predictions with functional ligand profiles. We determine that conformational ensembles accurately describe ligand responses based on observed population shifts. These studies provide a platform which will allow structural characterization of physiologically-relevant conformational ensembles, as well as provide the ability to design and predict transcriptional responses in novel ligands.
Collapse
Affiliation(s)
- Sabab Hasan Khan
- Department of Biochemistry and Molecular Biology, Pennsylvania State UniversityState CollegeUnited States
| | - Sean M Braet
- Department of Chemistry, Pennsylvania State UniversityState ParkUnited States
| | | | - Elizabeth Elacqua
- Department of Chemistry, Pennsylvania State UniversityState ParkUnited States
| | | | - C Denise Okafor
- Department of Biochemistry and Molecular Biology, Pennsylvania State UniversityState CollegeUnited States
- Department of Chemistry, Pennsylvania State UniversityState ParkUnited States
| |
Collapse
|
19
|
Yang Z, Zang Y, Wang H, Kang Y, Zhang J, Li X, Zhang L, Zhang S. Recognition between CD147 and cyclophilin A deciphered by accelerated molecular dynamics simulations. Phys Chem Chem Phys 2022; 24:18905-18914. [PMID: 35913096 DOI: 10.1039/d2cp01975b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
CD147 functions as the receptor of extracellular cyclophilin A (CypA) in various diseases, and CD147-CypA binding ulteriorly underlies the pathological process of various viral infections including HIV-1, SARS, and SARS-CoV-2. Although CyPA has been identified as a key intermediate pro-inflammatory factor, the mechanism by which CD147 cooperates with CypA in the development of the cytokine storm remains largely unknown, and the binding profile of CD147 with CypA remains to be elucidated as well. Here, we prepared three binding models of the CD147-CypA complex, including the active site of CypA severally binding to the groove bound by the Ig1 and Ig2 domains (model-0), P180-G181 (model-1), and P211 (model-2) of CD147, as well as introducing mutations P180A-G181A and P211A individually in each model. All systems were studied using accelerated molecular dynamics simulations and the molecular mechanics generalized Born surface area (MM/GBSA) method. For model-0, CypA bound to the ectodomain of CD147 with the highest binding affinity. Moreover, mutations P180A-G181A of CD147 in model-0 decreased the binding affinity and weakened the dynamic correlation between CD147 and CypA, which resulted in CypA shifting from the initial binding location. Other residue mutations of CD147 did not significantly affect the CD147-CypA binding, as reflected by the energy and structural analyses. Compared with surface plasmon resonance results and nuclear magnetic resonance shift signals, CypA should tend to reciprocally bind to the groove of CD147, and the binding process might be modulated by P180-G181 rather than P211. Besides, residue R201 of CD147 is critical for CD147-CypA binding and needs further experimental verification. These findings further our understanding of the recruitment between CD147 and CypA and its potential role in the development of inflammation and viral infection.
Collapse
Affiliation(s)
- Zhiwei Yang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Yongjian Zang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - He Wang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Ying Kang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Jianwen Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Xuhua Li
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Lei Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Shengli Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| |
Collapse
|
20
|
Dragan P, Atzei A, Sanmukh SG, Latek D. Computational and experimental approaches to probe GPCR activation and signaling. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2022; 193:1-36. [PMID: 36357073 DOI: 10.1016/bs.pmbts.2022.06.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.
Collapse
Affiliation(s)
- Paulina Dragan
- Faculty of Chemistry, University of Warsaw, Warsaw, Poland
| | | | | | - Dorota Latek
- Faculty of Chemistry, University of Warsaw, Warsaw, Poland.
| |
Collapse
|
21
|
Zang Y, Wang H, Kang Y, Zhang J, Li X, Zhang L, Yang Z, Zhang S. TAB1 binding induced p38α conformation change: an accelerated molecular dynamics simulation study. Phys Chem Chem Phys 2022; 24:10506-10513. [PMID: 35441632 DOI: 10.1039/d2cp00144f] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
p38α mitogen-activated protein kinase (MAPK) undergoes autophosphorylation induced by the binding of TGFβ-activated kinase 1 binding protein 1 (TAB1) in myocardial ischemia. Investigation of the conformational transformations in p38α triggered by TAB1 binding is motivated by the need to find selective p38α activation inhibitors to treat myocardial ischemia. Herein, the conformational transformations of p38α were studied via all-atom accelerated molecular dynamics simulations and principal component analysis. With the binding of TAB1, the conformational changes of p38α auto-activation were characterized by the movement of the activation loop (A-loop) away from the αG helix toward the αF, αE helixes and L16-loop. In addition, a diverse intermediate state with an extensional and phosphorylated A-loop different from the transition intermediate state was explored. The conformational changes, including the A-loop alpha-structure breaking and the stronger hydrogen bond network formation, are accompanied by the extension of the A-loop and more intramolecular interactions in p38α. TAB1 correlates with other regions of p38α that are distal from the TAB1-binding site, including the A-loop, αC helix, and L16-loop, which regulates the intramolecular correlation of p38α. And, the phosphorylation further enhances the correlations between the A-loop and the other regions of p38α. The correlation results imply the regulation process of p38α conformational transformations. These findings will improve our understanding of the autophosphorylation of kinase and facilitate the development of selective inhibitors for the treatment of ischemic injury.
Collapse
Affiliation(s)
- Yongjian Zang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - He Wang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Ying Kang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Jianwen Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Xuhua Li
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Lei Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Zhiwei Yang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Shengli Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| |
Collapse
|
22
|
Liao J, Nie X, Unarta IC, Ericksen SS, Tang W. In Silico Modeling and Scoring of PROTAC-Mediated Ternary Complex Poses. J Med Chem 2022; 65:6116-6132. [PMID: 35412837 DOI: 10.1021/acs.jmedchem.1c02155] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Proteolysis targeting chimeras (PROTACs) are molecules that induce protein degradation via formation of ternary complexes between an E3 ubiquitin ligase and a target protein. The rational design of PROTACs requires accurate knowledge of the native configuration of the PROTAC-induced ternary complex. This study demonstrates that native and non-native ternary complex poses can be distinguished based on the pose occupancy time in MD, where native poses exhibit longer occupancy times at both room and higher temperatures. Candidate poses are generated by MD sampling and pre-ranked by classic MM/GBSA. A specific heating scheme is then applied to accelerate ternary pose departure, with the pose occupancy time and fraction being measured. This scoring identifies the native pose in all systems tested. Its success is partially attributed to the dynamic nature of pose departure analyses, which accounts for entropic effects typically neglected in the faster static scoring methods, while entropy plays a greater role in protein-protein than in protein-ligand systems.
Collapse
Affiliation(s)
- Junzhuo Liao
- School of Pharmacy, University of Wisconsin-Madison, Madison, Wisconsin 53705, United States
| | - Xueqing Nie
- School of Pharmacy, University of Wisconsin-Madison, Madison, Wisconsin 53705, United States
| | - Ilona Christy Unarta
- Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Spencer S Ericksen
- Drug Development Core, UW Carbone Cancer Center, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin 53705, United States
| | - Weiping Tang
- School of Pharmacy, University of Wisconsin-Madison, Madison, Wisconsin 53705, United States.,Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States.,Drug Development Core, UW Carbone Cancer Center, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin 53705, United States
| |
Collapse
|
23
|
Zhang F, Yuan Y, Chen Y, Chen J, Guo Y, Pu X. Molecular insights into the allosteric coupling mechanism between an agonist and two different transducers for μ-opioid receptors. Phys Chem Chem Phys 2022; 24:5282-5293. [PMID: 35170592 DOI: 10.1039/d1cp05736g] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
G protein-coupled receptors (GPCRs) as the most important class of pharmacological targets regulate G-protein and β-arrestin-mediated signaling through allosteric interplay, which are responsible for different biochemical and physiological actions like therapeutic efficacy and side effects. However, the allosteric mechanism underlying preferentially recruiting one transducer versus the other has been poorly understood, limiting drug design. Motivated by this issue, we utilize accelerated molecular dynamics simulation coupled with potential of mean force (PMF), molecular mechanics Poisson Boltzmann surface area (MM/PBSA) and protein structure network (PSN) to study two ternary complex systems of a representative class A GPCR (μ-opioid receptor (μOR)) bound by an agonist and one specific transducer (G-protein and β-arrestin). The results show that no significant difference exists in the whole structure of μOR between two transducer couplings, but displays transducer-dependent changes in the intracellular binding region of μOR, where the β-arrestin coupling results in a narrower crevice with TM7 inward movement compared with the G-protein. In addition, both the G-protein and β-arrestin coupling can increase the binding affinity of the agonist to the receptor. However, the interactions between the agonist and μOR also exhibit transducer-specific changes, in particular for the interaction with ECL2 that plays an important role in recruiting β-arrestin. The allosteric network analysis further indicates that Y1483.33, F1523.37, F1563.41, N1914.49, T1603.45, Y1062.42, W2936.48, F2896.44, I2485.54 and Y2525.58 play important roles in equally activating G-protein and β-arrestin. In contrast, M1613.46 and R1653.50 devote important contributions to preferentially recruit G-protein while D1643.49 and R179ICL2 are revealed to be important for selectively activating β-arrestin. The observations provide useful information for understanding the biased activation mechanism.
Collapse
Affiliation(s)
- Fuhui Zhang
- Faculty of Chemistry, Sichuan University, Chengdu, Sichuan 610064, People's Republic of China.
| | - Yuan Yuan
- College of Management, Southwest University for Nationalities, Chengdu, Sichuan 610041, People's Republic of China
| | - Yichi Chen
- Faculty of Chemistry, Sichuan University, Chengdu, Sichuan 610064, People's Republic of China.
| | - Jianfang Chen
- Faculty of Chemistry, Sichuan University, Chengdu, Sichuan 610064, People's Republic of China.
| | - Yanzhi Guo
- Faculty of Chemistry, Sichuan University, Chengdu, Sichuan 610064, People's Republic of China.
| | - Xuemei Pu
- Faculty of Chemistry, Sichuan University, Chengdu, Sichuan 610064, People's Republic of China.
| |
Collapse
|
24
|
Li C, Liu K, Chen S, Han L, Han W. Gaussian Accelerated Molecular Dynamics Simulations Investigation on the Mechanism of Angiotensin-Converting Enzyme (ACE) C-Domain Inhibition by Dipeptides. Foods 2022; 11:foods11030327. [PMID: 35159478 PMCID: PMC8834632 DOI: 10.3390/foods11030327] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Revised: 12/29/2021] [Accepted: 01/01/2022] [Indexed: 02/06/2023] Open
Abstract
Angiotensin-converting enzyme (ACE)-inhibitory peptides extracted from food proteins can lower blood pressure by inhibiting ACE activity. A recent study showed that the inhibitory activity of IY (Ile-Tyr, a dipeptide derived from soybean protein) against ACE was much higher than that of LL (Leu-Leu), although they had similar hydrophobic and predicted activity values. It was difficult to reveal the deep molecular mechanism underlying this phenomenon by traditional experimental methods. The Apo and two complex systems (i.e., ACE–LL and ACE–IY) were therefore subjected to 1 μs long Gaussian accelerated molecular dynamics (GaMD) simulations. The results showed that the binding of IY can cause obvious contraction of the active site of ACE, mainly manifested by a significant lateral shift of α13, α14, and α15. In addition, hinge 2 and hinge 3 were more stable in the ACE–IY system, while these phenomena were not present in the ACE–LL system. Moreover, the α10 of the IY-bound ACE kept an inward state during the simulation progress, which facilitated the ACE to remain closed. However, for the LL-bound ACE, the α10 switched between two outward states. To sum up, our study provides detailed insights into inhibitor-induced conformational changes in ACE that may help in the design of specific inhibitors targeting ACE for the treatment of hypertension.
Collapse
|
25
|
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
|
26
|
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
|
27
|
Basciu A, Callea L, Motta S, Bonvin AM, Bonati L, Vargiu AV. No dance, no partner! A tale of receptor flexibility in docking and virtual screening. VIRTUAL SCREENING AND DRUG DOCKING 2022. [DOI: 10.1016/bs.armc.2022.08.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
28
|
Adelusi TI, Oyedele AQK, Boyenle ID, Ogunlana AT, Adeyemi RO, Ukachi CD, Idris MO, Olaoba OT, Adedotun IO, Kolawole OE, Xiaoxing Y, Abdul-Hammed M. Molecular modeling in drug discovery. INFORMATICS IN MEDICINE UNLOCKED 2022. [DOI: 10.1016/j.imu.2022.100880] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
|
29
|
Hernández González JE, Alberca LN, Masforrol González Y, Reyes Acosta O, Talevi A, Salas-Sarduy E. Tetracycline Derivatives Inhibit Plasmodial Cysteine Protease Falcipain-2 through Binding to a Distal Allosteric Site. J Chem Inf Model 2021; 62:159-175. [PMID: 34962803 DOI: 10.1021/acs.jcim.1c01189] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Allosteric inhibitors regulate enzyme activity from remote and usually specific pockets. As they promise an avenue for less toxic and safer drugs, the identification and characterization of allosteric inhibitors has gained great academic and biomedical interest in recent years. Research on falcipain-2 (FP-2), the major papain-like cysteine hemoglobinase of Plasmodium falciparum, might benefit from this strategy to overcome the low selectivity against human cathepsins shown by active site-directed inhibitors. Encouraged by our previous finding that methacycline inhibits FP-2 noncompetitively, here we assessed other five tetracycline derivatives against this target and characterized their inhibition mechanism. As previously shown for methacycline, tetracycline derivatives inhibited FP-2 in a noncompetitive fashion, with Ki values ranging from 121 to 190 μM. A possible binding to the S' side of the FP-2 active site, similar to that described by X-ray crystallography (PDB: 6SSZ) for the noncompetitive inhibitor E-chalcone 48 (EC48), was experimentally discarded by kinetic analysis using a large peptidyl substrate spanning the whole active site. By combining lengthy molecular dynamics (MD) simulations that allowed methacycline to diffuse from solution to different FP-2 surface regions and free energy calculations, we predicted the most likely binding mode of the ligand. Of note, the proposed binding pose explains the low differences in Ki values observed for the tested tetracycline derivatives and the calculated binding free energies match the experimental values. Overall, this study has implications for the design of novel allosteric inhibitors against FP-2 and sets the basis for further optimization of the tetracycline scaffold to produce more potent and selective inhibitors.
Collapse
Affiliation(s)
- Jorge Enrique Hernández González
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista Júlio de Mesquita Filho, Rua Cristóvão Colombo, 2265, Jardim Nazareth, São José do Rio Preto, São Paulo CEP 15054-000, Brazil
| | - Lucas N Alberca
- Laboratory of Bioactive Compounds Research and Development (LIDeB), Department of Biological Sciences, Exact Sciences College, Universidad Nacional de La Plata, La Plata B1900ADU, Argentina
| | | | - Osvaldo Reyes Acosta
- Chemistry and Physics Department, Center for Genetic Engineering and Biotechnology, Havana 10600, Cuba
| | - Alan Talevi
- Laboratory of Bioactive Compounds Research and Development (LIDeB), Department of Biological Sciences, Exact Sciences College, Universidad Nacional de La Plata, La Plata B1900ADU, Argentina
| | - Emir Salas-Sarduy
- Instituto de Investigaciones Biotecnológicas "Dr. Rodolfo Ugalde"─Universidad Nacional de San Martín─CONICET, San Martín B1650HMP, Buenos Aires, Argentina
| |
Collapse
|
30
|
Abrol R, Serrano E, Santiago LJ. Development of enhanced conformational sampling methods to probe the activation landscape of GPCRs. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2021; 128:325-359. [PMID: 35034722 DOI: 10.1016/bs.apcsb.2021.11.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
G protein-coupled receptors (GPCRs) make up the largest superfamily of integral membrane proteins and play critical signal transduction roles in many physiological processes. Developments in molecular biology, biophysical, biochemical, pharmacological, and computational techniques aimed at these important therapeutic targets are beginning to provide unprecedented details on the structural as well as functional basis of their pleiotropic signaling mediated by G proteins, β arrestins, and other transducers. This pleiotropy presents a pharmacological challenge as the same ligand-receptor interaction can cause a therapeutic effect as well as an undesirable on-target side-effect through different downstream pathways. GPCRs don't function as simple binary on-off switches but as finely tuned shape-shifting machines described by conformational ensembles, where unique subsets of conformations may be responsible for specific signaling cascades. X-ray crystallography and more recently cryo-electron microscopy are providing snapshots of some of these functionally-important receptor conformations bound to ligands and/or transducers, which are being utilized by computational methods to describe the dynamic conformational energy landscape of GPCRs. In this chapter, we review the progress in computational conformational sampling methods based on molecular dynamics and discrete sampling approaches that have been successful in complementing biophysical and biochemical studies on these receptors in terms of their activation mechanisms, allosteric effects, actions of biased ligands, and effects of pathological mutations. Some of the sampled simulation time scales are beginning to approach receptor activation time scales. The list of conformational sampling methods and example uses discussed is not exhaustive but includes representative examples that have pushed the limits of classical molecular dynamics and discrete sampling methods to describe the activation energy landscape of GPCRs.
Collapse
Affiliation(s)
- Ravinder Abrol
- Department of Chemistry and Biochemistry, California State University, Northridge, CA, United States.
| | - Erik Serrano
- Department of Chemistry and Biochemistry, California State University, Northridge, CA, United States
| | - Luis Jaimes Santiago
- Department of Chemistry and Biochemistry, California State University, Northridge, CA, United States
| |
Collapse
|
31
|
Zhang F, Chen X, Chen J, Xu Y, Li S, Guo Y, Pu X. Probing Allosteric Regulation Mechanism of W7.35 on Agonist-Induced Activity for μOR by Mutation Simulation. J Chem Inf Model 2021; 62:5120-5135. [PMID: 34779608 DOI: 10.1021/acs.jcim.1c00650] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The residue located at 15 positions before the most conserved residue in TM7 (7.35 of Ballesteros-Weinstein number) plays important roles in ligand binding and the receptor activity for class A GPCRs. Nevertheless, its regulation mechanism has not been clearly clarified in experiments, and some controversies also exist for its impact on μ-opioid receptors (μOR) bound by agonists. Thus, we chose the μ-opioid receptor (μOR) of class A GPCRs as a representative and utilized a microsecond accelerated molecular dynamics simulation (aMD) coupled with a protein structure network (PSN) to explore the effect of W3187.35 on its functional activity induced by the agonist endomorphin2 mainly by a comparison of the wild system and its W7.35A mutant. When endomorphin2 binds to the wild-type μOR, TM6 in μOR moves outward to form an open intracellular conformation that is beneficial to accommodating the β-arrestin transducer, rather than the G-protein transducer due to the clash with the α5 helix of G-protein, thus acting as a β-arrestin biased agonist. However, the W318A mutation induces the intracellular part of μOR to form a closed state, which disfavors coupling with either G-protein or β-arrestin. The allosteric pathway analysis further reveals that the binding of endomorphin2 to the wild-type μOR transmits more activation signals to the β-arrestin binding site while the W318A mutation induces more structural signals to transmit to common binding residues of the G protein and β-arrestin. More interestingly, the residue at the 7.35 position regulates the shortest allosteric pathway in indirect ways by influencing the interactions between other ligand-binding residues and endomorphin2. W2936.48 and F2896.44 are important for regulating the different activities of μOR induced either by the agonist or by the mutation. Y3367.53, F3438.50, and D3408.47 play crucial roles in activating the β-arrestin biased signal induced by the agonist endomorphin2, while L1583.43 and V2866.41 devote important contributions to the change in the activity of endomorphin2 from the β-arrestin biased agonist to the antagonist upon the W318A mutation.
Collapse
Affiliation(s)
- Fuhui Zhang
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Xin Chen
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Jianfang Chen
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Yanjiani Xu
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Shiqi Li
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Yanzhi Guo
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Xuemei Pu
- College of Chemistry, Sichuan University, Chengdu 610064, China
| |
Collapse
|
32
|
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
|
33
|
Conformational dynamics of androgen receptors bound to agonists and antagonists. Sci Rep 2021; 11:15887. [PMID: 34354111 PMCID: PMC8342701 DOI: 10.1038/s41598-021-94707-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2021] [Accepted: 07/13/2021] [Indexed: 11/09/2022] Open
Abstract
The androgen receptor (AR) is critical in the progression of prostate cancer (PCa). Small molecule antagonists that bind to the ligand binding domain (LBD) of the AR have been successful in treating PCa. However, the structural basis by which the AR antagonists manifest their therapeutic efficacy remains unclear, due to the lack of detailed structural information of the AR bound to the antagonists. We have performed accelerated molecular dynamics (aMD) simulations of LBDs bound to a set of ligands including a natural substrate (dihydrotestosterone), an agonist (RU59063) and three antagonists (bicalutamide, enzalutamide and apalutamide) as well as in the absence of ligand (apo). We show that the binding of AR antagonists at the substrate binding pocket alter the dynamic fluctuations of H12, thereby disrupting the structural integrity of the agonistic conformation of AR. Two antagonists, enzalutamide and apalutamide, induce considerable structural changes to the agonist conformation of LBD, when bound close to H12 of AR LBD. When the antagonists bind to the pocket with different orientations having close contact with H11, no significant conformational changes were observed, suggesting the AR remains in the functionally activated (agonistic) state. The simulations on a drug resistance mutant F876L bound to enzalutamide demonstrated that the mutation stabilizes the agonistic conformation of AR LBD, which compromises the efficacy of the antagonists. Principal component analysis (PCA) of the structural fluctuations shows that the binding of enzalutamide and apalutamide induce conformational fluctuations in the AR, which are markedly different from those caused by the agonist as well as another antagonist, bicalutamide. These fluctuations could only be observed with the use of aMD.
Collapse
|
34
|
Bhakat S. Pepsin-like aspartic proteases (PAPs) as model systems for combining biomolecular simulation with biophysical experiments. RSC Adv 2021; 11:11026-11047. [PMID: 35423571 PMCID: PMC8695779 DOI: 10.1039/d0ra10359d] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2020] [Accepted: 02/21/2021] [Indexed: 01/26/2023] Open
Abstract
Pepsin-like aspartic proteases (PAPs) are a class of aspartic proteases which shares tremendous structural similarity with human pepsin. One of the key structural features of PAPs is the presence of a β-hairpin motif otherwise known as flap. The biological function of the PAPs is highly dependent on the conformational dynamics of the flap region. In apo PAPs, the conformational dynamics of the flap is dominated by the rotational degrees of freedom associated with χ1 and χ2 angles of conserved Tyr (or Phe in some cases). However it is plausible that dihedral order parameters associated with several other residues might play crucial roles in the conformational dynamics of apo PAPs. Due to their size, complexities associated with conformational dynamics and clinical significance (drug targets for malaria, Alzheimer's disease etc.), PAPs provide a challenging testing ground for computational and experimental methods focusing on understanding conformational dynamics and molecular recognition in biomolecules. The opening of the flap region is necessary to accommodate substrate/ligand in the active site of the PAPs. The BIG challenge is to gain atomistic details into how reversible ligand binding/unbinding (molecular recognition) affects the conformational dynamics. Recent reports of kinetics (K i, K d) and thermodynamic parameters (ΔH, TΔS, and ΔG) associated with macro-cyclic ligands bound to BACE1 (belongs to PAP family) provide a perfect challenge (how to deal with big ligands with multiple torsional angles and select optimum order parameters to study reversible ligand binding/unbinding) for computational methods to predict binding free energies and kinetics beyond typical test systems e.g. benzamide-trypsin. In this work, i reviewed several order parameters which were proposed to capture the conformational dynamics and molecular recognition in PAPs. I further highlighted how machine learning methods can be used as order parameters in the context of PAPs. I then proposed some open ideas and challenges in the context of molecular simulation and put forward my case on how biophysical experiments e.g. NMR, time-resolved FRET etc. can be used in conjunction with biomolecular simulation to gain complete atomistic insights into the conformational dynamics of PAPs.
Collapse
Affiliation(s)
- Soumendranath Bhakat
- Division of Biophysical Chemistry, Center for Molecular Protein Science, Department of Chemistry, Lund University P. O. Box 124 SE-22100 Lund Sweden +46-769608418
| |
Collapse
|
35
|
Li C, Chen S, Huang T, Zhang F, Yuan J, Chang H, Li W, Han W. Conformational Changes of Glutamine 5'-Phosphoribosylpyrophosphate Amidotransferase for Two Substrates Analogue Binding: Insight from Conventional Molecular Dynamics and Accelerated Molecular Dynamics Simulations. Front Chem 2021; 9:640994. [PMID: 33718330 PMCID: PMC7953260 DOI: 10.3389/fchem.2021.640994] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Accepted: 01/19/2021] [Indexed: 12/31/2022] Open
Abstract
Glutamine 5′-phosphoribosylpyrophosphate amidotransferase (GPATase) catalyzes the synthesis of phosphoribosylamine, pyrophosphate, and glutamate from phosphoribosylpyrophosphate, as well as glutamine at two sites (i.e., glutaminase and phosphoribosylpyrophosphate sites), through a 20 Å NH3 channel. In this study, conventional molecular dynamics (cMD) simulations and enhanced sampling accelerated molecular dynamics (aMD) simulations were integrated to characterize the mechanism for coordination catalysis at two separate active sites in the enzyme. Results of cMD simulations illustrated the mechanism by which two substrate analogues, namely, DON and cPRPP, affect the structural stability of GPATase from the perspective of dynamic behavior. aMD simulations obtained several key findings. First, a comparison of protein conformational changes in the complexes of GPATase–DON and GPATase–DON–cPRPP showed that binding cPRPP to the PRTase flexible loop (K326 to L350) substantially effected the formation of the R73-DON salt bridge. Moreover, only the PRTase flexible loop in the GPATase–DON–cPRPP complex could remain closed and had sufficient space for cPRPP binding, indicating that binding of DON to the glutamine loop had an impact on the PRTase flexible loop. Finally, both DON and cPRPP tightly bonded to the two domains, thereby inducing the glutamine loop and the PRTase flexible loop to move close to each other. This movement facilitated the transfer of NH3 via the NH3 channel. These theoretical results are useful to the ongoing research on efficient inhibitors related to GPATase.
Collapse
Affiliation(s)
- Congcong Li
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Engineering Laboratory for AIDS Vaccine, School of Life Science, Jilin University, Changchun, China
| | - Siao Chen
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Engineering Laboratory for AIDS Vaccine, School of Life Science, Jilin University, Changchun, China
| | - Tianci Huang
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Engineering Laboratory for AIDS Vaccine, School of Life Science, Jilin University, Changchun, China
| | - Fangning Zhang
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Engineering Laboratory for AIDS Vaccine, School of Life Science, Jilin University, Changchun, China
| | - Jiawei Yuan
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Engineering Laboratory for AIDS Vaccine, School of Life Science, Jilin University, Changchun, China
| | - Hao Chang
- Jilin Province TeyiFood Biotechnology Company Limited, Changchun, China
| | - Wannan Li
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Engineering Laboratory for AIDS Vaccine, School of Life Science, Jilin University, Changchun, China
| | - Weiwei Han
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Engineering Laboratory for AIDS Vaccine, School of Life Science, Jilin University, Changchun, China
| |
Collapse
|
36
|
Hong ST, Su YC, Wang YJ, Cheng TL, Wang YT. Anti-TNF Alpha Antibody Humira with pH-dependent Binding Characteristics: A constant-pH Molecular Dynamics, Gaussian Accelerated Molecular Dynamics, and In Vitro Study. Biomolecules 2021; 11:334. [PMID: 33672169 PMCID: PMC7926962 DOI: 10.3390/biom11020334] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Revised: 02/18/2021] [Accepted: 02/20/2021] [Indexed: 12/17/2022] Open
Abstract
Humira is a monoclonal antibody that binds to TNF alpha, inactivates TNF alpha receptors, and inhibits inflammation. Neonatal Fc receptors can mediate the transcytosis of Humira-TNF alpha complex structures and process them toward degradation pathways, which reduces the therapeutic effect of Humira. Allowing the Humira-TNF alpha complex structures to dissociate to Humira and soluble TNF alpha in the early endosome to enable Humira recycling is crucial. We used the cytoplasmic pH (7.4), the early endosomal pH (6.0), and pKa of histidine side chains (6.0-6.4) to mutate the residues of complementarity-determining regions with histidine. Our engineered Humira (W1-Humira) can bind to TNF alpha in plasma at neutral pH and dissociate from the TNF alpha in the endosome at acidic pH. We used the constant-pH molecular dynamics, Gaussian accelerated molecular dynamics, two-dimensional potential mean force profiles, and in vitro methods to investigate the characteristics of W1-Humira. Our results revealed that the proposed Humira can bind TNF alpha with pH-dependent affinity in vitro. The W1-Humira was weaker than wild-type Humira at neutral pH in vitro, and our prediction results were close to the in vitro results. Furthermore, our approach displayed a high accuracy in antibody pH-dependent binding characteristics prediction, which may facilitate antibody drug design. Advancements in computational methods and computing power may further aid in addressing the challenges in antibody drug design.
Collapse
Affiliation(s)
- Shih-Ting Hong
- Graduate Institute of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan;
| | - Yu-Cheng Su
- Department of Biological Science and Technology, National Yang Ming Chiao Tung University, Hsin-Chu 300, Taiwan;
| | - Yu-Jen Wang
- Department of Mechanical and Electromechanical Engineering, National Sun Yat-sen University, Kaohsiung 80424, Taiwan;
| | - Tian-Lu Cheng
- Graduate Institute of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan;
- Drug Development and Value Creation Research Center, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
| | - Yeng-Tseng Wang
- Graduate Institute of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan;
- Drug Development and Value Creation Research Center, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
- School of Post-Baccalaureate Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
- Department of Medical Research, Kaohsiung Medical University Hospital, Kaohsiung 80708, Taiwan
| |
Collapse
|
37
|
Abstract
Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.
Collapse
|
38
|
Peng C, Wang J, Shi Y, Xu Z, Zhu W. Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis. J Chem Theory Comput 2020; 17:13-28. [PMID: 33351613 DOI: 10.1021/acs.jctc.0c00592] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Understanding conformational change at an atomic level is significant when determining a protein functional mechanism. Replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method to explore protein conformational space. However, REMD with an explicit solvent model requires huge computational resources, immensely limiting its application. In this study, a variation of parallel tempering metadynamics (PTMetaD) with the omission of solvent-solvent interactions in exchange attempts and the use of low-frequency modes calculated by normal-mode analysis (NMA) as collective variables (CVs), namely ossPTMetaD, is proposed with the aim to accelerate MD simulations simultaneously in temperature and geometrical spaces. For testing the performance of ossPTMetaD, five protein systems with diverse biological functions and motion patterns were selected, including large-scale domain motion (AdK), flap movement (HIV-1 protease and BACE1), and DFG-motif flip in kinases (p38α and c-Abl). The simulation results showed that ossPTMetaD requires much fewer numbers of replicas than temperature REMD (T-REMD) with a reduction of ∼70% to achieve a similar exchange ratio. Although it does not obey the detailed balance condition, ossPTMetaD provides consistent results with T-REMD and experimental data. The high accessibility of the large conformational change of protein systems by ossPTMetaD, especially in simulating the very challenging DFG-motif flip of protein kinases, demonstrated its high efficiency and robustness in the characterization of the large-scale protein conformational change pathway and associated free energy profile.
Collapse
Affiliation(s)
- Cheng Peng
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.,University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing, 100049, China
| | - Jinan Wang
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Yulong Shi
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.,University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing, 100049, China
| | - Zhijian Xu
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.,University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing, 100049, China
| | - Weiliang Zhu
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.,Open Studio for Druggability Research of Marine Lead Compounds, Qingdao National Laboratory for Marine Science and Technology, 1 Wenhai Road, Aoshanwei, Jimo, Qingdao 266237, China.,University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing, 100049, China
| |
Collapse
|
39
|
Oliva F, Flores-Canales JC, Pieraccini S, Morelli CF, Sironi M, Schiøtt B. Simulating Multiple Substrate-Binding Events by γ-Glutamyltransferase Using Accelerated Molecular Dynamics. J Phys Chem B 2020; 124:10104-10116. [PMID: 33112625 DOI: 10.1021/acs.jpcb.0c06907] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
γ-Glutamyltransferase (GGT) is an enzyme that uses γ-glutamyl compounds as substrates and catalyzes their transfer to a water molecule or an acceptor substrate with varied physiological function in bacteria, plants, and animals. Crystal structures of GGT are known for different species and in different states of the chemical reaction; however, the structural dynamics of the substrate binding to the catalytic site of GGT are unknown. Here, we modeled Escherichia coli GGT's glutamine binding by using a swarm of accelerated molecular dynamics (aMD) simulations. Characterization of multiple binding events identified three structural binding motifs composed of polar residues in the binding pocket that govern glutamine binding into the active site. Simulated open and closed conformations of a lid-loop protecting the binding cavity suggest its role as a gating element by allowing or blocking substrates entry into the binding pocket. Partially open states of the lid-loop are accessible within thermal fluctuations, while the estimated free energy cost of a complete open state is 2.4 kcal/mol. Our results suggest that both specific electrostatic interactions and GGT conformational dynamics dictate the molecular recognition of substrate-GGT complexes.
Collapse
Affiliation(s)
- Francesco Oliva
- Dipartimento di Chimica, Università degli studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Jose C Flores-Canales
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus, Denmark
| | - Stefano Pieraccini
- Dipartimento di Chimica, Università degli studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Carlo F Morelli
- Dipartimento di Chimica, Università degli studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Maurizio Sironi
- Dipartimento di Chimica, Università degli studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Birgit Schiøtt
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus, Denmark
| |
Collapse
|
40
|
Herrera-Zúñiga LD, Moreno-Vargas LM, Ballaud L, Correa-Basurto J, Prada-Gracia D, Pastré D, Curmi PA, Arrang JM, Maroun RC. Molecular dynamics of the histamine H3 membrane receptor reveals different mechanisms of GPCR signal transduction. Sci Rep 2020; 10:16889. [PMID: 33037273 PMCID: PMC7547658 DOI: 10.1038/s41598-020-73483-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2020] [Accepted: 09/11/2020] [Indexed: 02/02/2023] Open
Abstract
In this work, we studied the mechanisms of classical activation and inactivation of signal transduction by the histamine H3 receptor, a 7-helix transmembrane bundle G-Protein Coupled Receptor through long-time-scale atomistic molecular dynamics simulations of the receptor embedded in a hydrated double layer of dipalmitoyl phosphatidyl choline, a zwitterionic polysaturated ordered lipid. Three systems were prepared: the apo receptor, representing the constitutively active receptor; and two holo-receptors-the receptor coupled to the antagonist/inverse agonist ciproxifan, representing the inactive state of the receptor, and the receptor coupled to the endogenous agonist histamine and representing the active state of the receptor. An extensive analysis of the simulation showed that the three states of H3R present significant structural and dynamical differences as well as a complex behavior given that the measured properties interact in multiple and interdependent ways. In addition, the simulations described an unexpected escape of histamine from the orthosteric binding site, in agreement with the experimental modest affinities and rapid off-rates of agonists.
Collapse
Affiliation(s)
- Leonardo David Herrera-Zúñiga
- UMR-S U1204, Structure et Activité de Biomolécules Normales et Pathologiques, INSERM/Université d'Evry-Val d'Essonne/Université Paris-Saclay, 91000, Evry, France
- Laboratoire de Neurobiologie et Pharmacologie Moléculaire, INSERM U894, Centre de Psychiatrie et Neurosciences, 75014, Paris, France
- Área de Estudios de Posgrado e Investigación, Tecnológico de Estudios Superiores del Oriente del Estado de México, Los Reyes Acaquilpan, Mexico
| | - Liliana Marisol Moreno-Vargas
- Computational Biology and Drug Design Research Unit, Federico Gómez Children's Hospital of Mexico City, Mexico City, Mexico
- Laboratoire de Neurobiologie et Pharmacologie Moléculaire, INSERM U894, Centre de Psychiatrie et Neurosciences, 75014, Paris, France
| | - Luck Ballaud
- Laboratoire de Neurobiologie et Pharmacologie Moléculaire, INSERM U894, Centre de Psychiatrie et Neurosciences, 75014, Paris, France
| | - José Correa-Basurto
- UMR-S U1204, Structure et Activité de Biomolécules Normales et Pathologiques, INSERM/Université d'Evry-Val d'Essonne/Université Paris-Saclay, 91000, Evry, France
- Laboratorio de Modelado Molecular y Bioinformática, Escuela Superior de Medicina, Instituto Politécnico Nacional, Mexico City, Mexico
| | - Diego Prada-Gracia
- Computational Biology and Drug Design Research Unit, Federico Gómez Children's Hospital of Mexico City, Mexico City, Mexico
| | - David Pastré
- UMR-S U1204, Structure et Activité de Biomolécules Normales et Pathologiques, INSERM/Université d'Evry-Val d'Essonne/Université Paris-Saclay, 91000, Evry, France
| | - Patrick A Curmi
- UMR-S U1204, Structure et Activité de Biomolécules Normales et Pathologiques, INSERM/Université d'Evry-Val d'Essonne/Université Paris-Saclay, 91000, Evry, France
| | - Jean Michel Arrang
- Laboratoire de Neurobiologie et Pharmacologie Moléculaire, INSERM U894, Centre de Psychiatrie et Neurosciences, 75014, Paris, France
| | - Rachid C Maroun
- UMR-S U1204, Structure et Activité de Biomolécules Normales et Pathologiques, INSERM/Université d'Evry-Val d'Essonne/Université Paris-Saclay, 91000, Evry, France.
- Laboratoire de Neurobiologie et Pharmacologie Moléculaire, INSERM U894, Centre de Psychiatrie et Neurosciences, 75014, Paris, France.
| |
Collapse
|
41
|
Génin NEJ, Weinzierl ROJ. Nucleotide Loading Modes of Human RNA Polymerase II as Deciphered by Molecular Simulations. Biomolecules 2020; 10:biom10091289. [PMID: 32906795 PMCID: PMC7565877 DOI: 10.3390/biom10091289] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2020] [Revised: 08/31/2020] [Accepted: 09/03/2020] [Indexed: 01/01/2023] Open
Abstract
Mapping the route of nucleoside triphosphate (NTP) entry into the sequestered active site of RNA polymerase (RNAP) has major implications for elucidating the complete nucleotide addition cycle. Constituting a dichotomy that remains to be resolved, two alternatives, direct NTP delivery via the secondary channel (CH2) or selection to downstream sites in the main channel (CH1) prior to catalysis, have been proposed. In this study, accelerated molecular dynamics simulations of freely diffusing NTPs about RNAPII were applied to refine the CH2 model and uncover atomic details on the CH1 model that previously lacked a persuasive structural framework to illustrate its mechanism of action. Diffusion and binding of NTPs to downstream DNA, and the transfer of a preselected NTP to the active site, are simulated for the first time. All-atom simulations further support that CH1 loading is transcription factor IIF (TFIIF) dependent and impacts catalytic isomerization. Altogether, the alternative nucleotide loading systems may allow distinct transcriptional landscapes to be expressed.
Collapse
Affiliation(s)
- Nicolas E. J. Génin
- Institut de Chimie Organique et Analytique, Université d’Orléans, 45100 Orléans, France;
| | | |
Collapse
|
42
|
Torrens-Fontanals M, Stepniewski TM, Aranda-García D, Morales-Pastor A, Medel-Lacruz B, Selent J. How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs. Int J Mol Sci 2020; 21:E5933. [PMID: 32824756 PMCID: PMC7460635 DOI: 10.3390/ijms21165933] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 08/11/2020] [Accepted: 08/15/2020] [Indexed: 01/08/2023] Open
Abstract
G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.
Collapse
Affiliation(s)
- Mariona Torrens-Fontanals
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Tomasz Maciej Stepniewski
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
- InterAx Biotech AG, PARK innovAARE, 5234 Villigen, Switzerland
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, 02-093 Warsaw, Poland
| | - David Aranda-García
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Adrián Morales-Pastor
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Brian Medel-Lacruz
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Jana Selent
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| |
Collapse
|
43
|
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
|
44
|
Bhattarai A, Devkota S, Bhattarai S, Wolfe MS, Miao Y. Mechanisms of γ-Secretase Activation and Substrate Processing. ACS CENTRAL SCIENCE 2020; 6:969-983. [PMID: 32607444 PMCID: PMC7318072 DOI: 10.1021/acscentsci.0c00296] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2020] [Indexed: 05/14/2023]
Abstract
Amyloid β-peptide, the principal component of characteristic cerebral plaques of Alzheimer's disease (AD), is produced through intramembrane proteolysis of the amyloid precursor protein (APP) by γ-secretase. Despite the importance in the pathogenesis of AD, the mechanisms of intramembrane proteolysis and substrate processing by γ-secretase remain poorly understood. Here, complementary all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method and biochemical experiments were combined to investigate substrate processing of wildtype and mutant APP by γ-secretase. The GaMD simulations captured spontaneous activation of γ-secretase, with hydrogen bonded catalytic aspartates and water poised for proteolysis of APP at the ε cleavage site. Furthermore, GaMD simulations revealed that familial AD mutations I45F and T48P enhanced the initial ε cleavage between residues Leu49-Val50, while M51F mutation shifted the ε cleavage site to the amide bond between Thr48-Leu49. Detailed analysis of the GaMD simulations allowed us to identify distinct low-energy conformational states of γ-secretase, different secondary structures of the wildtype and mutant APP substrate, and important active-site subpockets for catalytic function of the enzyme. The simulation findings were highly consistent with experimental analyses of APP proteolytic products using mass spectrometry and Western blotting. Taken together, the GaMD simulations and biochemical experiments have enabled us to elucidate the mechanisms of γ-secretase activation and substrate processing, which should facilitate rational computer-aided drug design targeting this functionally important enzyme.
Collapse
Affiliation(s)
- Apurba Bhattarai
- Center
for Computational Biology and Department of Molecular Biosciences, and Department of
Medicinal Chemistry, School of Pharmacy, University of Kansas, Lawrence, Kansas 66047, United States
| | - Sujan Devkota
- Center
for Computational Biology and Department of Molecular Biosciences, and Department of
Medicinal Chemistry, School of Pharmacy, University of Kansas, Lawrence, Kansas 66047, United States
| | - Sanjay Bhattarai
- Center
for Computational Biology and Department of Molecular Biosciences, and Department of
Medicinal Chemistry, School of Pharmacy, University of Kansas, Lawrence, Kansas 66047, United States
| | - Michael S. Wolfe
- Center
for Computational Biology and Department of Molecular Biosciences, and Department of
Medicinal Chemistry, School of Pharmacy, University of Kansas, Lawrence, Kansas 66047, United States
- (M.S.W.)
| | - Yinglong Miao
- Center
for Computational Biology and Department of Molecular Biosciences, and Department of
Medicinal Chemistry, School of Pharmacy, University of Kansas, Lawrence, Kansas 66047, United States
- (Y.M.)
| |
Collapse
|
45
|
Probing the Molecular Mechanism of Rifampin Resistance Caused by the Point Mutations S456L and D441V on Mycobacterium tuberculosis RNA Polymerase through Gaussian Accelerated Molecular Dynamics Simulation. Antimicrob Agents Chemother 2020; 64:AAC.02476-19. [PMID: 32393493 DOI: 10.1128/aac.02476-19] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2019] [Accepted: 05/06/2020] [Indexed: 01/08/2023] Open
Abstract
Rifampin is the first-line antituberculosis drug, with Mycobacterium tuberculosis RNA polymerase as the molecular target. Unfortunately, M. tuberculosis strains that are resistant to rifampin have been identified in clinical settings, which limits its therapeutic effects. In clinical isolates, S531L and D516V (in Escherichia coli) are two common mutated codons in the gene rpoB, corresponding to S456L and D441V in M. tuberculosis However, the resistance mechanism at the molecular level is still elusive. In this work, Gaussian accelerated molecular dynamics simulations were performed to uncover the resistance mechanism of rifampin due to S456L and D441V mutations at the atomic level. The binding free energy analysis revealed that the reduction in the ability of two mutants to bind rifampin is mainly due to a decrease in electrostatic interaction, specifically, a decrease in the energy contribution of the R454 residue. R454 acts as an anchor and forms stable hydrogen bond interaction with rifampin, allowing rifampin to be stably incorporated in the center of the binding pocket. However, the disappearance of the hydrogen bond between R454 and the mutated residues increases the flexibility of the side chain of R454. The conformation of R454 changes, and the hydrogen bond interaction between it and rifampin is disrupted. As result, the rifampin molecule moves to the outside of the pocket, and the binding affinity decreases. Overall, these findings can provide useful information for understanding the drug resistance mechanism of rifampin and also can give theoretical guidance for further design of novel inhibitors to overcome the drug resistance.
Collapse
|
46
|
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
|
47
|
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
|
48
|
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
|
49
|
Limongelli V. Ligand binding free energy and kinetics calculation in 2020. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1455] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Vittorio Limongelli
- Faculty of Biomedical Sciences, Institute of Computational Science – Center for Computational Medicine in Cardiology Università della Svizzera italiana (USI) Lugano Switzerland
- Department of Pharmacy University of Naples “Federico II” Naples Italy
| |
Collapse
|
50
|
Li X, Dai J, Ni D, He X, Zhang H, Zhang J, Fu Q, Liu Y, Lu S. Insight into the mechanism of allosteric activation of PI3Kα by oncoprotein K-Ras4B. Int J Biol Macromol 2019; 144:643-655. [PMID: 31816384 DOI: 10.1016/j.ijbiomac.2019.12.020] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2019] [Revised: 12/02/2019] [Accepted: 12/03/2019] [Indexed: 12/16/2022]
Abstract
Ras is a key member in the superfamily of small GTPase. Transforming between GTP-bound active state and GDP-bound inactive state in response to exogenous signals, Ras serves as a binary switch in various signaling pathways. One of its downstream effectors is phosphatidylinositol-4,5-bisphosphate 3-kinase α (PI3Kα), which phosphorylates phosphatidylinositol 4,5-bisphosphate into phosphatidylinositol 3,4,5-trisphosphate in the PI3K/Akt/mTOR pathway and mediates an array of important cellular activities including cell growth, migration and survival. Hyperactivation of PI3Kα induced by the Ras isoform K-Ras4B has been unveiled as a key event during the oncogenesis of pancreatic ductal adenocarcinoma, but the underlying mechanism of how K-Ras4B allosterically activates PI3Kα still remains largely unsolved. Here, we employed accelerated molecular dynamic simulations and allosteric pathway analysis to explore into the activation process of PI3Kα by K-Ras4B and unraveled the underlying structural mechanisms. We found that K-Ras4B binding induced more conformational dynamics within PI3Kα and triggered its step-wise transition from a self-inhibited state towards an activated state. Moreover, K-Ras4B binding markedly disrupted the interactions along the p110/p85 interface, especially the ones between nSH2 in p85 and its nearby functional domains in p110 like C2, helical, and kinase domains. The altered inter-domain interactions exposed the kinase domain, which promoted the membrane association and substrate phosphorylation of PI3Kα, thereby facilitating its activation. In particular, the community networks and allosteric pathways analysis further revealed that in PI3Kα/K-Ras4B system, allosteric signaling regulating p110/p85 interaction was rewired from the helical domain to the kinase domain and several important residues and their related allosteric pathways mediating PI3Kα autoinhibition were bypassed. The obtained structural mechanisms provide an in-depth mechanistic insight into the allosteric activation of PI3Kα by K-Ras4B as well as shed light on its drug discovery.
Collapse
Affiliation(s)
- Xinyi Li
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China
| | - Jinyuan Dai
- Chemical Engineering and Technology, School of Chemical Engineering, East China University of Science and Technology, Shanghai 201424, China
| | - Duan Ni
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China
| | - Xinheng He
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China
| | - Hao Zhang
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China
| | - Jian Zhang
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China
| | - Qiang Fu
- Department of Orthopedics, Shanghai General Hospital, Shanghai Jiao Tong University, School of Medicine, Shanghai 200080, China.
| | - Yaqin Liu
- Medicinal Bioinformatics Center, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China.
| | - Shaoyong Lu
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China; Medicinal Bioinformatics Center, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China.
| |
Collapse
|