ABSTRACT
G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.
Subject(s)
Molecular Dynamics Simulation , Receptors, G-Protein-Coupled/chemistry , Software , Metabolome , Models, Molecular , Protein ConformationABSTRACT
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
ABSTRACT
GPCR-ModSim (http://open.gpcr-modsim.org) is a centralized and easy to use service dedicated to the structural modeling of G-protein Coupled Receptors (GPCRs). 3D molecular models can be generated from amino acid sequence by homology-modeling techniques, considering different receptor conformations. GPCR-ModSim includes a membrane insertion and molecular dynamics (MD) equilibration protocol, which can be used to refine the generated model or any GPCR structure uploaded to the server, including if desired non-protein elements such as orthosteric or allosteric ligands, structural waters or ions. We herein revise the main characteristics of GPCR-ModSim and present new functionalities. The templates used for homology modeling have been updated considering the latest structural data, with separate profile structural alignments built for inactive, partially-active and active groups of templates. We have also added the possibility to perform multiple-template homology modeling in a unique and flexible way. Finally, our new MD protocol considers a series of distance restraints derived from a recently identified conserved network of helical contacts, allowing for a smoother refinement of the generated models which is particularly advised when there is low homology to the available templates. GPCR- ModSim has been tested on the GPCR Dock 2013 competition with satisfactory results.
Subject(s)
Internet , Models, Molecular , Receptors, G-Protein-Coupled/chemistry , Software , Algorithms , Allosteric Regulation , Amino Acid Sequence , Humans , Ligands , Molecular Dynamics Simulation , Receptor, Angiotensin, Type 2/chemistryABSTRACT
Agonists of the angiotensin II receptor type 2 (AT2), a G-protein coupled receptor, promote tissue protective effects in cardiovascular and renal diseases, while antagonists reduce neuropathic pain. We here report detailed molecular models that explain the AT2 receptor selectivity of our recent series of non-peptide ligands. In addition, minor structural changes of these ligands that provoke different functional activity are rationalized at a molecular level, and related to the selectivity for the different receptor conformations. These findings should pave the way to structure based drug discovery of AT2 receptor ligands.
Subject(s)
Angiotensin II Type 1 Receptor Blockers/chemistry , Receptor, Angiotensin, Type 2/metabolism , Amino Acid Sequence , Angiotensin II Type 1 Receptor Blockers/metabolism , Binding Sites , Ligands , Molecular Docking Simulation , Protein Binding , Protein Structure, Tertiary , Receptor, Angiotensin, Type 1/chemistry , Receptor, Angiotensin, Type 1/metabolism , Receptor, Angiotensin, Type 2/agonists , Receptor, Angiotensin, Type 2/chemistry , Sequence AlignmentABSTRACT
Two series of N-(heteroaryl)thiophene sulfonamides, encompassing either a methylene imidazole group or a tert-butylimidazolylacetyl group in the meta position of the benzene ring, have been synthesized. An AT2R selective ligand with a Ki of 42 nM was identified in the first series and in the second series, six AT2R selective ligands with significantly improved binding affinities and Ki values of <5 nM were discovered. The binding modes to AT2R were explored by docking calculations combined with molecular dynamics simulations. Although some of the high affinity ligands exhibited fair stability in human liver microsomes, comparable to that observed with C21 undergoing clinical trials, most ligands displayed a very low metabolic stability with t½ of less than 10 min in human liver microsomes. The most promising ligand, with an AT2R Ki value of 4.9 nM and with intermediate stability in human hepatocytes (t½ = 77 min) caused a concentration-dependent vasorelaxation of pre-contracted mouse aorta.
Subject(s)
Receptor, Angiotensin, Type 2 , Sulfonamides , Mice , Humans , Animals , Receptor, Angiotensin, Type 2/metabolism , Ligands , Sulfonamides/chemistry , Thiophenes/chemistry , Aorta/metabolism , Angiotensin II/metabolismABSTRACT
Angiotensin II receptor type 1 and 2 (AT1R and AT2R) are two G-protein coupled receptors that mediate most biological functions of the octapeptide Angiotensin II (Ang II). AT2R is upregulated upon tissue damage and its activation by selective AT2R agonists has become a promising approach in the search for new classes of pharmaceutical agents. We herein analyzed the chemical evolution of AT2R agonists starting from octapeptides, through shorter peptides and peptidomimetics to the first drug-like AT2R-selective agonist, C21, which is in Phase II clinical trials and aimed for idiopathic pulmonary fibrosis. Based on the recent crystal structures of AT1R and AT2R in complex with sarile, we identified a common binding model for a series of 11 selected AT2R agonists, consisting of peptides and peptidomimetics of different length, affinity towards AT2R and selectivity versus AT1R. Subsequent molecular dynamics simulations and free energy perturbation (FEP) calculations of binding affinities allowed the identification of the bioactive conformation and common pharmacophoric points, responsible for the key interactions with the receptor, which are maintained by the drug-like agonists. The results of this study should be helpful and facilitate the search for improved and even more potent AT2R-selective drug-like agonists.
Subject(s)
Angiotensin II/pharmacology , Peptides/pharmacology , Peptidomimetics/pharmacology , Receptor, Angiotensin, Type 2/agonists , Angiotensin II/chemistry , Binding Sites , Molecular Docking Simulation , Molecular Dynamics Simulation , Peptides/chemistry , Peptidomimetics/chemistry , Receptor, Angiotensin, Type 1/agonists , Receptor, Angiotensin, Type 1/chemistry , Receptor, Angiotensin, Type 2/chemistry , ThermodynamicsABSTRACT
The extended classic ternary complex model predicts that a G protein-coupled receptor (GPCR) exists in only two interconvertible states: an inactive R, and an active R(*). However, different structural active R(*) complexes may exist in addition to a silent inactive R ground state (Rg). Here we demonstrate, in a cellular context, that several R(*) states of 5-hydroxytryptamine-4 (5-HT(4)) receptors involve different side-chain conformational toggle switches. Using site-directed mutagenesis and molecular modeling approaches, we show that the basal constitutive receptor (R(*)basal) results from stabilization of an obligatory double toggle switch (Thr3.36 from inactive g- to active g+ and Trp6.48 from inactive g+ to active t). Mutation of either threonine or tryptophan to alanine resulted in a lowering of the activity of the R(*)basal similar to the Rg. The T3.36A mutation shows that the Thr3.36 toggle switch plays a minor role in the stabilization of R(*) induced by 5-HT (R(*)-5-HT) and BIMU8 (R(*)-BIMU8) and is fully required in the stabilization of R(*) induced by (S)-zacopride, cisapride, and 1-(4-amino-5-chloro-2-methoxyphenyl)-3-(1-butyl-4-piperidinyl)-1-propanone (RS 67333) (R(*)-benzamides). Thus, benzamides stabilize R(*)-benzamides by forming a specific hydrogen bond with Thr3.36 in the active g+ conformation. Conversely, R(*)-BIMU8 was probably the result of a direct conformational transition of Trp6.48 from inactive g+ to active t by hydrogen bonding of this residue to a carboxyl group of BIMU8. We were surprised that the Trp6.48 toggle switch was not necessary for receptor activation by the natural agonist 5-HT. R(*)-5-HT is probably attained through other routes of activation. Thus, different conformational arrangements occur during stabilization of R(*)basal, R(*)-5-HT, R(*)-benzamides, and R(*)-BIMU8.
Subject(s)
Receptors, Serotonin, 5-HT4/chemistry , Serotonin 5-HT4 Receptor Agonists , Serotonin Receptor Agonists/pharmacology , Amino Acid Sequence/genetics , Animals , COS Cells , Chlorocebus aethiops , Dose-Response Relationship, Drug , Mice , Molecular Sequence Data , Protein Conformation/drug effects , Protein Structure, Secondary , Receptors, Serotonin, 5-HT4/physiologyABSTRACT
We here report on our continued studies of ligands binding to the promising drug target angiotensinâ II typeâ 2 receptor (AT2R). Two series of compounds were synthesized and investigated. The first series explored the effects of adding small substituents to the phenyl ring of the known selective nonpeptide AT2R antagonist C38, generating small but significant shifts in AT2R affinity. One compound in the first series was equipotent to C38 and showed similar kinetic solubility, and stability in both human and mouse liver microsomes. The second series was comprised of new bicyclic derivatives, amongst which one ligand exhibited a five-fold improved affinity to AT2R as compared to C38. The majority of the compounds in the second series, including the most potent ligand, were inferior to C38 with regard to stability in both human and mouse microsomes. In contrast to our previously reported findings, ligands with shorter carbamate alkyl chains only demonstrated slightly improved stability in microsomes. Based on data presented herein, a more adequate, tentative model of the binding modes of ligand analogues to the prototype AT2R antagonist C38 is proposed, as deduced from docking redefined by molecular dynamic simulations.
ABSTRACT
The recent increase in available G protein-coupled receptor structures now contributes decisively to the structure-based ligand design. In this context, computational approaches in combination with medicinal chemistry and pharmacology are extremely helpful. Here, we provide an update on our structure-based computational protocols, used to answer key questions related to GPCR-ligand binding. All combined, these techniques can shed light on ligand binding modes, determine the molecular basis of conformational selection, for agonists and antagonists, as well as of subtype selectivity. To illustrate each of these questions, we will consider examples from existing projects on three families of class A (rhodopsin-like) GPCRs: one small-molecule (nucleotide-like) family, i.e., the adenosine receptors, and two peptide-binding receptors: neuropeptide-Y and angiotensin II receptors. The successful application of the same computational protocols to investigate this diverse group of receptor families gives an idea of the general applicability of our methodology in the characterization of GPCR-ligand binding.
Subject(s)
Ligands , Models, Molecular , Receptors, G-Protein-Coupled/chemistry , Receptors, G-Protein-Coupled/metabolism , Algorithms , Binding Sites , Drug Design , Humans , Molecular Docking Simulation , Molecular Dynamics Simulation , Protein Binding , Protein Conformation , Structure-Activity RelationshipABSTRACT
The family of adenosine receptors (ARs) is focus of several medicinal chemistry programs aimed to find new potent and selective drugs. Each receptor subtype has been proposed as a relevant drug target in the treatment of, e.g., cardiovascular or inflammatory diseases, asthma or Parkinson's disease. Until recently, most of these efforts have been dominated by ligand-based or empirical approaches. However, the latest advances in G protein-coupled receptor (GPCR) crystallography allowed for a thorough structural characterization of the A2AAR subtype, which has been crystalized with a number of agonists and antagonists. Consequently, the ligand discovery of AR ligands has been enriched with a number of structure-based approaches. These include the generation of higher-confident homology models for the remaining AR subtypes, virtual screening identification of novel chemotypes, structure-based lead-optimization programs, rationalization of selectivity profiles, or the structural characterization of novel binding sites that enable the design of novel allosteric modulators. Computational methodologies have importantly contributed to the success of these structure-based approaches, and the recent advances in the field are also analyzed in this review. We conclude that the design of adenosine receptor ligands has improved dramatically with the consideration of structure- based approaches, which is paving the way to a better understanding of the biology and pharmacological modulation of this relevant family of receptors.
Subject(s)
Drug Design , Receptors, Purinergic P1/chemistry , Amino Acid Sequence , Computational Biology , Humans , Ligands , Molecular Structure , Sequence Homology, Amino Acid , Structure-Activity RelationshipABSTRACT
On the basis of our previously described pharmacophore model for serotonin 5-HT(6) receptor (5-HT(6)R) antagonists, we have designed, synthesized, and pharmacologically characterized a series of benzimidazole derivatives 1-20 that represent a new family of potent antagonists at the human 5-HT(6)R. Site-directed mutagenesis and a beta(2)-adrenoceptor-based homology model of the 5-HT(6)R were used to predict the mode of binding of antagonist SB-258585 and the new synthesized ligands. Substitution of W6.48, F6.52, or N6.55 by Ala fully impedes compound 4 to block 5-HT-induced activation. Thus, we propose that D3.32 in TM 3 anchors the protonated piperazine ring, the benzimidazole ring expands parallel to EL 2 to hydrogen bond N6.55 in TM 6, and the aromatic ring is placed between TMs 3 and 5 in CH(2)-containing compounds and between TMs 3 and 6 in CO-containing compounds. This combined experimental and computational study has permitted to propose the molecular mechanisms by which the new benzimidazole derivatives act as 5-HT(6)R antagonists.
Subject(s)
Benzimidazoles/chemistry , Benzimidazoles/pharmacology , Receptors, Serotonin/chemistry , Receptors, Serotonin/metabolism , Serotonin Antagonists/chemical synthesis , Animals , Benzimidazoles/chemical synthesis , Binding, Competitive , COS Cells , Cells, Cultured , Chlorocebus aethiops , Computer Simulation , Cyclic AMP/metabolism , Enzyme-Linked Immunosorbent Assay , Humans , Kidney/cytology , Kidney/drug effects , Molecular Structure , Mutagenesis, Site-Directed , Protein Conformation , Radioligand Assay , Receptors, Serotonin/genetics , Serotonin Antagonists/chemistry , Serotonin Antagonists/pharmacology , Structure-Activity RelationshipABSTRACT
We report the synthesis of a new set of compounds of general structure I (1-20) with structural modifications in the pharmacophoric elements of the previously reported lead UCM-5600. The new derivatives have been evaluated for binding affinity at 5-HT(7) and 5-HT(1A) receptors. The influence of the different structural features in terms of 5-HT(7)/5-HT(1A) receptor affinity and selectivity was analyzed by computational simulations of the complexes between compounds I and beta(2)-based 3-D models of these receptors. Compound 18 (HYD(1) = 1,3-dihydro-2H-indol-2-one; spacer = -(CH(2))(4)-; HYD(2) + HYD(3) = 3,4-dihydroisoquinolin-2(1H)-yl) exhibits high 5-HT(7)R affinity (K(i) = 7 nM) and selectivity over the 5-HT(1A)R (31-fold), and has been characterized as a partial agonist of the human 5-HT(7)R.