1
|
Bao H, Wang W, Sun H, Chen J. The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes. J Biomol Struct Dyn 2024; 42:3363-3381. [PMID: 37216340 DOI: 10.1080/07391102.2023.2213355] [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/23/2023] [Accepted: 05/04/2023] [Indexed: 05/24/2023]
Abstract
Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Huayin Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
2
|
Yu Z, Wang Z, Cui X, Cao Z, Zhang W, Sun K, Hu G. Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics. Molecules 2024; 29:645. [PMID: 38338389 PMCID: PMC10856033 DOI: 10.3390/molecules29030645] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2023] [Revised: 01/01/2024] [Accepted: 01/25/2024] [Indexed: 02/12/2024] Open
Abstract
The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.
Collapse
Affiliation(s)
- Zhiping Yu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| | - Zhen Wang
- Pingyin People’s Hospital, Jinan 250400, China; (Z.W.); (X.C.)
| | - Xiuzhen Cui
- Pingyin People’s Hospital, Jinan 250400, China; (Z.W.); (X.C.)
| | - Zanxia Cao
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| | - Wanyunfei Zhang
- School of Science, Xi’an Polytechnic University, Xi’an 710048, China; (W.Z.); (K.S.)
| | - Kunxiao Sun
- School of Science, Xi’an Polytechnic University, Xi’an 710048, China; (W.Z.); (K.S.)
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| |
Collapse
|
3
|
Madhu MK, Shewani K, Murarka RK. Biased Signaling in Mutated Variants of β 2-Adrenergic Receptor: Insights from Molecular Dynamics Simulations. J Chem Inf Model 2024; 64:449-469. [PMID: 38194225 DOI: 10.1021/acs.jcim.3c01481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
The molecular basis of receptor bias in G protein-coupled receptors (GPCRs) caused by mutations that preferentially activate specific intracellular transducers over others remains poorly understood. Two experimentally identified biased variants of β2-adrenergic receptors (β2AR), a prototypical GPCR, are a triple mutant (T68F, Y132A, and Y219A) and a single mutant (Y219A); the former bias the receptor toward the β-arrestin pathway by disfavoring G protein engagement, while the latter induces G protein signaling explicitly due to selection against GPCR kinases (GRKs) that phosphorylate the receptor as a prerequisite of β-arrestin binding. Though rigorous characterizations have revealed functional implications of these mutations, the atomistic origin of the observed transducer selectivity is not clear. In this study, we investigated the allosteric mechanism of receptor bias in β2AR using microseconds of all-atom Gaussian accelerated molecular dynamics (GaMD) simulations. Our observations reveal distinct rearrangements in transmembrane helices, intracellular loop 3, and critical residues R1313.50 and Y3267.53 in the conserved motifs D(E)RY and NPxxY for the mutant receptors, leading to their specific transducer interactions. Moreover, partial dissociation of G protein from the receptor core is observed in the simulations of the triple mutant in contrast to the single mutant and wild-type receptor. The reorganization of allosteric communications from the extracellular agonist BI-167107 to the intracellular receptor-transducer interfaces drives the conformational rearrangements responsible for receptor bias in the single and triple mutants. The molecular insights into receptor bias of β2AR presented here could improve the understanding of biased signaling in GPCRs, potentially opening new avenues for designing novel therapeutics with fewer side-effects and superior efficacy.
Collapse
Affiliation(s)
- Midhun K Madhu
- Department of Biological Sciences, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| | - Kunal Shewani
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| | - Rajesh K Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| |
Collapse
|
4
|
Pawnikar S, Magenheimer BS, Munoz EN, Haldane A, Maser RL, Miao Y. Activation of Polycystin-1 Signaling by Binding of Stalk-derived Peptide Agonists. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.06.574465. [PMID: 38260358 PMCID: PMC10802338 DOI: 10.1101/2024.01.06.574465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
Polycystin-1 (PC1) is the membrane protein product of the PKD1 gene whose mutation is responsible for 85% of the cases of autosomal dominant polycystic kidney disease (ADPKD). ADPKD is primarily characterized by the formation of renal cysts and potential kidney failure. PC1 is an atypical G protein-coupled receptor (GPCR) consisting of 11 transmembrane helices and an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal (NTF) and membrane-embedded C-terminal (CTF) fragments. Recently, signaling activation of the PC1 CTF was shown to be regulated by a stalk tethered agonist (TA), a distinct mechanism observed in the adhesion GPCR family. A novel allosteric activation pathway was elucidated for the PC1 CTF through a combination of Gaussian accelerated molecular dynamics (GaMD), mutagenesis and cellular signaling experiments. Here, we show that synthetic, soluble peptides with 7 to 21 residues derived from the stalk TA, in particular, peptides including the first 9 residues (p9), 17 residues (p17) and 21 residues (p21) exhibited the ability to re-activate signaling by a stalkless PC1 CTF mutant in cellular assays. To reveal molecular mechanisms of stalk peptide-mediated signaling activation, we have applied a novel Peptide GaMD (Pep-GaMD) algorithm to elucidate binding conformations of selected stalk peptide agonists p9, p17 and p21 to the stalkless PC1 CTF. The simulations revealed multiple specific binding regions of the stalk peptide agonists to the PC1 protein including an "intermediate" bound yet inactive state. Our Pep-GaMD simulation findings were consistent with the cellular assay experimental data. Binding of peptide agonists to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of the PC1 CTF signaling activation mechanism. Using sequence covariation analysis of PC1 homologs, we further showed that the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. Therefore, structural dynamic insights into the mechanisms of PC1 activation by stalk-derived peptide agonists have enabled an in-depth understanding of PC1 signaling. They will form a foundation for development of PC1 as a therapeutic target for the treatment of ADPKD.
Collapse
Affiliation(s)
- Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Brenda S. Magenheimer
- Departments of Biochemistry and Molecular Biology, Temple University, Philadelphia, PA 19122
- Clinical Laboratory Sciences, Temple University, Philadelphia, PA 19122
| | | | - Allan Haldane
- Dept of Physics, and Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA 19122
| | - Robin L. Maser
- Departments of Biochemistry and Molecular Biology, Temple University, Philadelphia, PA 19122
- Clinical Laboratory Sciences, Temple University, Philadelphia, PA 19122
- The Jared Grantham Kidney Institute, University of Kansas Medical Center, Kansas City, KS 66160
| | - Yinglong Miao
- Department of Pharmacology and Computational Medicine Program, University of North Carolina – Chapel Hill, Chapel Hill, NC 27599
| |
Collapse
|
5
|
Wu Z, Han Z, Tao L, Sun X, Su J, Hu J, Li C. Dynamic Insights into the Self-Activation Pathway and Allosteric Regulation of the Orphan G-Protein-Coupled Receptor GPR52. J Chem Inf Model 2023; 63:5847-5862. [PMID: 37651308 DOI: 10.1021/acs.jcim.3c00672] [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: 09/02/2023]
Abstract
Within over 800 members of G-protein-coupled receptors, there are numerous orphan receptors whose endogenous ligands are largely unknown, providing many opportunities for novel drug discovery. However, the lack of an in-depth understanding of the intrinsic working mechanism for orphan receptors severely limits the related rational drug design. The G-protein-coupled receptor 52 (GPR52) is a unique orphan receptor that constitutively increases cellular 5'-cyclic adenosine monophosphate (cAMP) levels without binding any exogenous agonists and has been identified as a promising therapeutic target for central nervous system disorders. Although recent structural biology studies have provided snapshots of both active and inactive states of GPR52, the mechanism of the conformational transition between these states remains unclear. Here, an acceptable self-activation pathway for GPR52 was proposed through 6 μs Gaussian accelerated molecular dynamics (GaMD) simulations, in which the receptor spontaneously transitions from the active state to that matching the inactive crystal structure. According to the three intermediate states of the receptor obtained by constructing a reweighted potential of mean force, how the allosteric regulation occurs between the extracellular orthosteric binding pocket and the intracellular G-protein-binding site is revealed. Combined with the independent gradient model, several important microswitch residues and the allosteric communication pathway that directly links the two regions are both identified. Transfer entropy calculations not only reveal the complex allosteric signaling within GPR52 but also confirm the unique role of ECL2 in allosteric regulation, which is mutually validated with the results of GaMD simulations. Overall, this work elucidates the allosteric mechanism of GPR52 at the atomic level, providing the most detailed information to date on the self-activation of the orphan receptor.
Collapse
Affiliation(s)
- Zhixiang Wu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Lianci Tao
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Xiaohan Sun
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Jingjie Su
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Jianping Hu
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu 610106, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| |
Collapse
|
6
|
Liu K, Guo F, Ma Y, Yu X, Fu X, Li W, Han W. Functionalized Fullerene Potentially Inhibits SARS-CoV-2 Infection by Modulating Spike Protein Conformational Changes. Int J Mol Sci 2023; 24:14471. [PMID: 37833919 PMCID: PMC10572755 DOI: 10.3390/ijms241914471] [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: 08/13/2023] [Revised: 09/18/2023] [Accepted: 09/20/2023] [Indexed: 10/15/2023] Open
Abstract
The disease of SARS-CoV-2 has caused considerable morbidity and mortality globally. Spike proteins on the surface of SARS-CoV-2 allow it to bind with human cells, leading to infection. Fullerenes and their derivatives are promising SARS-CoV-2 inhibitors and drug-delivery vehicles. In this study, Gaussian accelerated molecular dynamics simulations and the Markov state model were employed to delve into the inhibitory mechanism of Fullerene-linear-polyglycerol-b-amine sulfate (F-LGPS) on spike proteins. During the study, it was discovered that fullerene derivatives can operate at the interface of the receptor-binding domain (RBD) and the N-terminal domain (NTD), keeping structural domains in a downward conformation. It was also observed that F-LGPS demonstrated superior inhibitory effects on the XBB variant in comparison to the wild-type variant. This study yielded invaluable insights for the potential development of efficient therapeutics targeting the spike protein of SARS-CoV-2.
Collapse
Affiliation(s)
- Kaifeng Liu
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China; (K.L.); (Y.M.)
| | - Fangfang Guo
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (F.G.); (X.Y.); (X.F.)
| | - Yingying Ma
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China; (K.L.); (Y.M.)
| | - Xiangyu Yu
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (F.G.); (X.Y.); (X.F.)
| | - Xueqi Fu
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (F.G.); (X.Y.); (X.F.)
| | - Wannan Li
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (F.G.); (X.Y.); (X.F.)
| | - Weiwei Han
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China; (K.L.); (Y.M.)
| |
Collapse
|
7
|
Poudel B, Rajeshwar T R, Vanegas JM. Membrane mediated mechanical stimuli produces distinct active-like states in the AT1 receptor. Nat Commun 2023; 14:4690. [PMID: 37542033 PMCID: PMC10403497 DOI: 10.1038/s41467-023-40433-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Accepted: 07/27/2023] [Indexed: 08/06/2023] Open
Abstract
The Angiotensin II Type 1 (AT1) receptor is one of the most widely studied GPCRs within the context of biased signaling. While the AT1 receptor is activated by agonists such as the peptide AngII, it can also be activated by mechanical stimuli such as membrane stretch or shear in the absence of a ligand. Despite the importance of mechanical activation of the AT1 receptor in biological processes such as vasoconstriction, little is known about the structural changes induced by external physical stimuli mediated by the surrounding lipid membrane. Here, we present a systematic simulation study that characterizes the activation of the AT1 receptor under various membrane environments and mechanical stimuli. We show that stability of the active state is highly sensitive to membrane thickness and tension. Structural comparison of membrane-mediated vs. agonist-induced activation shows that the AT1 receptor has distinct active conformations. This is supported by multi-microsecond free energy calculations that show unique landscapes for the inactive and various active states. Our modeling results provide structural insights into the mechanical activation of the AT1 receptor and how it may produce different functional outcomes within the framework of biased agonism.
Collapse
Affiliation(s)
- Bharat Poudel
- Materials Science Graduate Program, The University of Vermont, Burlington, VT, 05405, USA
| | - Rajitha Rajeshwar T
- Department of Physics, The University of Vermont, Burlington, VT, 05405, USA
| | - Juan M Vanegas
- Materials Science Graduate Program, The University of Vermont, Burlington, VT, 05405, USA.
- Department of Physics, The University of Vermont, Burlington, VT, 05405, USA.
- Department of Biochemistry and Biophysics, Oregon State University, Corvallis, OR, 97330, USA.
| |
Collapse
|
8
|
Lu X, Shi X, Fan J, Li M, Zhang Y, Lu S, Xu G, Chen Z. Mechanistic Elucidation of Activation/Deactivation Signal Transduction within Neurotensin Receptor 1 Triggered by 'Driver Chemical Groups' of Modulators: A Comparative Molecular Dynamics Simulation. Pharmaceutics 2023; 15:2000. [PMID: 37514186 PMCID: PMC10385606 DOI: 10.3390/pharmaceutics15072000] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 07/16/2023] [Accepted: 07/18/2023] [Indexed: 07/30/2023] Open
Abstract
Small-molecule modulators of neurotensin receptor 1 (NTSR1), a class A G-protein-coupled receptor (GPCR), has emerged as promising therapeutic agent for psychiatric disorders and cancer. Interestingly, a chemical group substitution in NTSR1 modulators can launch different types of downstream regulation, highlighting the significance of deciphering the internal fine-tuning mechanism. Here, we conducted a synergistic application of a Gaussian accelerated molecular dynamics simulation, a conventional molecular dynamics simulation, and Markov state models (MSM) to investigate the underlying mechanism of 'driver chemical groups' of modulators triggering inverse signaling. The results indicated that the flexibility of the leucine moiety in NTSR1 agonists contributes to the inward displacement of TM7 through a loosely coupled allosteric pathway, while the rigidity of the adamantane moiety in NTSR1 antagonists leads to unfavorable downward transduction of agonistic signaling. Furthermore, we found that R3226.54, Y3196.51, F3537.42, R1483.32, S3567.45, and S3577.46 may play a key role in inducing the activation of NTSR1. Together, our findings not only highlight the ingenious signal transduction within class A GPCRs but also lay a foundation for the development of targeted drugs harboring different regulatory functions of NTSR1.
Collapse
Affiliation(s)
- Xun Lu
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
| | - Xinchao Shi
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
| | - Jigang Fan
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
| | - Mingyu Li
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
| | - Yuxiang Zhang
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
| | - Shaoyong Lu
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
| | - Guanghuan Xu
- Department of VIP Clinic, Changhai Hospital, Affiliated to Navy Medical University, Shanghai 200433, China
| | - Ziqiang Chen
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai 200433, China
| |
Collapse
|
9
|
Vuckovic Z, Wang J, Pham V, Mobbs JI, Belousoff MJ, Bhattarai A, Burger WAC, Thompson G, Yeasmin M, Nawaratne V, Leach K, van der Westhuizen ET, Khajehali E, Liang YL, Glukhova A, Wootten D, Lindsley CW, Tobin A, Sexton P, Danev R, Valant C, Miao Y, Christopoulos A, Thal DM. Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics. eLife 2023; 12:83477. [PMID: 37248726 DOI: 10.7554/elife.83477] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Accepted: 04/12/2023] [Indexed: 05/31/2023] Open
Abstract
Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M4 muscarinic acetylcholine receptor (M4 mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M4 mAChR. By determining the cryo-electron microscopy structures of the M4 mAChR, bound to a cognate Gi1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand-receptor-transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.
Collapse
Affiliation(s)
- Ziva Vuckovic
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
| | - Vi Pham
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Jesse I Mobbs
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Matthew J Belousoff
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
| | - Wessel A C Burger
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Geoff Thompson
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Mahmuda Yeasmin
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Vindhya Nawaratne
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Katie Leach
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Emma T van der Westhuizen
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Elham Khajehali
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Yi-Lynn Liang
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Alisa Glukhova
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Denise Wootten
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Craig W Lindsley
- Department of Pharmacology, Warren Center for Neuroscience Drug Discovery and Department of Chemistry, Warren Center for Neuroscience Drug Discovery, Vanderbilt University, Nashville, United States
| | - Andrew Tobin
- The Centre for Translational Pharmacology, Advanced Research Centre (ARC), College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, United Kingdom
| | - Patrick Sexton
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Radostin Danev
- Graduate School of Medicine, University of Tokyo, Tokyo, Japan
| | - Celine Valant
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
| | - Arthur Christopoulos
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- Neuromedicines Discovery Centre, Monash University, Parkville, Australia
| | - David M Thal
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| |
Collapse
|
10
|
Shao Q, Xiong M, Li J, Hu H, Su H, Xu Y. Unraveling the catalytic mechanism of SARS-CoV-2 papain-like protease with allosteric modulation of C270 mutation using multiscale computational approaches. Chem Sci 2023; 14:4681-4696. [PMID: 37181765 PMCID: PMC10171076 DOI: 10.1039/d3sc00166k] [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: 01/11/2023] [Accepted: 04/10/2023] [Indexed: 05/16/2023] Open
Abstract
Papain-like protease (PLpro) is a promising therapeutic target against SARS-CoV-2, but its restricted S1/S2 subsites pose an obstacle in developing active site-directed inhibitors. We have recently identified C270 as a novel covalent allosteric site for SARS-CoV-2 PLpro inhibitors. Here we present a theoretical investigation of the proteolysis reaction catalyzed by the wild-type SARS-CoV-2 PLpro as well as the C270R mutant. Enhanced sampling MD simulations were first performed to explore the influence of C270R mutation on the protease dynamics, and sampled thermodynamically favorable conformations were then submitted to MM/PBSA and QM/MM MD simulations for thorough characterization of the protease-substrate binding and covalent reactions. The disclosed proteolysis mechanism of PLpro, as characterized by the occurrence of proton transfer from the catalytic C111 to H272 prior to the substrate binding and with deacylation being the rate-determining step of the whole proteolysis process, is not completely identical to that of the 3C-like protease, another key cysteine protease of coronaviruses. The C270R mutation alters the structural dynamics of the BL2 loop that indirectly impairs the catalytic function of H272 and reduces the binding of the substrate with the protease, ultimately showing an inhibitory effect on PLpro. Together, these results provide a comprehensive understanding at the atomic level of the key aspects of SARS-CoV-2 PLpro proteolysis, including the catalytic activity allosterically regulated by C270 modification, which is crucial to the follow-up inhibitor design and development.
Collapse
Affiliation(s)
- Qiang Shao
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences Shanghai 201203 China
- University of Chinese Academy of Sciences Beijing 100049 China
| | - Muya Xiong
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences Shanghai 201203 China
- University of Chinese Academy of Sciences Beijing 100049 China
| | - Jiameng Li
- School of Chinese Materia Medica, Nanjing University of Chinese Medicine Nanjing 210023 China
| | - Hangchen Hu
- School of Pharmaceutical Science and Technology, Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences Hangzhou 310024 China
| | - Haixia Su
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences Shanghai 201203 China
| | - Yechun Xu
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences Shanghai 201203 China
- University of Chinese Academy of Sciences Beijing 100049 China
- School of Chinese Materia Medica, Nanjing University of Chinese Medicine Nanjing 210023 China
- School of Pharmaceutical Science and Technology, Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences Hangzhou 310024 China
| |
Collapse
|
11
|
Li H, Urs NM, Horenstein N. Computational insights into ligand-induced G protein and β-arrestin signaling of the dopamine D1 receptor. J Comput Aided Mol Des 2023; 37:227-244. [PMID: 37060492 DOI: 10.1007/s10822-023-00503-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2022] [Accepted: 03/28/2023] [Indexed: 04/16/2023]
Abstract
The dopamine D1 receptor (D1R), is a class A G protein coupled-receptor (GPCR) which has been a promising drug target for psychiatric and neurological disorders such as Parkinson's disease (PD). Previous studies have suggested that therapeutic effects can be realized by targeting the β-arrestin signaling pathway of dopamine receptors, while overactivation of the G protein-dependent pathways leads to side effects, such as dyskinesias. Therefore, it is highly desirable to develop a D1R ligand that selectively regulates the β-arrestin pathway. Currently, most D1R agonists are signaling-balanced and stimulate both G protein and β-arrestin pathways, with a few reports of G protein biased ligands. However, identification and characterization of β-arrestin biased D1R agonists has been a challenge thus far. In this study, we implemented Gaussian accelerated molecular dynamics (GaMD) simulations to provide valuable computational insights into the possible underlying molecular mechanism of the different signaling properties of two catechol and two non-catechol D1R agonists that are either G protein biased or signaling-balanced. Dynamic network analysis further identified critical residues in the allosteric signaling network of D1R for each ligand at different conformational or binding states. Some of these residues are crucial for G protein or arrestin signals of GPCRs based on previous studies. Finally, we provided a molecular design strategy which can be utilized by medicinal chemists to develop potential β-arrestin biased D1R ligands. The proposed hypotheses are experimentally testable and can guide the development of safer and more effective medications for a variety of CNS disorders.
Collapse
Affiliation(s)
- Haoxi Li
- Department of Chemistry, University of Florida, Gainesville, FL, 32611, USA
| | - Nikhil M Urs
- Department of Pharmacology and Therapeutics, University of Florida, Gainesville, FL, 32610, USA
| | - Nicole Horenstein
- Department of Chemistry, University of Florida, Gainesville, FL, 32611, USA.
| |
Collapse
|
12
|
Hasse T, Zhang Z, Huang YMM. Molecular dynamics study reveals key disruptors of MEIG1-PACRG interaction. Proteins 2023; 91:555-566. [PMID: 36444670 PMCID: PMC10374433 DOI: 10.1002/prot.26449] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Revised: 11/07/2022] [Accepted: 11/15/2022] [Indexed: 11/30/2022]
Abstract
Interactions between the meiosis-expressed gene 1 (MEIG1) and Parkin co-regulated gene (PACRG) protein are critical in the formation of mature sperm cells. Targeting either MEIG1 or PACRG protein could be a contraceptive strategy. The W50A and Y68A mutations on MEIG1 are known to interrupt the MEIG1-PACRG interactions resulting in defective sperm cells. However, the details about how the mutants disrupt the protein-protein binding are not clear. In this study, we reveal insights on MEIG1 and PACRG protein dynamics by applying Gaussian-accelerated molecular dynamics (GaMD) simulations and post-GaMD analysis. Our results show that the mutations destabilize the protein-protein interfacial interaction. The effect of the Y68A mutation is more significant than W50A as Y68 forms stronger polar interactions with PACRG. Because both human and mouse models demonstrate similar dynamic properties, the findings from mouse proteins can be applied to the human system. Moreover, we report a potential ligand binding pocket on the MEIG1 and PACRG interaction surface that could be a target for future drug design to inhibit the MEIG1-PACRG interaction. PACRG shows more qualified pockets along the protein-protein interface, implying that it is a better target than MEIG1. Our work provides a fundamental understanding of MEIG1 and PACRG protein dynamics, paving the way for drug discovery in male-based contraception.
Collapse
Affiliation(s)
- Timothy Hasse
- Department of Physics and Astronomy, Wayne State University, Detroit, Michigan, USA
| | - Zhibing Zhang
- Department of Physiology, Wayne State University School of Medicine, Detroit, Michigan, USA.,Department of Obstetrics and Gynecology, Wayne State University, Detroit, Michigan, USA
| | - Yu-Ming M Huang
- Department of Physics and Astronomy, Wayne State University, Detroit, Michigan, USA
| |
Collapse
|
13
|
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
|
14
|
Shi S, Zheng L, Ren Y, Wang Z. Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations. Molecules 2023; 28:molecules28072886. [PMID: 37049650 PMCID: PMC10095679 DOI: 10.3390/molecules28072886] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Revised: 03/19/2023] [Accepted: 03/20/2023] [Indexed: 04/14/2023] Open
Abstract
G12 mutations heavily affect conformational transformation and activity of KRAS. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the GDP-bound wild-type (WT), G12A, G12D, and G12R KRAS to probe mutation-mediated impacts on conformational alterations of KRAS. The results indicate that three G12 mutations obviously affect the structural flexibility and internal dynamics of the switch domains. The analyses of the free energy landscapes (FELs) suggest that three G12 mutations induce more conformational states of KRAS and lead to more disordered switch domains. The principal component analysis shows that three G12 mutations change concerted motions and dynamics behavior of the switch domains. The switch domains mostly overlap with the binding region of KRAS to its effectors. Thus, the high disorder states and concerted motion changes of the switch domains induced by G12 mutations affect the activity of KRAS. The analysis of interaction network of GDP with KRAS signifies that the instability in the interactions of GDP and magnesium ion with the switch domain SW1 drives the high disordered state of the switch domains. This work is expected to provide theoretical aids for understanding the function of KRAS.
Collapse
Affiliation(s)
- Shuhua Shi
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Linqi Zheng
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Yonglian Ren
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Ziyu Wang
- School of Science, Shandong Jianzhu University, Jinan 250101, China
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| |
Collapse
|
15
|
Zhao Y, Zhang J, Zhang H, Gu S, Deng Y, Tu Y, Hou T, Kang Y. Sigmoid Accelerated Molecular Dynamics: An Efficient Enhanced Sampling Method for Biosystems. J Phys Chem Lett 2023; 14:1103-1112. [PMID: 36700836 DOI: 10.1021/acs.jpclett.2c03688] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Gaussian accelerated molecular dynamics (GaMD) is recognized as a popular enhanced sampling method for tackling long-standing challenges in biomolecular simulations. Inspired by GaMD, Sigmoid accelerated molecular dynamics (SaMD) is proposed in this work by adding a Sigmoid boost potential to improve the balance between the highest acceleration and accurate reweighting. Compared with GaMD, SaMD extends the accessible time scale and improves the computational efficiency as tested in three tasks. In the alanine dipeptide task, SaMD can produce the free energy landscape with better accuracy and efficiency. In the chignolin folding task, the estimated Gibbs free energy difference can converge to the experimental value ∼30% faster. In the protein-ligand binding task, the bound conformations are closer to the crystal structure with a minimal ligand root-mean-square deviation of 1.7 Å. The binding of the ligand XK263 to the HIV protease is reproduced by SaMD in ∼60% less simulation time.
Collapse
Affiliation(s)
- Yihao Zhao
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou310058, Zhejiang, China
| | - Jintu Zhang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou310058, Zhejiang, China
- CarbonSilicon AI Technology Company, Ltd., Hangzhou310018, Zhejiang, China
| | - Haotian Zhang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou310058, Zhejiang, China
- CarbonSilicon AI Technology Company, Ltd., Hangzhou310018, Zhejiang, China
| | - Shukai Gu
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou310058, Zhejiang, China
| | - Yafeng Deng
- CarbonSilicon AI Technology Company, Ltd., Hangzhou310018, Zhejiang, China
| | - Yaoquan Tu
- Division of Theoretical Chemistry and Biology, Department of Chemistry, KTH Royal Institute of Technology, 114 28Stockholm, Sweden
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou310058, Zhejiang, China
| | - Yu Kang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou310058, Zhejiang, China
| |
Collapse
|
16
|
Bao HY, Wang W, Sun HB, Chen JZ. Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2023; 34:65-89. [PMID: 36762439 DOI: 10.1080/1062936x.2023.2165542] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 12/31/2022] [Indexed: 06/18/2023]
Abstract
Probing binding modes of GDP, GTP and GNP to NRAS are of significance for understanding the regulation mechanism on the activity of RAS proteins. Four separate Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the apo, GDP-, GTP- and GNP-bound NRAS. Dynamics analyses suggest that binding of three ligands highly affects conformational states of the switch domains from NRAS, which disturbs binding of NRAS to its effectors. The analyses of free energy landscapes (FELs) indicate that binding of GDP, GTP and GNP induces more energetic states of NRAS compared to the apo NRAS but the presence of GNP makes the switch domains more ordered than binding of GDP and GNP. The information of interaction networks of ligands with NRAS reveals that the π-π interaction of residue F28 and the salt bridge interactions of K16 and D119 with ligands stabilize binding of GDP, GTP and GNP to NRAS. Meanwhile magnesium ion plays a bridge role in interactions of ligands with NRAS, which is favourable for associations of GDP, GTP and GNP with NRAS. This work is expected to provide useful information for deeply understanding the function and activity of NRAS.
Collapse
Affiliation(s)
- H Y Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - J Z Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
17
|
Xiaoli A, Yuzhen N, Qiong Y, Yang L, Yao X, Bing Z. Investigating the Dynamic Binding Behavior of PMX53 Cooperating with Allosteric Antagonist NDT9513727 to C5a Anaphylatoxin Chemotactic Receptor 1 through Gaussian Accelerated Molecular Dynamics and Free-Energy Perturbation Simulations. ACS Chem Neurosci 2022; 13:3502-3511. [PMID: 36428153 DOI: 10.1021/acschemneuro.2c00556] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
Abstract
C5a anaphylatoxin chemotactic receptor 1 (C5aR1) is an important target in anti-inflammatory therapeutics. The cyclic peptide antagonist PMX53 binds to the orthosteric site located in the extracellular vestibule of C5aR1, and the non-peptide antagonist NDT9513727 binds to the allosteric site formed by the middle region of TM3 (trans-membrane helix), TM4, and TM5. We catch a sight of the variational binding mode of PMX53 during the Gaussian accelerated molecular dynamic (GaMD) simulations. In the binary complex of C5aR1 and PMX53, the PMX53 takes a dynamic binding mechanism during the simulation. Namely, the side chain of Arg6 of PMX53 extends to TM6-TM7 (pose 1) or swings to TM5 (pose 2), forming a salt bridge with Glu199. Meanwhile, in the ternary complex of C5aR1 with PMX53 and NDT9513727, the side chain of Arg6 of PMX53 swings to TM5 (pose 2) from extending to TM6-TM7 (pose 1) at the beginning of the GaMD simulation. In subsequent simulation, PMX53 stabilizes in the pose 2 binding mode by forming a stable salt bridge with Glu199. The free-energy perturbation (FEP) calculations demonstrate that pose 1 (ΔGbinding = -10.94 kcal/mol) is more stable in the binary complex and pose 2 (ΔGbinding = -7.91 kcal/mol) is unstable because of highly dynamic TM5. NDT9513727 interacts directly with TM4 and TM5 and stabilizes the hydrophobic stack between the extracellular sides of the two helices. Therefore, pose 2 (ΔGbinding = -16.27 kcal/mol) is notably stable than pose 1 (ΔGbinding = -9.78 kcal/mol) in the ternary complex. The identification of a novel binding mode of PMX53 and the detailed structural information of PMX53 interacting with a receptor obtained by GaMD simulations will be helpful in designing potent antagonists of C5aR1.
Collapse
Affiliation(s)
- An Xiaoli
- Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China
| | - Niu Yuzhen
- Shandong Laboratory of Yantai Advanced Materials and Green Manufacturing, Yantai 264006, China.,Yantai Zhongke Research Institute of Advanced Materials and Green Chemical Engineering, Yantai 264006, China
| | - Yang Qiong
- Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China
| | - Lei Yang
- Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China
| | - Xiaojun Yao
- Dr. Neher's Biophysics Laboratory for Innovative Drug Discovery, State Key Laboratory of Quality Research in Chinese Medicine, Macau Institute for Applied Research in Medicine and Health, Macau University of Science and Technology, Taipa, Macau 999078, China
| | - Zhitong Bing
- Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China.,Advanced Energy Science and Technology Guangdong Laboratory, Huizhou 516000, China
| |
Collapse
|
18
|
Chen X, Yuan Y, Chen Y, Yu J, Wang J, Chen J, Guo Y, Pu X. Biased Activation Mechanism Induced by GPCR Heterodimerization: Observations from μOR/δOR Dimers. J Chem Inf Model 2022; 62:5581-5600. [PMID: 36377848 DOI: 10.1021/acs.jcim.2c00962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
GPCRs regulate multiple intracellular signaling cascades. Biasedly activating one signaling pathway over the others provides additional clinical utility to optimize GPCR-based therapies. GPCR heterodimers possess different functions from their monomeric states, including their selectivity to different transducers. However, the biased signaling mechanism induced by the heterodimerization remains unclear. Motivated by the issue, we select an important GPCR heterodimer (μOR/δOR heterodimer) as a case and use microsecond Gaussian accelerated molecular dynamics simulation coupled with potential of mean force and protein structure network (PSN) to probe mechanisms regarding the heterodimerization-induced constitutive β-arrestin activity and efficacy change of the agonist DAMGO. The results show that only the lowest energy state of the μOR/δOR heterodimer, which adopts a slightly outward shift of TM6 and an ICL2 conformation close to the receptor core, can selectively accommodate β-arrestins. PSN further reveals important roles of H8, ICL1, and ICL2 in regulating the constitutive β-arrestin-biased activity for the apo μOR/δOR heterodimer. In addition, the heterodimerization can allosterically alter the binding mode of DAMGO mainly by means of W7.35. Consequently, DAMGO transmits the structural signal mainly through TM6 and TM7 in the dimer, rather than TM3 similar to the μOR monomer, thus changing the efficacy of DAMGO from a balanced agonist to the β-arrestin-biased one. On the other side, the binding of DAMGO to the heterodimer can stabilize μOR/δOR heterodimers through a stronger interaction of TM1/TM1 and H8/H8, accordingly enhancing the interaction of μOR with δOR and the binding affinity of the dimer to the β-arrestin. The agonist DAMGO does not change main compositions of the regulation network from the dimer interface to the transducer binding pocket of the μOR protomer, but induces an increase in the structural communication of the network, which should contribute to the enhanced β-arrestin coupling. Our observations, for the first time, reveal the molecular mechanism of the biased signaling induced by the heterodimerization for GPCRs, which should be beneficial to more comprehensively understand the GPCR bias signaling.
Collapse
Affiliation(s)
- Xin Chen
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Yuan Yuan
- College of Management, Southwest University for Nationalities, Chengdu610041, China
| | - Yichi Chen
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Jin Yu
- Department of Physics and Astronomy, University of California, Irvine, California92697, United States
| | - Jingzhou Wang
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Jianfang Chen
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Yanzhi Guo
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Xuemei Pu
- College of Chemistry, Sichuan University, Chengdu610064, China
| |
Collapse
|
19
|
The full activation mechanism of the adenosine A 1 receptor revealed by GaMD and Su-GaMD simulations. Proc Natl Acad Sci U S A 2022; 119:e2203702119. [PMID: 36215480 PMCID: PMC9586258 DOI: 10.1073/pnas.2203702119] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The full activation process of G protein-coupled receptor (GPCR) plays an important role in cellular signal transduction. However, it remains challenging to simulate the whole process in which the GPCR is recognized and activated by a ligand and then couples to the G protein on a reasonable simulation timescale. Here, we developed a molecular dynamics (MD) approach named supervised (Su) Gaussian accelerated MD (GaMD) by incorporating a tabu-like supervision algorithm into a standard GaMD simulation. By using this Su-GaMD method, from the active and inactive structure of adenosine A1 receptor (A1R), we successfully revealed the full activation mechanism of A1R, including adenosine (Ado)-A1R recognition, preactivation of A1R, and A1R-G protein recognition, in hundreds of nanoseconds of simulations. The binding of Ado to the extracellular side of A1R initiates conformational changes and the preactivation of A1R. In turn, the binding of Gi2 to the intracellular side of A1R causes a decrease in the volume of the extracellular orthosteric site and stabilizes the binding of Ado to A1R. Su-GaMD could be a useful tool to reconstruct or even predict ligand-protein and protein-protein recognition pathways on a short timescale. The intermediate states revealed in this study could provide more detailed complementary structural characterizations to facilitate the drug design of A1R in the future.
Collapse
|
20
|
Yu Z, Su H, Chen J, Hu G. Deciphering Conformational Changes of the GDP-Bound NRAS Induced by Mutations G13D, Q61R, and C118S through Gaussian Accelerated Molecular Dynamic Simulations. Molecules 2022; 27:molecules27175596. [PMID: 36080363 PMCID: PMC9457619 DOI: 10.3390/molecules27175596] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 08/24/2022] [Accepted: 08/25/2022] [Indexed: 11/29/2022] Open
Abstract
The conformational changes in switch domains significantly affect the activity of NRAS. Gaussian-accelerated molecular dynamics (GaMD) simulations of three separate replicas were performed to decipher the effects of G13D, Q16R, and C118S on the conformational transformation of the GDP-bound NRAS. The analyses of root-mean-square fluctuations and dynamics cross-correlation maps indicated that the structural flexibility and motion modes of the switch domains involved in the binding of NRAS to effectors are highly altered by the G13D, Q61R, and C118Smutations. The free energy landscapes (FELs) suggested that mutations induce more energetic states in NRAS than the GDP-bound WT NRAS and lead to high disorder in the switch domains. The FELs also indicated that the different numbers of sodium ions entering the GDP binding regions compensate for the changes in electrostatic environments caused by mutations, especially for G13D. The GDP–residue interactions revealed that the disorder in the switch domains was attributable to the unstable hydrogen bonds between GDP and two residues, V29 and D30. This work is expected to provide information on the energetic basis and dynamics of conformational changes in switch domains that can aid in deeply understanding the target roles of NRAS in anticancer treatment.
Collapse
Affiliation(s)
- Zhiping Yu
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| | - Hongyi Su
- Laoling People's Hospital, Dezhou 253600, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
- Laoling People's Hospital, Dezhou 253600, China
| |
Collapse
|
21
|
Chen J, Wang J, Zeng Q, Wang W, Sun H, Wei B. Exploring the deactivation mechanism of human β2 adrenergic receptor by accelerated molecular dynamic simulations. Front Mol Biosci 2022; 9:972463. [PMID: 36111136 PMCID: PMC9468641 DOI: 10.3389/fmolb.2022.972463] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2022] [Accepted: 08/05/2022] [Indexed: 11/13/2022] Open
Abstract
The β2 adrenergic receptor (β2AR), one of important members of the G protein coupled receptors (GPCRs), has been suggested as an important target for cardiac and asthma drugs. Two replicas of Gaussian accelerated molecular dynamics (GaMD) simulations are performed to explore the deactivation mechanism of the active β2AR bound by three different substrates, including the agonist (P0G), antagonist (JTZ) and inverse agonist (JRZ). The simulation results indicate that the Gs protein is needed to stabilize the active state of the β2AR. Without the Gs protein, the receptor could transit from the active state toward the inactive state. During the transition process, helix TM6 moves toward TM3 and TM5 in geometric space and TM5 shrinks upwards. The intermediate state is captured during the transition process of the active β2AR toward the inactive one, moreover the changes in hydrophobic interaction networks between helixes TM3, TM5, and TM6 and the formation of a salt bridge between residues Arg3.50 and Glu6.30 drive the transition process. We expect that this finding can provide energetic basis and molecular mechanism for further understanding the function and target roles of the β2AR.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
- *Correspondence: Jianzhong Chen, ; Benzheng Wei,
| | - Jian Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Qingkai Zeng
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Benzheng Wei
- Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao, China
- *Correspondence: Jianzhong Chen, ; Benzheng Wei,
| |
Collapse
|
22
|
Copeland M, Do HN, Votapka L, Joshi K, Wang J, Amaro RE, Miao Y. Gaussian Accelerated Molecular Dynamics in OpenMM. J Phys Chem B 2022; 126:5810-5820. [PMID: 35895977 PMCID: PMC9773147 DOI: 10.1021/acs.jpcb.2c03765] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Gaussian accelerated molecular dynamics (GaMD) is a computational technique that provides both unconstrained enhanced sampling and free energy calculations of biomolecules. Here, we present the implementation of GaMD in the OpenMM simulation package and validate it on model systems of alanine dipeptide and RNA folding. For alanine dipeptide, 30 ns GaMD production simulations reproduced free energy profiles of 1000 ns conventional molecular dynamics (cMD) simulations. In addition, GaMD simulations captured the folding pathways of three hyperstable RNA tetraloops (UUCG, GCAA, and CUUG) and binding of the rbt203 ligand to the HIV-1 Tar RNA, both of which involved critical electrostatic interactions such as hydrogen bonding and base stacking. Together with previous implementations, GaMD in OpenMM will allow for wider applications in simulations of proteins, RNA, and other biomolecules.
Collapse
Affiliation(s)
- Matthew Copeland
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Hung N. Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Lane Votapka
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, CA 92093
| | - Keya Joshi
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, CA 92093
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047,To whom correspondence should be addressed:
| |
Collapse
|
23
|
Su M, Paknejad N, Zhu L, Wang J, Do HN, Miao Y, Liu W, Hite RK, Huang XY. Structures of β 1-adrenergic receptor in complex with Gs and ligands of different efficacies. Nat Commun 2022; 13:4095. [PMID: 35835792 PMCID: PMC9283524 DOI: 10.1038/s41467-022-31823-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 07/01/2022] [Indexed: 11/09/2022] Open
Abstract
G-protein-coupled receptors (GPCRs) receive signals from ligands with different efficacies, and transduce to heterotrimeric G-proteins to generate different degrees of physiological responses. Previous studies revealed how ligands with different efficacies activate GPCRs. Here, we investigate how a GPCR activates G-proteins upon binding ligands with different efficacies. We report the cryo-EM structures of β1-adrenergic receptor (β1-AR) in complex with Gs (GαsGβ1Gγ2) and a partial agonist or a very weak partial agonist, and compare them to the β1-AR-Gs structure in complex with a full agonist. Analyses reveal similar overall complex architecture, with local conformational differences. Cellular functional studies with mutations of β1-AR residues show effects on the cellular signaling from β1-AR to the cAMP response initiated by the three different ligands, with residue-specific functional differences. Biochemical investigations uncover that the intermediate state complex comprising β1-AR and nucleotide-free Gs is more stable when binding a full agonist than a partial agonist. Molecular dynamics simulations support the local conformational flexibilities and different stabilities among the three complexes. These data provide insights into the ligand efficacy in the activation of GPCRs and G-proteins.
Collapse
Affiliation(s)
- Minfei Su
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, NY, 10065, USA
| | - Navid Paknejad
- Structural Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
| | - Lan Zhu
- Cancer Center and Department of Pharmacology and Toxicology, Medical College of Wisconsin, Milwaukee, WI, 53226, USA
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Hung Nguyen Do
- 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
| | - Wei Liu
- Cancer Center and Department of Pharmacology and Toxicology, Medical College of Wisconsin, Milwaukee, WI, 53226, USA
| | - Richard K Hite
- Structural Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA.
| | - Xin-Yun Huang
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, NY, 10065, USA.
| |
Collapse
|
24
|
Shahoei R, Pangeni S, Sanders MA, Zhang H, Mladenovic-Lucas L, Roush WR, Halvorsen G, Kelly CV, Granneman JG, Huang YMM. Molecular Modeling of ABHD5 Structure and Ligand Recognition. Front Mol Biosci 2022; 9:935375. [PMID: 35836935 PMCID: PMC9274090 DOI: 10.3389/fmolb.2022.935375] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Accepted: 06/10/2022] [Indexed: 02/06/2023] Open
Abstract
Alpha/beta hydrolase domain-containing 5 (ABHD5), also termed CGI-58, is the key upstream activator of adipose triglyceride lipase (ATGL), which plays an essential role in lipid metabolism and energy storage. Mutations in ABHD5 disrupt lipolysis and are known to cause the Chanarin-Dorfman syndrome. Despite its importance, the structure of ABHD5 remains unknown. In this work, we combine computational and experimental methods to build a 3D structure of ABHD5. Multiple comparative and machine learning-based homology modeling methods are used to obtain possible models of ABHD5. The results from Gaussian accelerated molecular dynamics and experimental data of the apo models and their mutants are used to select the most likely model. Moreover, ensemble docking is performed on representative conformations of ABHD5 to reveal the binding mechanism of ABHD5 and a series of synthetic ligands. Our study suggests that the ABHD5 models created by deep learning-based methods are the best candidate structures for the ABHD5 protein. The mutations of E41, R116, and G328 disturb the hydrogen bonding network with nearby residues and suppress membrane targeting or ATGL activation. The simulations also reveal that the hydrophobic interactions are responsible for binding sulfonyl piperazine ligands to ABHD5. Our work provides fundamental insight into the structure of ABHD5 and its ligand-binding mode, which can be further applied to develop ABHD5 as a therapeutic target for metabolic disease and cancer.
Collapse
Affiliation(s)
- Rezvan Shahoei
- Department of Physics and Astronomy, Wayne State University, Detroit, MI, United States
| | - Susheel Pangeni
- Department of Physics and Astronomy, Wayne State University, Detroit, MI, United States
| | - Matthew A. Sanders
- Center for Molecular Medicine and Genetics, School of Medicine, Wayne State University, Detroit, MI, United States
| | - Huamei Zhang
- Center for Molecular Medicine and Genetics, School of Medicine, Wayne State University, Detroit, MI, United States
| | - Ljiljana Mladenovic-Lucas
- Center for Molecular Medicine and Genetics, School of Medicine, Wayne State University, Detroit, MI, United States
| | - William R. Roush
- Department of Chemistry, Scripps Florida, Jupiter, FL, United States
| | - Geoff Halvorsen
- Department of Chemistry, Scripps Florida, Jupiter, FL, United States
| | - Christopher V. Kelly
- Department of Physics and Astronomy, Wayne State University, Detroit, MI, United States
| | - James G. Granneman
- Center for Molecular Medicine and Genetics, School of Medicine, Wayne State University, Detroit, MI, United States,Center for Integrative Metabolic and Endocrine Research, School of Medicine, Wayne State University, Detroit, MI, United States
| | - Yu-ming M. Huang
- Department of Physics and Astronomy, Wayne State University, Detroit, MI, United States,*Correspondence: Yu-ming M. Huang,
| |
Collapse
|
25
|
Li M, Bao Y, Xu R, La H, Guo J. Critical Extracellular Ca 2+ Dependence of the Binding between PTH1R and a G-Protein Peptide Revealed by MD Simulations. ACS Chem Neurosci 2022; 13:1666-1674. [PMID: 35543321 DOI: 10.1021/acschemneuro.2c00176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
The parathyroid hormone type 1 receptor (PTH1R), a canonical class B GPCR, is regulated by a positive allosteric modulator, extracellular Ca2+. Calcium ions prolong the residence time of PTH on the PTH1R, leading to increased receptor activation and duration of cAMP signaling. But the essential mechanism of the allosteric behavior of PTH1R is not fully understood. Here, extensive molecular dynamics (MD) simulations are performed for the PTH1R-G-protein combinations with and without Ca2+ to describe how calcium ions allosterically engage receptor-G-protein coupling. We find that the binding of Ca2+ stabilizes the conformation of the PTH1R-PTH-spep (the α5 helix of Gs protein) complex, especially the extracellular loop 1 (ECL1). Moreover, the MM-GBSA result indicates that Ca2+ allosterically promotes the interaction between PTH1R and spep, consistent with the observation of steered molecular dynamics (SMD) simulations. We further illuminate the possible allosteric signaling pathway from the stable Ca2+-coupling site to the intracellular G-protein binding site. These results unveil structural determinants for Ca2+ allosterism in the PTH1R-PTH-spep complex and give insights into pluridimensional GPCR signaling regulated by calcium ions.
Collapse
Affiliation(s)
- Mengrong Li
- College of Life Sciences, Nanjing Agricultural University, Nanjing 210095, China
| | - Yiqiong Bao
- College of Life Sciences, Nanjing Agricultural University, Nanjing 210095, China
| | - Ran Xu
- College of Life Sciences, Nanjing Agricultural University, Nanjing 210095, China
| | - Honggui La
- College of Life Sciences, Nanjing Agricultural University, Nanjing 210095, China
| | - Jingjing Guo
- College of Life Sciences, Nanjing Agricultural University, Nanjing 210095, China
- Engineering Research Centre of Applied Technology on Machine Translation and Artificial Intelligence, Faculty of Applied Science, Macao Polytechnic University, Macao 999078, China
| |
Collapse
|
26
|
Pawnikar S, Miao Y. Mechanism of Peptide Agonist Binding in CXCR4 Chemokine Receptor. Front Mol Biosci 2022; 9:821055. [PMID: 35359589 PMCID: PMC8963245 DOI: 10.3389/fmolb.2022.821055] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 02/14/2022] [Indexed: 01/07/2023] Open
Abstract
Chemokine receptors are key G-protein-coupled receptors (GPCRs) that control cell migration in immune system responses, development of cardiovascular and central nervous systems, and numerous diseases. In particular, the CXCR4 chemokine receptor promotes metastasis, tumor growth and angiogenesis in cancers. CXCR4 is also used as one of the two co-receptors for T-tropic HIV-1 entry into host cells. Therefore, CXCR4 serves as an important therapeutic target for treating cancers and HIV infection. Apart from the CXCL12 endogenous peptide agonist, previous studies suggested that the first 17 amino acids of CXCL12 are sufficient to activate CXCR4. Two 17-residue peptides with positions 1-4 mutated to RSVM and ASLW functioned as super and partial agonists of CXCR4, respectively. However, the mechanism of peptide agonist binding in CXCR4 remains unclear. Here, we have investigated this mechanism through all-atom simulations using a novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method. The Pep-GaMD simulations have allowed us to explore representative binding conformations of each peptide and identify critical low-energy states of CXCR4 activated by the super versus partial peptide agonists. Our simulations have provided important mechanistic insights into peptide agonist binding in CXCR4, which are expected to facilitate rational design of new peptide modulators of CXCR4 and other chemokine receptors.
Collapse
|
27
|
Wang J, Miao Y. Protein-Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD): Characterization of Protein Binding Thermodynamics and Kinetics. J Chem Theory Comput 2022; 18:1275-1285. [PMID: 35099970 DOI: 10.1021/acs.jctc.1c00974] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Protein-protein interactions (PPIs) play key roles in many fundamental biological processes such as cellular signaling and immune responses. However, it has proven challenging to simulate repetitive protein association and dissociation in order to calculate binding free energies and kinetics of PPIs due to long biological timescales and complex protein dynamics. To address this challenge, we have developed a new computational approach to all-atom simulations of PPIs based on a robust Gaussian accelerated molecular dynamics (GaMD) technique. The method, termed "PPI-GaMD", selectively boosts interaction potential energy between protein partners to facilitate their slow dissociation. Meanwhile, another boost potential is applied to the remaining potential energy of the entire system to effectively model the protein's flexibility and rebinding. PPI-GaMD has been demonstrated on a model system of the ribonuclease barnase interactions with its inhibitor barstar. Six independent 2 μs PPI-GaMD simulations have captured repetitive barstar dissociation and rebinding events, which enable calculations of the protein binding thermodynamics and kinetics simultaneously. The calculated binding free energies and kinetic rate constants agree well with the experimental data. Furthermore, PPI-GaMD simulations have provided mechanistic insights into barstar binding to barnase, which involves long-range electrostatic interactions and multiple binding pathways, being consistent with previous experimental and computational findings of this model system. In summary, PPI-GaMD provides a highly efficient and easy-to-use approach for binding free energy and kinetics calculations of PPIs.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
28
|
Do HN, Wang J, Bhattarai A, Miao Y. GLOW: A Workflow Integrating Gaussian-Accelerated Molecular Dynamics and Deep Learning for Free Energy Profiling. J Chem Theory Comput 2022; 18:1423-1436. [PMID: 35200019 PMCID: PMC9773012 DOI: 10.1021/acs.jctc.1c01055] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We introduce a Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and free energy profiling workflow (GLOW) to predict molecular determinants and map free energy landscapes of biomolecules. All-atom GaMD-enhanced sampling simulations are first performed on biomolecules of interest. Structural contact maps are then calculated from GaMD simulation frames and transformed into images for building DL models using a convolutional neural network. Important structural contacts are further determined from DL models of attention maps of the structural contact gradients, which allow us to identify the system reaction coordinates. Finally, free energy profiles are calculated for the selected reaction coordinates through energetic reweighting of the GaMD simulations. We have also successfully demonstrated GLOW for the characterization of activation and allosteric modulation of a G protein-coupled receptor, using the adenosine A1 receptor (A1AR) as a model system. GLOW findings are highly consistent with previous experimental and computational studies of the A1AR, while also providing further mechanistic insights into the receptor function. In summary, GLOW provides a systematic approach to mapping free energy landscapes of biomolecules. The GLOW workflow and its user manual can be downloaded at http://miaolab.org/GLOW.
Collapse
Affiliation(s)
- Hung N. Do
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Jinan Wang
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Apurba Bhattarai
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Yinglong Miao
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047,Corresponding author:
| |
Collapse
|
29
|
Madhu MK, Debroy A, Murarka RK. Molecular Insights into Phosphorylation-Induced Allosteric Conformational Changes in a β 2-Adrenergic Receptor. J Phys Chem B 2022; 126:1917-1932. [PMID: 35196859 DOI: 10.1021/acs.jpcb.1c08610] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The large conformational flexibility of G protein-coupled receptors (GPCRs) has been a puzzle in structural and pharmacological studies for the past few decades. Apart from structural rearrangements induced by ligands, enzymatic phosphorylations by GPCR kinases (GRKs) at the carboxy-terminal tail (C-tail) of a GPCR also make conformational alterations to the transmembrane helices and facilitates the binding of one of its transducer proteins named β-arrestin. The phosphorylation-induced conformational transition of the receptor that causes specific binding to β-arrestin but prevents the association of other transducers such as G proteins lacks atomistic understanding and is elusive to experimental studies. Using microseconds of all-atom conventional and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigate the allosteric mechanism of phosphorylation induced-conformational changes in β2-adrenergic receptor, a well-characterized GPCR model system. Free energy profiles reveal that the phosphorylated receptor samples a new conformational state in addition to the canonical active state corroborating with recent nuclear magnetic resonance experimental findings. The new state has a smaller intracellular cavity that is likely to accommodate β-arrestin better than G protein. Using contact map and inter-residue interaction energy calculations, we found the phosphorylated C-tail adheres to the cytosolic surface of the transmembrane domain of the receptor. Transfer entropy calculations show that the C-tail residues drive the correlated motions of TM residues, and the allosteric signal is relayed via several residues at the cytosolic surface. Our results also illustrate how the redistribution of inter-residue nonbonding interaction couples with the allosteric communication from the phosphorylated C-tail to the transmembrane. Atomistic insight into phosphorylation-induced β-arrestin specific conformation is therapeutically important to design drugs with higher efficacy and fewer side effects. Our results, therefore, open novel opportunities to fine-tune β-arrestin bias in GPCR signaling.
Collapse
Affiliation(s)
- Midhun K Madhu
- Department of Biological Sciences, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| | - Annesha Debroy
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| | - Rajesh K Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| |
Collapse
|
30
|
Mansoor S, Kayık G, Durdagi S, Sensoy O. Mechanistic insight into the impact of a bivalent ligand on the structure and dynamics of a GPCR oligomer. Comput Struct Biotechnol J 2022; 20:925-936. [PMID: 35242285 PMCID: PMC8861583 DOI: 10.1016/j.csbj.2022.01.016] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2021] [Revised: 12/25/2021] [Accepted: 01/17/2022] [Indexed: 12/25/2022] Open
Abstract
Bivalent ligand; modulates conformational preferences of and correlations among microswitches, strengthens interaction between G protein and the receptor and also water channel formation.
Development of effective bivalent ligands has become the focus of intensive research toward modulation of G protein-coupled receptor (GPCR) oligomers, particularly in the field of GPCR pharmacology. Experimental studies have shown that they increased binding affinity and signaling potency compared to their monovalent counterparts, yet underlying molecular mechanism remains elusive. To address this, we performed accelerated molecular dynamics simulations on bivalent-ligand bound Adenosine 2A receptor (A2AR) dimer in the context of a modeled tetramer, which consists of A2AR and dopamine 2 receptor (D2R) homodimers and their cognate G proteins. Our results demonstrate that bivalent ligand impacted interactions between pharmacophore groups and ligand binding residues, thus modulating allosteric communication network and water channel formed within the receptor. Moreover, it also strengthens contacts between receptor and G protein, by modulating the volume of ligand binding pocket and intracellular domain of the receptor. Importantly, we showed that impact evoked by the bivalent ligand on A2AR dimer was also transmitted to apo D2R, which is part of the neighboring D2R dimer. To the best of our knowledge, this is the first study that provides a mechanistic insight into the impact of a bivalent ligand on dynamics of a GPCR oligomer. Consequently, this will pave the way for development of effective ligands for modulation of GPCR oligomers and hence treatment of crucial diseases such as Parkinson’s disease and cancer.
Collapse
Affiliation(s)
- Samman Mansoor
- School of Engineering and Natural Sciences, Department of Biomedical Engineering and Bioinformatics, Istanbul Medipol University, Istanbul 34810, Turkey
| | - Gülru Kayık
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University, Istanbul, Turkey
| | - Serdar Durdagi
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University, Istanbul, Turkey
| | - Ozge Sensoy
- Regenerative and Restorative Medicine Research Center (REMER), Research Institute for Health Sciencesand Technologies (SABITA), Istanbul Medipol University, 34810 Istanbul, Turkey
- School of Engineering and Natural Sciences, Department of Computer Engineering, Istanbul Medipol University, Turkey
- Corresponding author at: Regenerative and Restorative Medicine Research Center (REMER), Research Institute for Health Sciencesand Technologies (SABITA), Istanbul Medipol University, 34810 Istanbul, Turkey.
| |
Collapse
|
31
|
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
|
32
|
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
|
33
|
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
|
34
|
Higashida R, Matsunaga Y. Enhanced Conformational Sampling of Nanobody CDR H3 Loop by Generalized Replica-Exchange with Solute Tempering. Life (Basel) 2021; 11:life11121428. [PMID: 34947959 PMCID: PMC8706460 DOI: 10.3390/life11121428] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Revised: 12/13/2021] [Accepted: 12/14/2021] [Indexed: 12/22/2022] Open
Abstract
The variable domains of heavy-chain antibodies, known as nanobodies, are potential substitutes for IgG antibodies. They have similar affinities to antigens as antibodies, but are more heat resistant. Their small size allows us to exploit computational approaches for structural modeling or design. Here, we investigate the applicability of an enhanced sampling method, a generalized replica-exchange with solute tempering (gREST) for sampling CDR-H3 loop structures of nanobodies. In the conventional replica-exchange methods, temperatures of only a whole system or scaling parameters of a solute molecule are selected for temperature or parameter exchange. In gREST, we can flexibly select a part of a solute molecule and a part of the potential energy terms as a parameter exchange region. We selected the CDR-H3 loop and investigated which potential energy term should be selected for the efficient sampling of the loop structures. We found that the gREST with dihedral terms can explore a global conformational space, but the relaxation to the global equilibrium is slow. On the other hand, gREST with all the potential energy terms can sample the equilibrium distribution, but the structural exploration is slower than with dihedral terms. The lessons learned from this study can be applied to future studies of loop modeling.
Collapse
|
35
|
Wang J, Lan L, Wu X, Xu L, Miao Y. Mechanism of RNA recognition by a Musashi RNA-binding protein. Curr Res Struct Biol 2021; 4:10-20. [PMID: 34988468 PMCID: PMC8695263 DOI: 10.1016/j.crstbi.2021.12.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Revised: 10/31/2021] [Accepted: 12/07/2021] [Indexed: 12/21/2022] Open
Abstract
The Musashi RNA-binding proteins (RBPs) regulate translation of target mRNAs and maintenance of cell stemness and tumorigenesis. Musashi-1 (MSI1), long considered as an intestinal and neural stem cell marker, has been more recently found to be over expressed in many cancers. It has served as an important drug target for treating acute myeloid leukemia and solid tumors such as ovarian, colorectal and bladder cancer. One of the reported binding targets of MSI1 is Numb, a negative regulator of the Notch signaling. However, the dynamic mechanism of Numb RNA binding to MSI1 remains unknown, largely hindering effective drug design targeting this critical interaction. Here, we have performed extensive all-atom microsecond-timescale simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method, which successfully captured multiple times of spontaneous and highly accurate binding of the Numb RNA from bulk solvent to the MSI1 protein target site. GaMD simulations revealed that Numb RNA binding to MSI1 involved largely induced fit in both the RNA and protein. The simulations also identified important low-energy intermediate conformational states during RNA binding, in which Numb interacted mainly with the β2-β3 loop and C terminus of MSI1. The mechanistic understanding of RNA binding obtained from our GaMD simulations is expected to facilitate rational structure-based drug design targeting MSI1 and other RBPs.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology, University of Kansas, Lawrence, KS, 66047, USA
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Lan Lan
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Xiaoqing Wu
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Liang Xu
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
- Department of Radiation Oncology, The University of Kansas Cancer Center, Kansas City, KS, 66160, USA
| | - Yinglong Miao
- Center for Computational Biology, University of Kansas, Lawrence, KS, 66047, USA
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| |
Collapse
|
36
|
Chen J, Liu J, Yuan Y, Chen X, Zhang F, Pu X. Molecular Mechanisms of Diverse Activation Stimulated by Different Biased Agonists for the β2-Adrenergic Receptor. J Chem Inf Model 2021; 62:5175-5192. [PMID: 34802238 DOI: 10.1021/acs.jcim.1c01016] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
β2AR is an important drug target protein involving many diseases. Biased drugs induce specific signaling and provide additional clinical utility to optimize β2AR-based therapies. However, the biased signaling mechanism has not been elucidated. Motivated by the issue, we chose four agonists with divergent bias (balanced agonist, G-protein-biased agonist, and β-arrestin-biased agonists) and utilized Gaussian accelerated molecular dynamics simulation coupled with a dynamic network to probe the molecular mechanisms of distinct biased activation induced by the structural differences between the four agonists. Our simulations reveal that the G-protein-biased agonist induces an open conformation with the outward shifts of TM6 and TM7 for the intracellular domain, which will be beneficial to couple G protein. In contrast, the β-arrestin-biased agonists regulate an occluded conformation with a slightly outward movement of TM6 and an inward shift of TM7, which should favor β-arrestin signaling. The balanced agonist does not induce an observable outward shift for TM6 but, along with a slight tilt for TM7, leads to an inactive-like conformation. In addition, our results reveal the first time that ICL3 presents specific conformations with different agonists. The G-protein-biased agonist drives ICL3 to open so that the G protein-binding pocket can be available, while the β-arrestin-biased agonists induce ICL3 to form a closed conformation with a stable local α-helix. MM/PBSA analysis further reveals that the hydroxyl groups in the resorcinol of the G-protein-biased agonist form strong interactions with Y5.38 and S5.42, thus preventing tilting of the TM5 extracellular end. The catechol of the balanced agonist and the β-arrestin-biased ones induces the rearrangement of two hydrophobic residues F6.52 and W6.48. However, different from the balanced agonist, the ethyl substituent of β-arrestin-biased agonists forms additional hydrophobic interactions with W6.48 and F6.51 after the rearrangement, which should contribute to the β-arrestin bias. The shortest pathway analysis further reveals that the three residues Y7.43, N7.45, and N7.49 are crucial for allosterically regulating G-protein-biased signaling, while the two residues W6.48 and F6.44 make an important contribution to regulate β-arrestin-biased signaling. For the balanced agonist NE, the allosteric regulation pathway simultaneously involves the residue associated with G-protein-biased signaling like S5.46 and the residues related to β-arrestin-biased signaling like W6.48 and F6.44, thus producing unbiased signaling. The observations could advance our understanding of the biased activation mechanism on class A GPCRs and provide a useful guideline for the design of biased drugs.
Collapse
Affiliation(s)
- Jianfang Chen
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Jiangting Liu
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Yuan Yuan
- College of Management, Southwest University for Nationalities, Chengdu 610041, China
| | - Xin Chen
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Fuhui Zhang
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Xuemei Pu
- College of Chemistry, Sichuan University, Chengdu 610064, China
| |
Collapse
|
37
|
Höring C, Conrad M, Söldner CA, Wang J, Sticht H, Strasser A, Miao Y. Specific Engineered G Protein Coupling to Histamine Receptors Revealed from Cellular Assay Experiments and Accelerated Molecular Dynamics Simulations. Int J Mol Sci 2021; 22:10047. [PMID: 34576210 PMCID: PMC8467750 DOI: 10.3390/ijms221810047] [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: 08/12/2021] [Revised: 09/15/2021] [Accepted: 09/15/2021] [Indexed: 01/29/2023] Open
Abstract
G protein-coupled receptors (GPCRs) are targets of extracellular stimuli and hence occupy a key position in drug discovery. By specific and not yet fully elucidated coupling profiles with α subunits of distinct G protein families, they regulate cellular responses. The histamine H2 and H4 receptors (H2R and H4R) are prominent members of Gs- and Gi-coupled GPCRs. Nevertheless, promiscuous G protein and selective Gi signaling have been reported for the H2R and H4R, respectively, the molecular mechanism of which remained unclear. Using a combination of cellular experimental assays and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigated the coupling profiles of the H2R and H4R to engineered mini-G proteins (mG). We obtained coupling profiles of the mGs, mGsi, or mGsq proteins to the H2R and H4R from the mini-G protein recruitment assays using HEK293T cells. Compared to H2R-mGs expressing cells, histamine responses were weaker (pEC50, Emax) for H2R-mGsi and -mGsq. By contrast, the H4R selectively bound to mGsi. Similarly, in all-atom GaMD simulations, we observed a preferential binding of H2R to mGs and H4R to mGsi revealed by the structural flexibility and free energy landscapes of the complexes. Although the mG α5 helices were consistently located within the HR binding cavity, alternative binding orientations were detected in the complexes. Due to the specific residue interactions, all mG α5 helices of the H2R complexes adopted the Gs-like orientation toward the receptor transmembrane (TM) 6 domain, whereas in H4R complexes, only mGsi was in the Gi-like orientation toward TM2, which was in agreement with Gs- and Gi-coupled GPCRs structures resolved by X-ray/cryo-EM. These cellular and molecular insights support (patho)physiological profiles of the histamine receptors, especially the hitherto little studied H2R function in the brain, as well as of the pharmacological potential of H4R selective drugs.
Collapse
Affiliation(s)
- Carina Höring
- Institute of Pharmacy, Faculty of Chemistry and Pharmacy, University of Regensburg, 93040 Regensburg, Germany
| | - Marcus Conrad
- Bioinformatik, Institut für Biochemie, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Fahrstraße 17, 91054 Erlangen, Germany
| | - Christian A Söldner
- Bioinformatik, Institut für Biochemie, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Fahrstraße 17, 91054 Erlangen, Germany
| | - Jinan Wang
- Department of Computational Biology and Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| | - Heinrich Sticht
- Bioinformatik, Institut für Biochemie, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Fahrstraße 17, 91054 Erlangen, Germany
- Erlangen National High Performance Computing Center (NHR@FAU), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), 91058 Erlangen, Germany
| | - Andrea Strasser
- Institute of Pharmacy, Faculty of Chemistry and Pharmacy, University of Regensburg, 93040 Regensburg, Germany
| | - Yinglong Miao
- Department of Computational Biology and Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA
| |
Collapse
|
38
|
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: 84] [Impact Index Per Article: 28.0] [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
|
39
|
Positive allosteric mechanisms of adenosine A 1 receptor-mediated analgesia. Nature 2021; 597:571-576. [PMID: 34497422 PMCID: PMC8711093 DOI: 10.1038/s41586-021-03897-2] [Citation(s) in RCA: 63] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Accepted: 08/11/2021] [Indexed: 02/08/2023]
Abstract
The adenosine A1 receptor (A1R) is a promising therapeutic target for non-opioid analgesic agents to treat neuropathic pain1,2. However, development of analgesic orthosteric A1R agonists has failed because of a lack of sufficient on-target selectivity as well as off-tissue adverse effects3. Here we show that [2-amino-4-(3,5-bis(trifluoromethyl)phenyl)thiophen-3-yl)(4-chlorophenyl)methanone] (MIPS521), a positive allosteric modulator of the A1R, exhibits analgesic efficacy in rats in vivo through modulation of the increased levels of endogenous adenosine that occur in the spinal cord of rats with neuropathic pain. We also report the structure of the A1R co-bound to adenosine, MIPS521 and a Gi2 heterotrimer, revealing an extrahelical lipid-detergent-facing allosteric binding pocket that involves transmembrane helixes 1, 6 and 7. Molecular dynamics simulations and ligand kinetic binding experiments support a mechanism whereby MIPS521 stabilizes the adenosine-receptor-G protein complex. This study provides proof of concept for structure-based allosteric drug design of non-opioid analgesic agents that are specific to disease contexts.
Collapse
|
40
|
An X, Bai Q, Bing Z, Liu H, Yao X. Insights into the molecular mechanism of positive cooperativity between partial agonist MK-8666 and full allosteric agonist AP8 of hGPR40 by Gaussian accelerated molecular dynamics (GaMD) simulations. Comput Struct Biotechnol J 2021; 19:3978-3989. [PMID: 34377364 PMCID: PMC8313488 DOI: 10.1016/j.csbj.2021.07.008] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Revised: 07/07/2021] [Accepted: 07/08/2021] [Indexed: 10/29/2022] Open
Abstract
Activation of human free fatty acid receptor 1 (FFAR1, also called hGPR40) enhances insulin secretion in a glucose-dependent manner. Hence, the development of selective agonist targeting hGPR40 has been proposed as a therapeutic strategy of type 2 diabetes mellitus. Some agonists targeting hGPR40 were reported. The radioligand-binding studies and the crystal structures reveal that there are multiple sites on GPR40, and there exists positive binding cooperativity between the partial agonist MK-8666 and full allosteric agonist (AgoPAM) AP8. In this work, we carried out long-time Gaussian accelerated molecular dynamics (GaMD) simulations on hGPR40 to shed light on the mechanism of the cooperativity between the two agonists at different sites. Our results reveal that the induced-fit conformational coupling is bidirectional between the two sites. The movements and rotations of TM3, TM4, TM5 and TM6 due to their inherent flexibility are crucial in coupling the conformational changes of the two agonists binding sites. These helices adopt similar conformational states upon alternative ligand or both ligands binding. The Leu1384.57, Leu1865.42 and Leu1905.46 play roles in coordinating the rearrangements of residues in the two pockets, which makes the movements of residues in the two sites like gear movements. These results provide detailed information at the atomic level about the conformational coupling between different sites of GPR40, and also provide the structural information for further design of new agonists of GPR40. In addition, these results suggest that it is necessary by considering the effect of other site bound in structure-based ligands discovery.
Collapse
Affiliation(s)
- Xiaoli An
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry, Lanzhou University, Lanzhou 730000, China
| | - Qifeng Bai
- School of Basic Medical Science, Lanzhou University, Lanzhou, China
| | - Zhitong Bing
- Institute of Modern Physics of Chinese Academy of Sciences, Gansu Province, Lanzhou, China
| | - Huanxiang Liu
- School of Pharmacy, Lanzhou University, Lanzhou 730000, China
| | - Xiaojun Yao
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry, Lanzhou University, Lanzhou 730000, China.,State Key Laboratory of Quality Research in Chinese Medicine, Macau Institute for Applied Research in Medicine and Health, Macau University of Science and Technology, Taipa, Macau, China
| |
Collapse
|
41
|
Structural basis for chemokine recognition and receptor activation of chemokine receptor CCR5. Nat Commun 2021; 12:4151. [PMID: 34230484 PMCID: PMC8260604 DOI: 10.1038/s41467-021-24438-5] [Citation(s) in RCA: 44] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 06/18/2021] [Indexed: 02/06/2023] Open
Abstract
The chemokine receptor CCR5 plays a vital role in immune surveillance and inflammation. However, molecular details that govern its endogenous chemokine recognition and receptor activation remain elusive. Here we report three cryo-electron microscopy structures of Gi1 protein-coupled CCR5 in a ligand-free state and in complex with the chemokine MIP-1α or RANTES, as well as the crystal structure of MIP-1α-bound CCR5. These structures reveal distinct binding modes of the two chemokines and a specific accommodate pattern of the chemokine for the distal N terminus of CCR5. Together with functional data, the structures demonstrate that chemokine-induced rearrangement of toggle switch and plasticity of the receptor extracellular region are critical for receptor activation, while a conserved tryptophan residue in helix II acts as a trigger of receptor constitutive activation.
Collapse
|
42
|
Identification of pyrogallol as a warhead in design of covalent inhibitors for the SARS-CoV-2 3CL protease. Nat Commun 2021; 12:3623. [PMID: 34131140 PMCID: PMC8206078 DOI: 10.1038/s41467-021-23751-3] [Citation(s) in RCA: 95] [Impact Index Per Article: 31.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2021] [Accepted: 05/14/2021] [Indexed: 12/11/2022] Open
Abstract
The ongoing pandemic of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) urgently needs an effective cure. 3CL protease (3CLpro) is a highly conserved cysteine proteinase that is indispensable for coronavirus replication, providing an attractive target for developing broad-spectrum antiviral drugs. Here we describe the discovery of myricetin, a flavonoid found in many food sources, as a non-peptidomimetic and covalent inhibitor of the SARS-CoV-2 3CLpro. Crystal structures of the protease bound with myricetin and its derivatives unexpectedly revealed that the pyrogallol group worked as an electrophile to covalently modify the catalytic cysteine. Kinetic and selectivity characterization together with theoretical calculations comprehensively illustrated the covalent binding mechanism of myricetin with the protease and demonstrated that the pyrogallol can serve as an electrophile warhead. Structure-based optimization of myricetin led to the discovery of derivatives with good antiviral activity and the potential of oral administration. These results provide detailed mechanistic insights into the covalent mode of action by pyrogallol-containing natural products and a template for design of non-peptidomimetic covalent inhibitors against 3CLpros, highlighting the potential of pyrogallol as an alternative warhead in design of targeted covalent ligands. SARS-CoV-2 3CL protease (3CLpro) is essential for coronavirus replication and of great interest as an antiviral drug target. Here, the authors show that the naturally occurring flavonoid myricetin is a non-peptidomimetic and covalent inhibitor of 3CLpro, and they solve crystal structures of 3CLpro with myricetin and derivatives, which reveal that the pyrogallol group covalently modifies the catalytic cysteine.
Collapse
|
43
|
Glaser M, Bruce NJ, Han SB, Wade RC. Simulation of the Positive Inotropic Peptide S100A1ct in Aqueous Environment by Gaussian Accelerated Molecular Dynamics. J Phys Chem B 2021; 125:4654-4666. [PMID: 33944558 DOI: 10.1021/acs.jpcb.1c00902] [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
The S100A1ct peptide, consisting of the C-terminal 20 residues of the S100A1 protein fused to an N-terminal 6-residue hydrophilic tag, has been found to exert a positive inotropic effect, resulting in improved contractile performance of failing cardiac and skeletal muscle without arrhythmic side-effects. The S100A1ct peptide thus has high potential for the treatment of acute heart failure. As a step toward understanding its molecular mechanism of action, and to provide a basis for peptidomimetic design to optimize its properties, we here describe de novo structure predictions and molecular dynamics simulations to characterize the conformational landscape of S100A1ct in aqueous environment. In S100A1, the C-terminal 20 residues form an α-helix, but de novo peptide structure predictions indicate that other conformations are also possible. Conventional molecular dynamics simulations in implicit and explicit solvent corroborated this finding. To ensure adequate sampling, we performed simulations of a tagged 10-residue segment of S100A1ct, and we carried out Gaussian accelerated molecular dynamics simulations of the peptides. These simulations showed that although the helical conformation of S100A1ct was the most energetically stable, the peptide can adopt a range of kinked conformations, suggesting that its activity may be related to its ability to act as a conformational switch.
Collapse
Affiliation(s)
- Manuel Glaser
- Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Informatics for Life, Heidelberg, Germany
| | - Neil J Bruce
- Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Informatics for Life, Heidelberg, Germany
| | - Sungho Bosco Han
- Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, 69120 Heidelberg, Germany
| | - Rebecca C Wade
- Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Informatics for Life, Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, 69120 Heidelberg, Germany.,Zentrum für Molekulare Biologie der Universität Heidelberg (ZMBH), DKFZ-ZMBH Alliance and Interdisciplinary Center for Scientific Computing (IWR), 69120 Heidelberg, Germany
| |
Collapse
|
44
|
Nierzwicki Ł, Arantes PR, Saha A, Palermo G. Establishing the allosteric mechanism in CRISPR-Cas9. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2021; 11:e1503. [PMID: 34322166 PMCID: PMC8315640 DOI: 10.1002/wcms.1503] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Accepted: 08/26/2020] [Indexed: 12/18/2022]
Abstract
Allostery is a fundamental property of proteins, which regulates biochemical information transfer between spatially distant sites. Here, we report on the critical role of molecular dynamics (MD) simulations in discovering the mechanism of allosteric communication within CRISPR-Cas9, a leading genome editing machinery with enormous promises for medicine and biotechnology. MD revealed how allostery intervenes during at least three steps of the CRISPR-Cas9 function: affecting DNA recognition, mediating the cleavage and interfering with the off-target activity. An allosteric communication that activates concerted DNA cleavages was found to led through the L1/L2 loops, which connect the HNH and RuvC catalytic domains. The identification of these "allosteric transducers" inspired the development of novel variants of the Cas9 protein with improved specificity, opening a new avenue for controlling the CRISPR-Cas9 activity. Discussed studies also highlight the critical role of the recognition lobe in the conformational activation of the catalytic HNH domain. Specifically, the REC3 region was found to modulate the dynamics of HNH by sensing the formation of the RNA:DNA hybrid. The role of REC3 was revealed to be particularly relevant in the presence of DNA mismatches. Indeed, interference of REC3 with the RNA:DNA hybrid containing mismatched pairs at specific positions resulted in locking HNH in an inactive "conformational checkpoint" conformation, thereby hampering off-target cleavages. Overall, MD simulations established the fundamental mechanisms underlying the allosterism of CRISPR-Cas9, aiding engineering strategies to develop new CRISPR-Cas9 variants for improved genome editing.
Collapse
Affiliation(s)
- Łukasz Nierzwicki
- Department of Bioengineering, University of California Riverside, Riverside, California
| | - Pablo Ricardo Arantes
- Department of Bioengineering, University of California Riverside, Riverside, California
| | - Aakash Saha
- Department of Bioengineering, University of California Riverside, Riverside, California
| | - Giulia Palermo
- Department of Bioengineering and Department of Chemistry, University of California Riverside, Riverside, California
| |
Collapse
|
45
|
Do HN, Akhter S, Miao Y. Pathways and Mechanism of Caffeine Binding to Human Adenosine A 2A Receptor. Front Mol Biosci 2021; 8:673170. [PMID: 33987207 PMCID: PMC8111288 DOI: 10.3389/fmolb.2021.673170] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Accepted: 03/24/2021] [Indexed: 11/13/2022] Open
Abstract
Caffeine (CFF) is a common antagonist to the four subtypes of adenosine G-protein-coupled receptors (GPCRs), which are critical drug targets for treating heart failure, cancer, and neurological diseases. However, the pathways and mechanism of CFF binding to the target receptors remain unclear. In this study, we have performed all-atom-enhanced sampling simulations using a robust Gaussian-accelerated molecular dynamics (GaMD) method to elucidate the binding mechanism of CFF to human adenosine A2A receptor (A2AAR). Multiple 500–1,000 ns GaMD simulations captured both binding and dissociation of CFF in the A2AAR. The GaMD-predicted binding poses of CFF were highly consistent with the x-ray crystal conformations with a characteristic hydrogen bond formed between CFF and residue N6.55 in the receptor. In addition, a low-energy intermediate binding conformation was revealed for CFF at the receptor extracellular mouth between ECL2 and TM1. While the ligand-binding pathways of the A2AAR were found similar to those of other class A GPCRs identified from previous studies, the ECL2 with high sequence divergence serves as an attractive target site for designing allosteric modulators as selective drugs of the A2AAR.
Collapse
Affiliation(s)
- Hung N Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| | - Sana Akhter
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| |
Collapse
|
46
|
Ma N, Nivedha AK, Vaidehi N. Allosteric communication regulates ligand-specific GPCR activity. FEBS J 2021; 288:2502-2512. [PMID: 33738925 PMCID: PMC9805801 DOI: 10.1111/febs.15826] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Revised: 03/14/2021] [Accepted: 03/17/2021] [Indexed: 01/11/2023]
Abstract
G protein-coupled receptors (GPCRs) are membrane-bound proteins that are ubiquitously expressed in many cell types and take part in mediating multiple signaling pathways. GPCRs are dynamic proteins and exist in an equilibrium between an ensemble of conformational states such as inactive and fully active states. This dynamic nature of GPCRs is one of the factors that confers their basal activity even in the absence of any ligand-mediated activation. Ligands selectively bind and stabilize a subset of the conformations from the ensemble leading to a shift in the equilibrium toward the inactive or the active state depending on the nature of the ligand. This ligand-selective effect is achieved through allosteric communication between the ligand binding site and G protein or β-arrestin coupling site. Similarly, the G protein coupling to the receptor exerts the allosteric effect on the ligand binding region leading to increased binding affinity for agonists and decreased affinity for antagonists or inverse agonists. In this review, we enumerate the current state of our understanding of the mechanism of allosteric communication in GPCRs with a specific focus on the critical role of computational methods in delineating the residues involved in allosteric communication. Analyzing allosteric communication mechanism using molecular dynamics simulations has revealed (a) a structurally conserved mechanism of allosteric communication that regulates the G protein coupling, (b) a rational structure-based approach to designing selective ligands, and (c) an approach to designing allosteric GPCR mutants that are either ligand and G protein or β-arrestin selective.
Collapse
Affiliation(s)
- Ning Ma
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA 91010
| | - Anita K. Nivedha
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA 91010
| | - Nagarajan Vaidehi
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA 91010,to whom correspondence should be addressed:
| |
Collapse
|
47
|
Manna M, Murarka RK. Polyunsaturated Fatty Acid Modulates Membrane-Bound Monomeric α-Synuclein by Modulating Membrane Microenvironment through Preferential Interactions. ACS Chem Neurosci 2021; 12:675-688. [PMID: 33538574 DOI: 10.1021/acschemneuro.0c00694] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
There is ample evidence that both native functions and pathogenic aggregation of α-synuclein are intimately dependent on lipid interactions and fatty acid type; the regulatory mechanism however remains unclear. In the present work, using extensive atomistic molecular dynamics simulations and enhanced-sampling, we have focused on exploring the mechanism of fatty acid dependent regulation of monomeric α-Syn100 in a native synaptic vesicle-like membrane. Our results show that α-Syn100 spontaneously binds to the membrane through its N-terminal region (residues 1-34), where the depth of membrane insertion, the structure, and orientation of the membrane-bound α-Syn100 and its impact on membrane structure are modulated by docosahexaenoic acid (DHA). DHA is a polyunsaturated fatty acid abundantly found in the brain and known to promote the oligomerization of α-synuclein. We found that DHA exhibits marked propensity to interact with monomeric α-Syn100 and modulates the microenvironment of the protein by preferentially sorting DHA-containing phospholipids, depleting other phospholipids and cholesterol as well as increasing the proportion of anionic to neutral lipids in the immediate vicinity of the protein. Owing to the unique conformational flexibility, DHA chains form more lipid-packing defects in the membrane and efficiently coat the membrane-embedded surface of the protein, compared to the saturated and monounsaturated fatty acids. DHA thus makes the bilayer more amiable to protein adsorption and less prone to α-synuclein-induced perturbation associated with cytotoxicity. Indeed, in the absence of DHA, we observed significant thinning of the local bilayer membrane induced by α-Syn100. Though α-Syn100 is predominantly α-helical in membranes studied here, in the presence of DHA we observe formation of β-sheet/β-strands in the C-terminal region (residues 35-100) of α-Syn100, which is extended out from the membrane surface. Notably, DHA induces β structure in the NAC domain of α-Syn100 and promotes extended conformations as well as large solvent exposure of this hydrophobic domain, properties that are known to facilitate self-assembly of α-synuclein. To the best of our knowledge, this study for the first time provides the atomistic insights into DHA-induced regulatory mechanism of monomeric α-synuclein, having implications in protein structure and its physiological/pathological functions.
Collapse
Affiliation(s)
- Moutusi Manna
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462 066, India
- Applied Phycology and Biotechnology Division, CSIR−Central Salt & Marine Chemicals Research Institute (CSIR−CSMCRI), Gijubhai Badheka Marg, Bhavnagar, Gujarat 364002, India
| | - Rajesh K. Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462 066, India
| |
Collapse
|
48
|
Jonniya NA, Sk MF, Kar P. Characterizing an allosteric inhibitor-induced inactive state in with-no-lysine kinase 1 using Gaussian accelerated molecular dynamics simulations. Phys Chem Chem Phys 2021; 23:7343-7358. [DOI: 10.1039/d0cp05733a] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The binding of an allosteric inhibitor in WNK1 leads to the inactive state.
Collapse
Affiliation(s)
- Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, MP
- India
| | - Md Fulbabu Sk
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, MP
- India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, MP
- India
| |
Collapse
|
49
|
Sevy AM, Chen MT, Castor M, Sylvia T, Krishnamurthy H, Ishchenko A, Hsieh CM. Structure- and sequence-based design of synthetic single-domain antibody libraries. Protein Eng Des Sel 2020; 33:6042250. [DOI: 10.1093/protein/gzaa028] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2020] [Revised: 10/07/2020] [Accepted: 10/30/2020] [Indexed: 12/30/2022] Open
Abstract
Abstract
Single-domain antibody fragments known as VHH have emerged in the pharmaceutical industry as useful biotherapeutics. These molecules, which are naturally produced by camelids, share the characteristics of high affinity and specificity with traditional human immunoglobulins, while consisting of only a single heavy chain. Currently, the most common method for generating VHH is via animal immunization, which can be costly and time-consuming. Here we describe the development of a synthetic VHH library for in vitro selection of single domain binders. We combine structure-based design and next-generation sequencing analysis to build a library with characteristics that closely mimic the natural repertoire. To validate the performance of our synthetic library, we isolated VHH against three model antigens (soluble mouse PD-1 ectodomain, amyloid-β peptide, and MrgX1 GPCR) of different sizes and characteristics. We were able to isolate diverse binders targeting different epitopes with high affinity (as high as 5 nM) against all three targets. We then show that anti-mPD-1 binders have functional activity in a receptor blocking assay.
Collapse
Affiliation(s)
| | - Ming-Tang Chen
- Discovery Biologics, Merck & Co., Inc., Boston, MA 02115, USA
| | - Michelle Castor
- Discovery Biologics, Merck & Co., Inc., Boston, MA 02115, USA
| | - Tyler Sylvia
- Discovery Biologics, Merck & Co., Inc., Boston, MA 02115, USA
| | - Harini Krishnamurthy
- Computational and Structural Chemistry, Merck & Co., Inc., West Point, PA 19486, USA
| | - Andrii Ishchenko
- Computational and Structural Chemistry, Merck & Co., Inc., West Point, PA 19486, USA
| | | |
Collapse
|
50
|
Mondal D, Kolev V, Warshel A. Exploring the activation pathway and G i-coupling specificity of the μ-opioid receptor. Proc Natl Acad Sci U S A 2020; 117:26218-26225. [PMID: 33020275 PMCID: PMC7585030 DOI: 10.1073/pnas.2013364117] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Understanding the activation mechanism of the μ-opioid receptor (μ-OR) and its selective coupling to the inhibitory G protein (Gi) is vital for pharmaceutical research aimed at finding treatments for the opioid overdose crisis. Many attempts have been made to understand the mechanism of the μ-OR activation, following the elucidation of new crystal structures such as the antagonist- and agonist-bound μ-OR. However, the focus has not been placed on the underlying energetics and specificity of the activation process. An energy-based picture would not only help to explain this coupling but also help to explore why other possible options are not common. For example, one would like to understand why μ-OR is more selective to Gi than a stimulatory G protein (Gs). Our study used homology modeling and a coarse-grained model to generate all of the possible "end states" of the thermodynamic cycle of the activation of μ-OR. The end points were further used to generate reasonable intermediate structures of the receptor and the Gi to calculate two-dimensional free energy landscapes. The results of the landscape calculations helped to propose a plausible sequence of conformational changes in the μ-OR and Gi system and for exploring the path that leads to its activation. Furthermore, in silico alanine scanning calculations of the last 21 residues of the C terminals of Gi and Gs were performed to shed light on the selective binding of Gi to μ-OR. Overall, the present work appears to demonstrate the potential of multiscale modeling in exploring the action of G protein-coupled receptors.
Collapse
Affiliation(s)
- Dibyendu Mondal
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089
| | - Vesselin Kolev
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089
| |
Collapse
|