G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called “molecular switches” buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homo- and heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

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 3D molecular structures of GPCRs (3D-GPCRome) during 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 to explore 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, analyse and share GPCR MD data.GPCRmd originates from a community-driven effort to create the first open, interactive, and standardized database of GPCR MD simulations.However, static high-resolution structures provide little information on the intrinsic 71 flexibility of GPCRs, a key aspect to fully understand their function. Important advances 72

The years 2000 and 2007 witnessed milestones in current understanding of G protein-coupled receptor (GPCR) structural biology. In 2000 the first GPCR, bovine rhodopsin, was crystallized and the structure was solved, while in 2007 the structure of β2-adrenergic receptor, the first GPCR with diffusible ligands, was determined owing to advances in microcrystallization and an insertion of the fast-folding lysozyme into the receptor. In parallel with those crystallographic studies, the biological and biochemical characterization of GPCRs has advanced considerably because those receptors are molecular targets for many of currently used drugs. Therefore, the mechanisms of activation and signal transduction to the cell interior deduced from known GPCRs structures are of the highest importance for drug discovery. These proteins are the most diversified membrane receptors encoded by hundreds of genes in our genome. They participate in processes responsible for vision, smell, taste and neuronal transmission in response to photons or binding of ions, hormones, peptides, chemokines and other factors. Although the GPCRs share a common seven-transmembrane α-helical bundle structure their binding sites can accommodate thousands of different ligands. The ligands, including agonists, antagonists or inverse agonists change the structure of the receptor. With bound agonists they can form a complex with a suitable G protein, be phosphorylated by kinases or bind arrestin. The discovered signaling cascades invoked by arrestin independently of G proteins makes the GPCR activating scheme more complex such that a ligand acting as an antagonist for G protein signaling can also act as an agonist in arrestin-dependent signaling. Additionally, the existence of multiple ligand-dependent partial activation states as well as dimerization of GPCRs result in a ‘microprocessor-like’ action of these receptors rather than an ‘on-off’ switch as was commonly believed only a decade ago.

A key component to success in structure-based drug design is reliable information on protein-ligand interactions. Recent development in NMR techniques has accelerated this process by overcoming some of the limitations of X-ray crystallography and computational protein-ligand docking. In this work we present a new scoring protocol based on NMR-derived interligand INPHARMA NOEs to guide the selection of computationally generated docking modes. We demonstrate the performance in a range of scenarios, encompassing traditionally difficult cases such as docking to homology models and ligand dependent domain rearrangements. Ambiguities associated with sparse experimental information are lifted by searching a consensus solution based on simultaneously fitting multiple ligand pairs. This study provides a previously unexplored integration between molecular modeling and experimental data, in which interligand NOEs represent the key element in the rescoring algorithm. The presented protocol should be widely applicable for protein-ligand docking also in a different context from drug design and highlights the important role of NMR-based approaches to describe intermolecular ligand-receptor interactions.

Sphingosine 1-phosphate (S1P) is a lysophospholipid mediator which activates G protein–coupled sphingosine 1-phosphate receptors and thus evokes a variety of cell and tissue responses including lymphocyte trafficking, endothelial development, integrity, and maturation. We performed five all-atom 700 ns molecular dynamics simulations of the sphingosine 1-phosphate receptor 1 (S1P1) based on recently released crystal structure of that receptor with an antagonist. We found that the initial movements of amino acid residues occurred in the area of highly conserved W2696.48 in TM6 which is close to the ligand binding location. Those residues located in the central part of the receptor and adjacent to kinks of TM helices comprise of a transmission switch. Side chains movements of those residues were coupled to the movements of water molecules inside the receptor which helped in the gradual opening of intracellular part of the receptor. The most stable parts of the protein were helices TM1 and TM2, while the largest movement was observed for TM7, possibly due to the short intracellular part starting with a helix kink at P7.50, which might be the first helix to move at the intracellular side. We show for the first time the detailed view of the concerted action of the transmission switch and Trp (W6.48) rotamer toggle switch leading to redirection of water molecules flow in the central part of the receptor. That event is a prerequisite for subsequent changes in intracellular part of the receptor involving water influx and opening of the receptor structure.

G-protein coupled receptors (GPCRs) are targets of nearly one third of the drugs at the current pharmaceutical market. Despite their importance in many cellular processes the crystal structures are available for less than 20 unique GPCRs of the Rhodopsin-like class. Fortunately, even though involved in different signaling cascades, this large group of membrane proteins has preserved a uniform structure comprising seven transmembrane helices that allows quite reliable comparative modeling. Nevertheless, low sequence similarity between the GPCR family members is still a serious obstacle not only in template selection but also in providing theoretical models of acceptable quality. An additional level of difficulty is the prediction of kinks and bulges in transmembrane helices. Usage of multiple templates and generation of alignments based on sequence profiles may increase the rate of success in difficult cases of comparative modeling in which the sequence similarity between GPCRs is exceptionally low. Here, we present GPCRM, a novel method for fast and accurate generation of GPCR models using averaging of multiple template structures and profile-profile comparison. In particular, GPCRM is the first GPCR structure predictor incorporating two distinct loop modeling techniques: Modeller and Rosetta together with the filtering of models based on the Z-coordinate. We tested our approach on all unique GPCR structures determined to date and report its performance in comparison with other computational methods targeting the Rhodopsin-like class. We also provide a database of precomputed GPCR models of the human receptors from that class.AvailabilityGPCRM server and database: http://gpcrm.biomodellab.eu

A disturbance of glucose homeostasis leading to type 2 diabetes mellitus (T2DM) is one of the severe side effects that may occur during a prolonged use of many drugs currently available on the market. In this manuscript we describe the most common cases of drug-induced T2DM, discuss available pharmacotherapies and propose new ones. Among various pharmacotherapies of T2DM, incretin therapies have recently focused attention due to the newly determined crystal structure of incretin hormone receptor GLP1R. Incretin hormone receptors: GLP1R and GIPR together with the glucagon receptor GCGR regulate food intake and insulin and glucose secretion. Our study showed that incretin hormone receptors, named also gut hormone receptors as they are expressed in the gastrointestinal tract, could potentially act as unintended targets (off-targets) for orally administrated drugs. Such off-target interactions, depending on their effect on the receptor (stimulation or inhibition), could be beneficial, like in the case of incretin mimetics, or unwanted if they cause, e.g., decreased insulin secretion. In this in silico study we examined which well-known pharmaceuticals could potentially interact with gut hormone receptors in the off-target way. We observed that drugs with the strongest binding affinity for gut hormone receptors were also reported in the medical information resources as the least disturbing the glucose homeostasis among all drugs in their class. We suggested that those strongly binding molecules could potentially stimulate GIPR and GLP1R and/or inhibit GCGR which could lead to increased insulin secretion and decreased hepatic glucose production. Such positive effect on the glucose homeostasis could compensate for other, adverse effects of pharmacotherapy which lead to drug-induced T2DM. In addition, we also described several top hits as potential substitutes of peptidic incretin mimetics which were discovered in the drug repositioning screen using gut hormone receptors structures against the ZINC15 compounds subset.

The prolonged use of many currently available drugs results in the severe side effect of the disruption of glucose metabolism leading to type 2 diabetes mellitus (T2DM. Gut hormone receptors including glucagon receptor (GCGR) and the incretin hormone receptors: glucagon-like peptide 1 receptor (GLP1R) and gastric inhibitory polypeptide receptor (GIPR) are important drug targets for the treatment of T2DM, as they play roles in the regulation of glucose and insulin levels and of food intake. In this study, we hypothesized that we could compensate for the negative influences of specific drugs on glucose metabolism by the positive incretin effect enhanced by the off-target interactions with incretin GPCR receptors. As a test case, we chose to examine beta-blockers because beta-adrenergic receptors and incretin receptors are expressed in a similar location, making off-target interactions possible. The binding affinity of drugs for incretin receptors was approximated by using two docking scoring functions of Autodock VINA (GUT-DOCK) and Glide (Schrodinger) and juxtaposing these values with the medical information on drug-induced T2DM. We observed that beta-blockers with the highest theoretical binding affinities for gut hormone receptors were reported as the least harmful to glucose homeostasis in clinical trials. Notably, a recently discovered beta-blocker compound 15 ([4-((2S)-3-(((S)-3-(3-bromophenyl)-1-(methylamino)-1-oxopropan-2-yl)amino)-2-(2-cyclohexyl-2-phenylacetamido)-3-oxopropyl)benzamide was among the top-scoring drugs, potentially supporting its use in the treatment of hypertension in diabetic patients. Our recently developed web service GUT-DOCK (gut-dock.miningmembrane.com) allows for the execution of similar studies for any drug-like molecule. Specifically, users can compute the binding affinities for various class B GPCRs, gut hormone receptors, VIPR1 and PAC1R.