1
|
Okwei E, Smith ST, Bender BJ, Allison B, Ganguly S, Geanes A, Zhang X, Ledwitch K, Meiler J. Rosetta's Predictive Ability for Low-Affinity Ligand Binding in Fragment-Based Drug Discovery. Biochemistry 2023; 62:700-709. [PMID: 36626571 DOI: 10.1021/acs.biochem.2c00649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da which weakly bind to the protein of interest. This technique is challenging for computational docking methods as binding is determined by only a few specific interactions. Inaccuracies in the energy function or slight deviations in the docking pose can lead to the prediction of incorrect binding or difficulties in ranking fragments in in silico screening. Here, we test RosettaLigand by docking a series of fragments to a cysteine-depleted variant of the TIM-barrel protein, HisF (UniProtKB Q9X0C6). We compare the computational results with experimental NMR spectroscopy screens. NMR spectroscopy gives details on binding affinities of individual ligands, which allows assessment of the ligand-ranking ability using RosettaLigand and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigand's ability to identify plausible binding poses. From a library screen of 3456 fragments, we identified a set of 31 ligands with intrinsic affinities to HisF with dissociation constants as low as 400 μM. The same library of fragments was blindly screened in silico. RosettaLigand was able to rank binders before non-binders with an area under the curve of the receiver operating characteristics of 0.74. The docking poses observed for binders agreed with the binding pocket identified by NMR chemical shift perturbations for all fragments. Taken together, these results provide a baseline performance of RosettaLigand in a fragment-based drug discovery setting.
Collapse
Affiliation(s)
- Elleansar Okwei
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee37235, United States.,Center for Structural Biology, Vanderbilt University, Nashville, Tennessee37240, United States
| | - Shannon T Smith
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee37240, United States.,Program in Chemical and Physical Biology, Vanderbilt University, Nashville, Tennessee37240, United States
| | - Brian J Bender
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee37240, United States.,Department of Pharmacology, Vanderbilt University, Nashville, Tennessee37240, United States
| | - Brittany Allison
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee37235, United States.,Center for Structural Biology, Vanderbilt University, Nashville, Tennessee37240, United States
| | - Soumya Ganguly
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee37235, United States
| | - Alexander Geanes
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee37235, United States
| | - Xuan Zhang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee37235, United States
| | - Kaitlyn Ledwitch
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee37235, United States.,Center for Structural Biology, Vanderbilt University, Nashville, Tennessee37240, United States
| | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee37235, United States.,Center for Structural Biology, Vanderbilt University, Nashville, Tennessee37240, United States.,Program in Chemical and Physical Biology, Vanderbilt University, Nashville, Tennessee37240, United States.,Department of Pharmacology, Vanderbilt University, Nashville, Tennessee37240, United States.,Institute of Drug Discovery, Faculty of Medicine, University of Leipzig, 04103Leipzig, Germany
| |
Collapse
|
2
|
Bhunia SS, Saxena AK. Efficiency of Homology Modeling Assisted Molecular Docking in G-protein Coupled Receptors. Curr Top Med Chem 2021; 21:269-294. [PMID: 32901584 DOI: 10.2174/1568026620666200908165250] [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] [Received: 06/15/2020] [Revised: 08/20/2020] [Accepted: 09/01/2020] [Indexed: 11/22/2022]
Abstract
BACKGROUND Molecular docking is in regular practice to assess ligand affinity on a target protein crystal structure. In the absence of protein crystal structure, the homology modeling or comparative modeling is the best alternative to elucidate the relationship details between a ligand and protein at the molecular level. The development of accurate homology modeling (HM) and its integration with molecular docking (MD) is essential for successful, rational drug discovery. OBJECTIVE The G-protein coupled receptors (GPCRs) are attractive therapeutic targets due to their immense role in human pharmacology. The GPCRs are membrane-bound proteins with the complex constitution, and the understanding of their activation and inactivation mechanisms is quite challenging. Over the past decade, there has been a rapid expansion in the number of solved G-protein-coupled receptor (GPCR) crystal structures; however, the majority of the GPCR structures remain unsolved. In this context, HM guided MD has been widely used for structure-based drug design (SBDD) of GPCRs. METHODS The focus of this review is on the recent (i) developments on HM supported GPCR drug discovery in the absence of GPCR crystal structures and (ii) application of HM in understanding the ligand interactions at the binding site, virtual screening, determining receptor subtype selectivity and receptor behaviour in comparison with GPCR crystal structures. RESULTS The HM in GPCRs has been extremely challenging due to the scarcity in template structures. In such a scenario, it is difficult to get accurate HM that can facilitate understanding of the ligand-receptor interactions. This problem has been alleviated to some extent by developing refined HM based on incorporating active /inactive ligand information and inducing protein flexibility. In some cases, HM proteins were found to outscore crystal structures. CONCLUSION The developments in HM have been highly operative to gain insights about the ligand interaction at the binding site and receptor functioning at the molecular level. Thus, HM guided molecular docking may be useful for rational drug discovery for the GPCRs mediated diseases.
Collapse
Affiliation(s)
- Shome S Bhunia
- Global Institute of Pharmaceutical Education and Research, Kashipur, Uttarakhand, India
| | - Anil K Saxena
- Division of Medicinal and Process Chemistry, CSIR-CDRI, Lucknow 226031, India
| |
Collapse
|
3
|
Bender BJ, Marlow B, Meiler J. Improving homology modeling from low-sequence identity templates in Rosetta: A case study in GPCRs. PLoS Comput Biol 2020; 16:e1007597. [PMID: 33112852 PMCID: PMC7652349 DOI: 10.1371/journal.pcbi.1007597] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2019] [Revised: 11/09/2020] [Accepted: 08/11/2020] [Indexed: 01/16/2023] Open
Abstract
As sequencing methodologies continue to advance, the availability of protein sequences far outpaces the ability of structure determination. Homology modeling is used to bridge this gap but relies on high-identity templates for accurate model building. G-protein coupled receptors (GPCRs) represent a significant target class for pharmaceutical therapies in which homology modeling could fill the knowledge gap for structure-based drug design. To date, only about 17% of druggable GPCRs have had their structures characterized at atomic resolution. However, modeling of the remaining 83% is hindered by the low sequence identity between receptors. Here we test key inputs in the model building process using GPCRs as a focus to improve the pipeline in two critical ways: Firstly, we use a blended sequence- and structure-based alignment that accounts for structure conservation in loop regions. Secondly, by merging multiple template structures into one comparative model, the best possible template for every region of a target can be used expanding the conformational space sampled in a meaningful way. This optimization allows for accurate modeling of receptors using templates as low as 20% sequence identity, which accounts for nearly the entire druggable space of GPCRs. A model database of all non-odorant GPCRs is made available at www.rosettagpcr.org. Additionally, all protocols are made available with insights into modifications that may improve accuracy at new targets.
Collapse
Affiliation(s)
- Brian Joseph Bender
- Department of Pharmacology, Department of Chemistry, and Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
| | - Brennica Marlow
- Department of Pharmacology, Department of Chemistry, and Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
| | - Jens Meiler
- Department of Pharmacology, Department of Chemistry, and Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
- Institute for Drug Discovery, Leipzig University Medical School, Leipzig, SAC, Germany
| |
Collapse
|
4
|
Sisco NJ, Helsell CVM, Van Horn WD. Competitive Interactions between PIRT, the Cold Sensing Ion Channel TRPM8, and PIP 2 Suggest a Mechanism for Regulation. Sci Rep 2019; 9:14128. [PMID: 31575973 PMCID: PMC6773951 DOI: 10.1038/s41598-019-49912-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2018] [Accepted: 09/02/2019] [Indexed: 01/18/2023] Open
Abstract
TRPM8 is a member of the transient receptor potential ion channel family where it functions as a cold and pain sensor in humans and other higher organisms. Previous studies show that TRPM8 requires the signaling phosphoinositide lipid PIP2 to function. TRPM8 function is further regulated by other diverse mechanisms, including the small modulatory membrane protein PIRT (phosphoinositide regulator of TRP). Like TRPM8, PIRT also binds PIP2 and behavioral studies have shown that PIRT is required for normal TRPM8-mediated cold-sensing. To better understand the molecular mechanism of PIRT regulation of TRPM8, solution nuclear magnetic resonance (NMR) spectroscopy was used to assign the backbone resonances of full-length human PIRT and investigate the direct binding of PIRT to PIP2 and the human TRPM8 S1-S4 transmembrane domain. Microscale thermophoresis (MST) binding studies validate the NMR results and identify a competitive PIRT interaction between PIP2 and the TRPM8 S1-S4 domain. Computational PIP2 docking to a human TRPM8 comparative model was performed to help localize where PIRT may bind TRPM8. Taken together, our data suggest a mechanism where TRPM8, PIRT, and PIP2 form a regulatory complex and PIRT modulation of TRPM8 arises, at least in part, by regulating local concentrations of PIP2 accessible to TRPM8.
Collapse
Affiliation(s)
- Nicholas J Sisco
- School of Molecular Sciences, Arizona State University, Tempe, AZ, 85287, USA
- The Biodesign Institute, Arizona State University, Tempe, AZ, 85281, USA
- The Virginia G. Piper Center for Personalized Diagnostics, Arizona State University, Tempe, AZ, 85281, USA
- The Magnetic Resonance Research Center, Arizona State University, Tempe, AZ, 85287, USA
| | - Cole V M Helsell
- School of Molecular Sciences, Arizona State University, Tempe, AZ, 85287, USA
- The Biodesign Institute, Arizona State University, Tempe, AZ, 85281, USA
- The Virginia G. Piper Center for Personalized Diagnostics, Arizona State University, Tempe, AZ, 85281, USA
- The Magnetic Resonance Research Center, Arizona State University, Tempe, AZ, 85287, USA
| | - Wade D Van Horn
- School of Molecular Sciences, Arizona State University, Tempe, AZ, 85287, USA.
- The Biodesign Institute, Arizona State University, Tempe, AZ, 85281, USA.
- The Virginia G. Piper Center for Personalized Diagnostics, Arizona State University, Tempe, AZ, 85281, USA.
- The Magnetic Resonance Research Center, Arizona State University, Tempe, AZ, 85287, USA.
| |
Collapse
|
5
|
Razali NN, Hashim NH, Leow ATC, Salleh AB. Conformational Design and Characterisation of a Truncated Diamine Oxidase from Arthrobacter globiformis. High Throughput 2018; 7:ht7030021. [PMID: 30149644 PMCID: PMC6163217 DOI: 10.3390/ht7030021] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2018] [Revised: 08/22/2018] [Accepted: 08/23/2018] [Indexed: 11/16/2022] Open
Abstract
A functional mini protein can be developed by miniaturising its size. The minimisation technique provides an excellent model system for studying native enzymes, especially in creating an alternative novel biocatalyst. Miniaturised proteins may have enhanced stability, a crucial characteristic for large-scale production and industrial applications. In this study, a huge enzyme molecule, known as diamine oxidase (DAO, comprising 700 amino acids), was selected to undergo the process. By retaining the arrangement of the original functional sites of DAO in the fourth domain, a mini DAO can be designed via homology modelling. After several downsizing processes, a final configuration of 220 amino acids displayed high binding affinity towards histamine, a short-chain substrate that was catalysed by the parental DAO. The configuration also showed enhanced affinity towards a long-chain substrate known as spermidine. The gene for the designed protein was cloned and expressed in pET102/TOPO vector and overexpressed in E. coli BL21 (DE3). The new mini DAO had similar temperature tolerance and versatile substrates specificity characteristics as its parental protein. An active mini-protein with these characteristics is potentially useful for several applications such as detecting biogenic amines in the biological fluids and the environment that may give rise to health issues.
Collapse
Affiliation(s)
- Nur Nadia Razali
- Laboratory of Molecular Biomedicine, Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Malaysia.
- Enzyme and Microbial Technology Research Center, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia, Serdang 43400, Malaysia.
| | - Nur Hafizah Hashim
- Enzyme and Microbial Technology Research Center, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia, Serdang 43400, Malaysia.
- Department of Cell and Molecular Biology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Malaysia.
| | - Adam Thean Chor Leow
- Laboratory of Molecular Biomedicine, Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Malaysia.
- Enzyme and Microbial Technology Research Center, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia, Serdang 43400, Malaysia.
- Department of Cell and Molecular Biology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Malaysia.
| | - Abu Bakar Salleh
- Laboratory of Molecular Biomedicine, Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Malaysia.
- Enzyme and Microbial Technology Research Center, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia, Serdang 43400, Malaysia.
| |
Collapse
|
6
|
Fu DY, Meiler J. Predictive Power of Different Types of Experimental Restraints in Small Molecule Docking: A Review. J Chem Inf Model 2018; 58:225-233. [PMID: 29286651 DOI: 10.1021/acs.jcim.7b00418] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Incorporating experimental restraints is a powerful method of increasing accuracy in computational protein small molecule docking simulations. Different algorithms integrate distinct forms of biochemical data during the docking and/or scoring stages. These so-called hybrid methods make use of receptor-based information such as nuclear magnetic resonance (NMR) restraints or small molecule-based information such as structure-activity relationships (SARs). A third class of methods directly interrogates contacts between the protein receptor and the small molecule. This work reviews the current state of using such restraints in docking simulations, evaluates their feasibility across broad systems, and identifies potential areas of algorithm development.
Collapse
Affiliation(s)
- Darwin Y Fu
- Department of Chemistry Vanderbilt University Nashville, Tennessee 37235, United States
| | - Jens Meiler
- Department of Chemistry Vanderbilt University Nashville, Tennessee 37235, United States
| |
Collapse
|
7
|
The molecular basis of subtype selectivity of human kinin G-protein-coupled receptors. Nat Chem Biol 2018; 14:284-290. [PMID: 29334381 DOI: 10.1038/nchembio.2551] [Citation(s) in RCA: 64] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2017] [Accepted: 11/21/2017] [Indexed: 12/16/2022]
Abstract
G-protein-coupled receptors (GPCRs) are the most important signal transducers in higher eukaryotes. Despite considerable progress, the molecular basis of subtype-specific ligand selectivity, especially for peptide receptors, remains unknown. Here, by integrating DNP-enhanced solid-state NMR spectroscopy with advanced molecular modeling and docking, the mechanism of the subtype selectivity of human bradykinin receptors for their peptide agonists has been resolved. The conserved middle segments of the bound peptides show distinct conformations that result in different presentations of their N and C termini toward their receptors. Analysis of the peptide-receptor interfaces reveals that the charged N-terminal residues of the peptides are mainly selected through electrostatic interactions, whereas the C-terminal segments are recognized via both conformations and interactions. The detailed molecular picture obtained by this approach opens a new gateway for exploring the complex conformational and chemical space of peptides and peptide analogs for designing GPCR subtype-selective biochemical tools and drugs.
Collapse
|
8
|
Wang T, Liu H, Duan Y. Assessment of the transmembrane domain structures in GPCR Dock 2013 models. J Struct Biol 2017; 201:210-220. [PMID: 29174398 DOI: 10.1016/j.jsb.2017.11.007] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2017] [Revised: 11/15/2017] [Accepted: 11/20/2017] [Indexed: 01/08/2023]
Abstract
The community-wide blind prediction of G-protein coupled receptor (GPCR) structures and ligand docking has been conducted three times and the quality of the models was primarily assessed by the accuracy of ligand binding modes. The seven transmembrane (TM) helices of the receptors were taken as a whole; thus the model quality within the 7TM domains has not been evaluated. Here we evaluate the 7TM domain structures in the models submitted for the last round of prediction - GPCR Dock 2013. Applying the 7 × 7 RMSD matrix analysis described in our prior work, we show that the models vary widely in prediction accuracy of the 7TM structures, exhibiting diverse structural differences from the targets. For the prediction of the 5-hydroxytryptamine receptors, the top 7TM models are rather close to the targets, which however are not ranked top by ligand-docking. On the other hand, notable deviations of the TMs are found in in the previously identified top docking models that closely resemble other receptors. We further reveal reasons of success and failure in ligand docking for the models. This current assessment not only complements the previous assessment, but also provides important insights into the current status of GPCR modeling and ligand docking.
Collapse
Affiliation(s)
- Ting Wang
- Genome Center, 451 East Health Science Drive, University of California, Davis, CA 95616, United States; Sichuan University of Science and Engineering, 180 Xueyuan Street, Huixing Road, Zigong 643000, Sichuan Province, China.
| | - Haiguang Liu
- Complex Systems Division, Beijing Computational Science Research Center, 10 W. Dongbeiwang Rd, Haidian District, Beijing 100193, China
| | - Yong Duan
- Genome Center, 451 East Health Science Drive, University of California, Davis, CA 95616, United States
| |
Collapse
|
9
|
Coudrat T, Simms J, Christopoulos A, Wootten D, Sexton PM. Improving virtual screening of G protein-coupled receptors via ligand-directed modeling. PLoS Comput Biol 2017; 13:e1005819. [PMID: 29131821 PMCID: PMC5708846 DOI: 10.1371/journal.pcbi.1005819] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2017] [Revised: 11/30/2017] [Accepted: 10/12/2017] [Indexed: 11/22/2022] Open
Abstract
G protein-coupled receptors (GPCRs) play crucial roles in cell physiology and pathophysiology. There is increasing interest in using structural information for virtual screening (VS) of libraries and for structure-based drug design to identify novel agonist or antagonist leads. However, the sparse availability of experimentally determined GPCR/ligand complex structures with diverse ligands impedes the application of structure-based drug design (SBDD) programs directed to identifying new molecules with a select pharmacology. In this study, we apply ligand-directed modeling (LDM) to available GPCR X-ray structures to improve VS performance and selectivity towards molecules of specific pharmacological profile. The described method refines a GPCR binding pocket conformation using a single known ligand for that GPCR. The LDM method is a computationally efficient, iterative workflow consisting of protein sampling and ligand docking. We developed an extensive benchmark comparing LDM-refined binding pockets to GPCR X-ray crystal structures across seven different GPCRs bound to a range of ligands of different chemotypes and pharmacological profiles. LDM-refined models showed improvement in VS performance over origin X-ray crystal structures in 21 out of 24 cases. In all cases, the LDM-refined models had superior performance in enriching for the chemotype of the refinement ligand. This likely contributes to the LDM success in all cases of inhibitor-bound to agonist-bound binding pocket refinement, a key task for GPCR SBDD programs. Indeed, agonist ligands are required for a plethora of GPCRs for therapeutic intervention, however GPCR X-ray structures are mostly restricted to their inactive inhibitor-bound state. G protein-coupled receptors (GPCRs) are a major target for drug discovery. These receptors are highly dynamic membrane proteins, and have had limited tractability using with biophysical screens that are widely adopted for globular protein targets. Thus, structure-based virtual screening (SBVS) holds great promise as a complement to physical screening for rational design of novel drugs. Indeed, the increasing number of atomic-detail GPCR X-ray crystal structures has coincided with an increase in prospective SBVS studies that have identified novel compounds. However, experimentally solved GPCR structures do not meet the full demand for SBVS, as the GPCR structural landscape is incomplete, lacking both in coverage of available GPCRs, and diversity in both receptor conformations and the chemistry of co-crystalised ligands. Here we present a novel computational GPCR binding pocket refinement method that can generate predictive GPCR/ligand complexes with improved SBVS performance. This ligand-directed modeling workflow uses parallel processing and efficient algorithms to search the GPCR/ligand conformational space faster and more efficiently than the widely used protein refinement method molecular dynamics. In this study, the resulting models are evaluated both structurally, and in retrospective SBVS. We demonstrate improved performance of refined models over their starting structures in the majority of our test cases.
Collapse
Affiliation(s)
- Thomas Coudrat
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - John Simms
- School of Life and Health Sciences, Aston University, Birmingham, United Kingdom
| | - Arthur Christopoulos
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Denise Wootten
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
- * E-mail: (DW); (PMS)
| | - Patrick M. Sexton
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
- * E-mail: (DW); (PMS)
| |
Collapse
|
10
|
Computational design of ligand-binding membrane receptors with high selectivity. Nat Chem Biol 2017; 13:715-723. [PMID: 28459439 PMCID: PMC5478435 DOI: 10.1038/nchembio.2371] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2016] [Accepted: 02/02/2017] [Indexed: 12/13/2022]
Abstract
Accurate modeling and design of protein-ligand
interactions have broad applications in cell, synthetic
biology and drug discovery but remain challenging without
experimental protein structures. Here we developed an
integrated protein homology modeling-ligand docking-protein
design approach that reconstructs protein-ligand binding
sites from homolog protein structures in the presence of
protein-bound ligand poses to capture conformational
selection and induced fit modes of ligand binding. In
structure modeling tests, we blindly predicted near-atomic
accuracy ligand conformations bound to G protein-coupled
receptors (GPCRs) that were rarely identified by traditional
approaches. We also quantitatively predicted the binding
selectivity of diverse ligands to
structurally-uncharacterized GPCRs. We then applied the
technique to design functional human dopamine receptors with
novel ligand binding selectivity. Most blindly predicted
ligand binding specificities closely agreed with
experimental validations. Our method should prove useful in
ligand discovery approaches and in reprogramming the ligand
binding profile of membrane receptors that remain difficult
to crystallize.
Collapse
|
11
|
Abstract
Docking, a molecular modelling method, has wide applications in identification and optimization in modern drug discovery. This chapter addresses the recent advances in the docking methodologies like fragment docking, covalent docking, inverse docking, post processing, hybrid techniques, homology modeling etc. and its protocol like searching and scoring functions. Advances in scoring functions for e.g. consensus scoring, quantum mechanics methods, clustering and entropy based methods, fingerprinting, etc. are used to overcome the limitations of the commonly used force-field, empirical and knowledge based scoring functions. It will cover crucial necessities and different algorithms of docking and scoring. Further different aspects like protein flexibility, ligand sampling and flexibility, and the performance of scoring function will be discussed.
Collapse
Affiliation(s)
- Ashwani Kumar
- Guru Jambheshwar University of Science and Technology, India
| | - Ruchika Goyal
- Guru Jambheshwar University of Science and Technology, India
| | - Sandeep Jain
- Guru Jambheshwar University of Science and Technology, India
| |
Collapse
|
12
|
Ekins S, Liebler J, Neves BJ, Lewis WG, Coffee M, Bienstock R, Southan C, Andrade CH. Illustrating and homology modeling the proteins of the Zika virus. F1000Res 2016; 5:275. [PMID: 27746901 DOI: 10.12688/f1000research.8213.1] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 02/29/2016] [Indexed: 12/28/2022] Open
Abstract
The Zika virus (ZIKV) is a flavivirus of the family Flaviviridae, which is similar to dengue virus, yellow fever and West Nile virus. Recent outbreaks in South America, Latin America, the Caribbean and in particular Brazil have led to concern for the spread of the disease and potential to cause Guillain-Barré syndrome and microcephaly. Although ZIKV has been known of for over 60 years there is very little in the way of knowledge of the virus with few publications and no crystal structures. No antivirals have been tested against it either in vitro or in vivo. ZIKV therefore epitomizes a neglected disease. Several suggested steps have been proposed which could be taken to initiate ZIKV antiviral drug discovery using both high throughput screens as well as structure-based design based on homology models for the key proteins. We now describe preliminary homology models created for NS5, FtsJ, NS4B, NS4A, HELICc, DEXDc, peptidase S7, NS2B, NS2A, NS1, E stem, glycoprotein M, propeptide, capsid and glycoprotein E using SWISS-MODEL. Eleven out of 15 models pass our model quality criteria for their further use. While a ZIKV glycoprotein E homology model was initially described in the immature conformation as a trimer, we now describe the mature dimer conformer which allowed the construction of an illustration of the complete virion. By comparing illustrations of ZIKV based on this new homology model and the dengue virus crystal structure we propose potential differences that could be exploited for antiviral and vaccine design. The prediction of sites for glycosylation on this protein may also be useful in this regard. While we await a cryo-EM structure of ZIKV and eventual crystal structures of the individual proteins, these homology models provide the community with a starting point for structure-based design of drugs and vaccines as well as a for computational virtual screening.
Collapse
Affiliation(s)
- Sean Ekins
- Collaborations in Chemistry, Fuquay-Varina, NC, USA; Collaborations Pharmaceuticals Inc., Fuquay-Varina, NC, USA; Collaborative Drug Discovery Inc, Burlingame, CA, USA
| | | | - Bruno J Neves
- LabMol - Laboratory for Molecular Modeling and Drug Design, Faculty of Pharmacy, Federal University of Goias, GO, Brazil
| | - Warren G Lewis
- Department of Medicine, Washington University School of Medicine, St Louis, MO, USA
| | - Megan Coffee
- The International Rescue Committee, New York, NY, USA
| | | | | | - Carolina H Andrade
- LabMol - Laboratory for Molecular Modeling and Drug Design, Faculty of Pharmacy, Federal University of Goias, GO, Brazil
| |
Collapse
|
13
|
Bender BJ, Cisneros A, Duran AM, Finn JA, Fu D, Lokits AD, Mueller BK, Sangha AK, Sauer MF, Sevy AM, Sliwoski G, Sheehan JH, DiMaio F, Meiler J, Moretti R. Protocols for Molecular Modeling with Rosetta3 and RosettaScripts. Biochemistry 2016; 55:4748-63. [PMID: 27490953 PMCID: PMC5007558 DOI: 10.1021/acs.biochem.6b00444] [Citation(s) in RCA: 137] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
![]()
Previously, we published an article
providing an overview of the
Rosetta suite of biomacromolecular modeling software and a series
of step-by-step tutorials [Kaufmann, K. W., et al. (2010) Biochemistry 49, 2987–2998]. The overwhelming positive
response to this publication we received motivates us to here share
the next iteration of these tutorials that feature de novo folding, comparative modeling, loop construction, protein docking,
small molecule docking, and protein design. This updated and expanded
set of tutorials is needed, as since 2010 Rosetta has been fully redesigned
into an object-oriented protein modeling program Rosetta3. Notable
improvements include a substantially improved energy function, an
XML-like language termed “RosettaScripts” for flexibly
specifying modeling task, new analysis tools, the addition of the
TopologyBroker to control conformational sampling, and support for
multiple templates in comparative modeling. Rosetta’s ability
to model systems with symmetric proteins, membrane proteins, noncanonical
amino acids, and RNA has also been greatly expanded and improved.
Collapse
Affiliation(s)
- Brian J Bender
- Department of Pharmacology, Vanderbilt University , Nashville, Tennessee 37232-6600, United States.,Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States
| | - Alberto Cisneros
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Amanda M Duran
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Jessica A Finn
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Pathology, Microbiology and Immunology, Vanderbilt University , Nashville, Tennessee 37232-2561, United States
| | - Darwin Fu
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Alyssa D Lokits
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Neuroscience Program, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Benjamin K Mueller
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Amandeep K Sangha
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Marion F Sauer
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Alexander M Sevy
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Gregory Sliwoski
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Jonathan H Sheehan
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States
| | - Frank DiMaio
- Department of Biochemistry, University of Washington , Seattle, Washington 98195, United States
| | - Jens Meiler
- Department of Pharmacology, Vanderbilt University , Nashville, Tennessee 37232-6600, United States.,Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States.,Department of Pathology, Microbiology and Immunology, Vanderbilt University , Nashville, Tennessee 37232-2561, United States.,Neuroscience Program, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Rocco Moretti
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| |
Collapse
|
14
|
Ekins S, Liebler J, Neves BJ, Lewis WG, Coffee M, Bienstock R, Southan C, Andrade CH. Illustrating and homology modeling the proteins of the Zika virus. F1000Res 2016; 5:275. [PMID: 27746901 PMCID: PMC5040154 DOI: 10.12688/f1000research.8213.2] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 08/30/2016] [Indexed: 12/20/2022] Open
Abstract
The Zika virus (ZIKV) is a flavivirus of the family
Flaviviridae, which is similar to dengue virus, yellow fever and West Nile virus. Recent outbreaks in South America, Latin America, the Caribbean and in particular Brazil have led to concern for the spread of the disease and potential to cause Guillain-Barré syndrome and microcephaly. Although ZIKV has been known of for over 60 years there is very little in the way of knowledge of the virus with few publications and no crystal structures. No antivirals have been tested against it either
in vitro or
in vivo. ZIKV therefore epitomizes a neglected disease. Several suggested steps have been proposed which could be taken to initiate ZIKV antiviral drug discovery using both high throughput screens as well as structure-based design based on homology models for the key proteins. We now describe preliminary homology models created for NS5, FtsJ, NS4B, NS4A, HELICc, DEXDc, peptidase S7, NS2B, NS2A, NS1, E stem, glycoprotein M, propeptide, capsid and glycoprotein E using SWISS-MODEL. Eleven out of 15 models pass our model quality criteria for their further use. While a ZIKV glycoprotein E homology model was initially described in the immature conformation as a trimer, we now describe the mature dimer conformer which allowed the construction of an illustration of the complete virion. By comparing illustrations of ZIKV based on this new homology model and the dengue virus crystal structure we propose potential differences that could be exploited for antiviral and vaccine design. The prediction of sites for glycosylation on this protein may also be useful in this regard. While we await a cryo-EM structure of ZIKV and eventual crystal structures of the individual proteins, these homology models provide the community with a starting point for structure-based design of drugs and vaccines as well as a for computational virtual screening.
Collapse
Affiliation(s)
- Sean Ekins
- Collaborations in Chemistry, Fuquay-Varina, NC, USA; Collaborations Pharmaceuticals Inc., Fuquay-Varina, NC, USA; Collaborative Drug Discovery Inc, Burlingame, CA, USA
| | | | - Bruno J Neves
- LabMol - Laboratory for Molecular Modeling and Drug Design, Faculty of Pharmacy, Federal University of Goias, GO, Brazil
| | - Warren G Lewis
- Department of Medicine, Washington University School of Medicine, St Louis, MO, USA
| | - Megan Coffee
- The International Rescue Committee, New York, NY, USA
| | | | | | - Carolina H Andrade
- LabMol - Laboratory for Molecular Modeling and Drug Design, Faculty of Pharmacy, Federal University of Goias, GO, Brazil
| |
Collapse
|
15
|
Mutation-Guided Unbiased Modeling of the Fat Sensor GPR119 for High-Yield Agonist Screening. Structure 2015; 23:2377-2386. [PMID: 26526849 DOI: 10.1016/j.str.2015.09.014] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2015] [Revised: 09/14/2015] [Accepted: 09/15/2015] [Indexed: 12/14/2022]
Abstract
Recent benchmark studies have demonstrated the difficulties in obtaining accurate predictions of ligand binding conformations to comparative models of G-protein-coupled receptors. We have developed a data-driven optimization protocol, which integrates mutational data and structural information from multiple X-ray receptor structures in combination with a fully flexible ligand docking protocol to determine the binding conformation of AR231453, a small-molecule agonist, in the GPR119 receptor. Resulting models converge to one conformation that explains the majority of data from mutation studies and is consistent with the structure-activity relationship for a large number of AR231453 analogs. Another key property of the refined models is their success in separating active ligands from decoys in a large-scale virtual screening. These results demonstrate that mutation-guided receptor modeling can provide predictions of practical value for describing receptor-ligand interactions and drug discovery.
Collapse
|
16
|
Spyrakis F, Cavasotto CN. Open challenges in structure-based virtual screening: Receptor modeling, target flexibility consideration and active site water molecules description. Arch Biochem Biophys 2015; 583:105-19. [DOI: 10.1016/j.abb.2015.08.002] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2015] [Revised: 08/03/2015] [Accepted: 08/03/2015] [Indexed: 01/05/2023]
|
17
|
Engelstoft MS, Norn C, Hauge M, Holliday ND, Elster L, Lehmann J, Jones RM, Frimurer TM, Schwartz TW. Structural basis for constitutive activity and agonist-induced activation of the enteroendocrine fat sensor GPR119. Br J Pharmacol 2015; 171:5774-89. [PMID: 25117266 DOI: 10.1111/bph.12877] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2014] [Revised: 08/04/2014] [Accepted: 08/06/2014] [Indexed: 12/25/2022] Open
Abstract
BACKGROUND AND PURPOSE GPR119 is a Gαs-coupled 7TM receptor activated by endogenous lipids such as oleoylethanolamide (OEA) and by the dietary triglyceride metabolite 2-monoacylglycerol. GPR119 stimulates enteroendocrine hormone and insulin secretion. But despite massive drug discovery efforts in the field, very little is known about the basic molecular pharmacology of GPR119. EXPERIMENTAL APPROACH GPR119 receptor signalling was studied in transfected cells. Mutational mapping (30 mutations in 23 positions) was performed on residues required for ligand-independent and agonist-induced GPR119 activation (AR231453 and OEA). Novel Rosetta-based receptor modelling was applied, using a composite template approach with segments from different X-ray structures and fully flexible ligand docking. KEY RESULTS The increased signalling induced by increasing the cell surface expression of GPR119 in the absence of agonist and the inhibitory effect of two synthetic inverse agonists demonstrated that GRP119 signals with a high degree of constitutive activity through the Gαs pathway. The mutational maps for AR231453 and OEA were very similar and, surprisingly, also similar to the mutational map for residues affecting the constitutive signalling - albeit with key differences. Surprisingly, almost all residues in extracellular loop-2b were important for the constitutive activity. The molecular modelling and docking demonstrated that AR231453 binds in a 'vertical' pocket in between mutational hits reaching from the centre of the receptor out to extracellular loop-2b. CONCLUSIONS AND IMPLICATIONS The high constitutive activity of GPR119 should be taken into account in future drug discovery efforts, which can now be guided by the detailed knowledge of the physiochemical properties of the extended ligand-binding pocket.
Collapse
Affiliation(s)
- M S Engelstoft
- Novo Nordisk Foundation Center for Basic Metabolic Research, Section for Metabolic Receptology and Enteroendocrinology, University of Copenhagen, Copenhagen, Denmark; Laboratory for Molecular Pharmacology, Department of Neuroscience and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| | | | | | | | | | | | | | | | | |
Collapse
|
18
|
Topiol S, Sabio M. The role of experimental and computational structural approaches in 7TM drug discovery. Expert Opin Drug Discov 2015. [DOI: 10.1517/17460441.2015.1072166] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
|
19
|
Computational fishing of new DNA methyltransferase inhibitors from natural products. J Mol Graph Model 2015; 60:43-54. [PMID: 26099696 DOI: 10.1016/j.jmgm.2015.04.010] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2015] [Revised: 03/28/2015] [Accepted: 04/22/2015] [Indexed: 12/31/2022]
Abstract
DNA methyltransferase inhibitors (DNMTis) have become an alternative for cancer therapies. However, only two DNMTis have been approved as anticancer drugs, although with some restrictions. Natural products (NPs) are a promising source of drugs. In order to find NPs with novel chemotypes as DNMTis, 47 compounds with known activity against these enzymes were used to build a LDA-based QSAR model for active/inactive molecules (93% accuracy) based on molecular descriptors. This classifier was employed to identify potential DNMTis on 800 NPs from NatProd Collection. 447 selected compounds were docked on two human DNA methyltransferase (DNMT) structures (PDB codes: 3SWR and 2QRV) using AutoDock Vina and Surflex-Dock, prioritizing according to their score values, contact patterns at 4 Å and molecular diversity. Six consensus NPs were identified as virtual hits against DNMTs, including 9,10-dihydro-12-hydroxygambogic, phloridzin, 2',4'-dihydroxychalcone 4'-glucoside, daunorubicin, pyrromycin and centaurein. This method is an innovative computational strategy for identifying DNMTis, useful in the identification of potent and selective anticancer drugs.
Collapse
|
20
|
Yuriev E, Holien J, Ramsland PA. Improvements, trends, and new ideas in molecular docking: 2012-2013 in review. J Mol Recognit 2015; 28:581-604. [PMID: 25808539 DOI: 10.1002/jmr.2471] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2014] [Revised: 01/16/2015] [Accepted: 02/05/2015] [Indexed: 12/11/2022]
Abstract
Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher-level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods.
Collapse
Affiliation(s)
- Elizabeth Yuriev
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, 3052, Australia
| | - Jessica Holien
- ACRF Rational Drug Discovery Centre and Structural Biology Laboratory, St. Vincent's Institute of Medical Research, Fitzroy, Victoria, 3065, Australia
| | - Paul A Ramsland
- Centre for Biomedical Research, Burnet Institute, Melbourne, Victoria, 3004, Australia.,Department of Surgery Austin Health, University of Melbourne, Melbourne, Victoria, 3084, Australia.,Department of Immunology, Monash University, Alfred Medical Research and Education Precinct, Melbourne, Victoria, 3004, Australia.,School of Biomedical Sciences, CHIRI Biosciences, Curtin University, Perth, Western Australia, 6845, Australia
| |
Collapse
|
21
|
Helmer D, Rink I, Dalton JAR, Brahm K, Jöst M, Nargang TM, Blum W, Wadhwani P, Brenner-Weiss G, Rapp BE, Giraldo J, Schmitz K. Rational design of a peptide capture agent for CXCL8 based on a model of the CXCL8:CXCR1 complex. RSC Adv 2015. [DOI: 10.1039/c4ra13749c] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
A CXCL8-binding peptide designed from the interaction sites of CXCR1 with CXCL8 serves as a capture agent and inhibits neutrophil migration.
Collapse
|
22
|
Beuming T, Lenselink B, Pala D, McRobb F, Repasky M, Sherman W. Docking and Virtual Screening Strategies for GPCR Drug Discovery. Methods Mol Biol 2015; 1335:251-76. [PMID: 26260606 DOI: 10.1007/978-1-4939-2914-6_17] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Progress in structure determination of G protein-coupled receptors (GPCRs) has made it possible to apply structure-based drug design (SBDD) methods to this pharmaceutically important target class. The quality of GPCR structures available for SBDD projects fall on a spectrum ranging from high resolution crystal structures (<2 Å), where all water molecules in the binding pocket are resolved, to lower resolution (>3 Å) where some protein residues are not resolved, and finally to homology models that are built using distantly related templates. Each GPCR project involves a distinct set of opportunities and challenges, and requires different approaches to model the interaction between the receptor and the ligands. In this review we will discuss docking and virtual screening to GPCRs, and highlight several refinement and post-processing steps that can be used to improve the accuracy of these calculations. Several examples are discussed that illustrate specific steps that can be taken to improve upon the docking and virtual screening accuracy. While GPCRs are a unique target class, many of the methods and strategies outlined in this review are general and therefore applicable to other protein families.
Collapse
Affiliation(s)
- Thijs Beuming
- Schrödinger, Inc., 120 West 45th Street, New York, NY, 10036, USA,
| | | | | | | | | | | |
Collapse
|
23
|
Opportunities and challenges in the discovery of allosteric modulators of GPCRs for treating CNS disorders. Nat Rev Drug Discov 2014; 13:692-708. [PMID: 25176435 DOI: 10.1038/nrd4308] [Citation(s) in RCA: 203] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Novel allosteric modulators of G protein-coupled receptors (GPCRs) are providing fundamental advances in the development of GPCR ligands with high subtype selectivity and novel modes of efficacy that have not been possible with traditional approaches. As new allosteric modulators are advancing as drug candidates, we are developing an increased understanding of the major advantages and broad range of activities that can be achieved with these agents through selective modulation of specific signalling pathways, differential effects on GPCR homodimers versus heterodimers, and other properties. This understanding creates exciting opportunities, as well as unique challenges, in the optimization of novel therapeutic agents for disorders of the central nervous system.
Collapse
|
24
|
Valentin-Hansen L, Park M, Huber T, Grunbeck A, Naganathan S, Schwartz TW, Sakmar TP. Mapping substance P binding sites on the neurokinin-1 receptor using genetic incorporation of a photoreactive amino acid. J Biol Chem 2014; 289:18045-54. [PMID: 24831006 DOI: 10.1074/jbc.m113.527085] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Substance P (SP) is a neuropeptide that mediates numerous physiological responses, including transmission of pain and inflammation through the neurokinin-1 (NK1) receptor, a G protein-coupled receptor. Previous mutagenesis studies and photoaffinity labeling using ligand analogues suggested that the binding site for SP includes multiple domains in the N-terminal (Nt) segment and the second extracellular loop (ECLII) of NK1. To map precisely the NK1 residues that interact with SP, we applied a novel receptor-based targeted photocross-linking approach. We used amber codon suppression to introduce the photoreactive unnatural amino acid p-benzoyl-l-phenylalanine (BzF) at 11 selected individual positions in the Nt tail (residues 11-21) and 23 positions in the ECLII (residues 170(C-10)-193(C+13)) of NK1. The 34 NK1 variants were expressed in mammalian HEK293 cells and retained the ability to interact with a fluorescently labeled SP analog. Notably, 10 of the receptor variants with BzF in the Nt tail and 4 of those with BzF in ECLII cross-linked efficiently to SP, indicating that these 14 sites are juxtaposed to SP in the ligand-bound receptor. These results show that two distinct regions of the NK1 receptor possess multiple determinants for SP binding and demonstrate the utility of genetically encoded photocross-linking to map complex multitopic binding sites on G protein-coupled receptors in a cell-based assay format.
Collapse
Affiliation(s)
- Louise Valentin-Hansen
- From the Laboratory for Molecular Pharmacology, Department of Neuroscience and Pharmacology, University of Copenhagen, The Panum Institute, Blegdamsvej 3, 2200 Copenhagen, Denmark, Section for Metabolic Receptology and Enteroendocrinology, Novo Nordisk Foundation Center for Basic Metabolic Research, University of Copenhagen, Blegdamsvej 3, 2200 Copenhagen, Denmark, and
| | - Minyoung Park
- Section for Metabolic Receptology and Enteroendocrinology, Novo Nordisk Foundation Center for Basic Metabolic Research, University of Copenhagen, Blegdamsvej 3, 2200 Copenhagen, Denmark, and Laboratory of Chemical Biology and Signal Transduction, The Rockefeller University, New York, New York 10065
| | - Thomas Huber
- Laboratory of Chemical Biology and Signal Transduction, The Rockefeller University, New York, New York 10065
| | - Amy Grunbeck
- Laboratory of Chemical Biology and Signal Transduction, The Rockefeller University, New York, New York 10065
| | - Saranga Naganathan
- Laboratory of Chemical Biology and Signal Transduction, The Rockefeller University, New York, New York 10065
| | - Thue W Schwartz
- From the Laboratory for Molecular Pharmacology, Department of Neuroscience and Pharmacology, University of Copenhagen, The Panum Institute, Blegdamsvej 3, 2200 Copenhagen, Denmark, Section for Metabolic Receptology and Enteroendocrinology, Novo Nordisk Foundation Center for Basic Metabolic Research, University of Copenhagen, Blegdamsvej 3, 2200 Copenhagen, Denmark, and
| | - Thomas P Sakmar
- Section for Metabolic Receptology and Enteroendocrinology, Novo Nordisk Foundation Center for Basic Metabolic Research, University of Copenhagen, Blegdamsvej 3, 2200 Copenhagen, Denmark, and Laboratory of Chemical Biology and Signal Transduction, The Rockefeller University, New York, New York 10065
| |
Collapse
|
25
|
Dalton JAR, Gómez-Santacana X, Llebaria A, Giraldo J. Computational analysis of negative and positive allosteric modulator binding and function in metabotropic glutamate receptor 5 (in)activation. J Chem Inf Model 2014; 54:1476-87. [PMID: 24793143 DOI: 10.1021/ci500127c] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Metabotropic glutamate receptors (mGluRs) are high-profile G-protein coupled receptors drug targets because of their involvement in several neurological disease states, and mGluR5 in particular is a subtype whose controlled allosteric modulation, both positive and negative, can potentially be useful for the treatment of schizophrenia and relief of chronic pain, respectively. Here we model mGluR5 with a collection of positive and negative allosteric modulators (PAMs and NAMs) in both active and inactive receptor states, in a manner that is consistent with experimental information, using a specialized protocol that includes homology to increase docking accuracy, and receptor relaxation to generate an individual induced fit with each allosteric modulator. Results implicate two residues in particular for NAM and PAM function: NAM interaction with W785 for receptor inactivation, and NAM/PAM H-bonding with S809 for receptor (in)activation. Models suggest the orientation of the H-bond between allosteric modulator and S809, controlled by PAM/NAM chemistry, influences the position of TM7, which in turn influences the shape of the allosteric site, and potentially the receptor state. NAM-bound and PAM-bound mGluR5 models also reveal that although PAMs and NAMs bind in the same pocket and share similar binding modes, they have distinct effects on the conformation of the receptor. Our models, together with the identification of a possible activation mechanism, may be useful in the rational design of new allosteric modulators for mGluR5.
Collapse
Affiliation(s)
- James A R Dalton
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Universitat Autònoma de Barcelona , 08193 Bellaterra, Barcelona, Spain
| | | | | | | |
Collapse
|
26
|
Gregory KJ, Nguyen ED, Malosh C, Mendenhall JL, Zic JZ, Bates BS, Noetzel MJ, Squire EF, Turner EM, Rook JM, Emmitte KA, Stauffer SR, Lindsley CW, Meiler J, Conn PJ. Identification of specific ligand-receptor interactions that govern binding and cooperativity of diverse modulators to a common metabotropic glutamate receptor 5 allosteric site. ACS Chem Neurosci 2014; 5:282-95. [PMID: 24528109 DOI: 10.1021/cn400225x] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
A common metabotropic glutamate receptor 5 (mGlu5) allosteric site is known to accommodate diverse chemotypes. However, the structural relationship between compounds from different scaffolds and mGlu5 is not well understood. In an effort to better understand the molecular determinants that govern allosteric modulator interactions with mGlu5, we employed a combination of site-directed mutagenesis and computational modeling. With few exceptions, six residues (P654, Y658, T780, W784, S808, and A809) were identified as key affinity determinants across all seven allosteric modulator scaffolds. To improve our interpretation of how diverse allosteric modulators occupy the common allosteric site, we sampled the wealth of mGlu5 structure-activity relationship (SAR) data available by docking 60 ligands (actives and inactives) representing seven chemical scaffolds into our mGlu5 comparative model. To spatially and chemically compare binding modes of ligands from diverse scaffolds, the ChargeRMSD measure was developed. We found a common binding mode for the modulators that placed the long axes of the ligands parallel to the transmembrane helices 3 and 7. W784 in TM6 not only was identified as a key NAM cooperativity determinant across multiple scaffolds, but also caused a NAM to PAM switch for two different scaffolds. Moreover, a single point mutation in TM5, G747V, altered the architecture of the common allosteric site such that 4-nitro-N-(1,3-diphenyl-1H-pyrazol-5-yl)benzamide (VU29) was noncompetitive with the common allosteric site. Our findings highlight the subtleties of allosteric modulator binding to mGlu5 and demonstrate the utility in incorporating SAR information to strengthen the interpretation and analyses of docking and mutational data.
Collapse
Affiliation(s)
- Karen J. Gregory
- Drug Discovery
Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC 3052, Australia
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
27
|
Wu H, Wang C, Gregory KJ, Han GW, Cho HP, Xia Y, Niswender CM, Katritch V, Meiler J, Cherezov V, Conn PJ, Stevens RC. Structure of a class C GPCR metabotropic glutamate receptor 1 bound to an allosteric modulator. Science 2014; 344:58-64. [PMID: 24603153 DOI: 10.1126/science.1249489] [Citation(s) in RCA: 397] [Impact Index Per Article: 39.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The excitatory neurotransmitter glutamate induces modulatory actions via the metabotropic glutamate receptors (mGlus), which are class C G protein-coupled receptors (GPCRs). We determined the structure of the human mGlu1 receptor seven-transmembrane (7TM) domain bound to a negative allosteric modulator, FITM, at a resolution of 2.8 angstroms. The modulator binding site partially overlaps with the orthosteric binding sites of class A GPCRs but is more restricted than most other GPCRs. We observed a parallel 7TM dimer mediated by cholesterols, which suggests that signaling initiated by glutamate's interaction with the extracellular domain might be mediated via 7TM interactions within the full-length receptor dimer. A combination of crystallography, structure-activity relationships, mutagenesis, and full-length dimer modeling provides insights about the allosteric modulation and activation mechanism of class C GPCRs.
Collapse
Affiliation(s)
- Huixian Wu
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
28
|
Affiliation(s)
- Ingemar André
- Department of Biochemistry and Structural Biology, Lund University, Lund, Sweden
| | - Jacob Corn
- Department of Early Discovery Biochemistry, Genentech Inc., South San Francisco, California, United States of America
- * E-mail:
| |
Collapse
|