Search

G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.

The large amount of data that has been collected so far for G protein-coupled receptors requires machine learning (ML) approaches to fully exploit its potential. Our previous ML model based on gradient boosting used for prediction of drug affinity and selectivity for a receptor subtype was compared with explicit information on ligand-receptor interactions from induced-fit docking. Both methods have proved their usefulness in drug response predictions. Yet, their successful combination still requires allosteric/orthosteric assignment of ligands from datasets. Our ligand datasets included activities of two members of the secretin receptor family: GCGR and GLP-1R. Simultaneous activation of two or three receptors of this family by dual or triple agonists is not a typical kind of information included in compound databases. A precise allosteric/orthosteric ligand assignment requires a continuous update based on new structural and biological data. This data incompleteness remains the main obstacle for current ML methods applied to class B GPCR drug discovery. Even so, for these two class B receptors, our ligand-based ML model demonstrated high accuracy (5-fold cross-validation Q2 &gt, 0.63 and Q2 &gt, 0.67 for GLP-1R and GCGR, respectively). In addition, we performed a ligand annotation using recent cryogenic-electron microscopy (cryo-EM) and X-ray crystallographic data on small-molecule complexes of GCGR and GLP-1R. As a result, we assigned GLP-1R and GCGR actives deposited in ChEMBL to four small-molecule binding sites occupied by positive and negative allosteric modulators and a full agonist. Annotated compounds were added to our recently released repository of GPCR data.

The selective targeting of the cannabinoid receptor 2 (CB2) is crucial for the development of peripheral system-acting cannabinoid analgesics. This work aimed at computer-assisted identification of prospective CB2-selective compounds among the constituents of Cannabis Sativa. The molecular structures and corresponding binding affinities to CB1 and CB2 receptors were collected from ChEMBL. The molecular structures of Cannabis Sativa constituents were collected from a phytochemical database. The collected records were curated and applied for the development of quantitative structure-activity relationship (QSAR) models with a machine learning approach. The validated models predicted the affinities of Cannabis Sativa constituents. Four structures of CB2 were acquired from the Protein Data Bank (PDB) and the discriminatory ability of CB2-selective ligands and two sets of decoys were tested. We succeeded in developing the QSAR model by achieving Q2 5-CV &gt, 0.62. The QSAR models helped to identify three prospective CB2-selective molecules that are dissimilar to already tested compounds. In a complementary structure-based virtual screening study that used available PDB structures of CB2, the agonist-bound, Cryogenic Electron Microscopy structure of CB2 showed the best statistical performance in discriminating between CB2-active and non-active ligands. The same structure also performed best in discriminating between CB2-selective ligands from non-selective ligands.

The selective targeting of the cannabinoid receptor 2 (CB2) is crucial for the development of peripheral system-acting cannabinoid analgesics. This work aimed at computer-assisted identification of prospective CB2-selective compounds among the constituents of Cannabis Sativa. The molecular structures and corresponding binding affinities to CB1 and CB2 receptors were collected from ChEMBL. The molecular structures of Cannabis Sativa constituents were collected from a phytochemical database. The collected records were curated and applied for the development of quantitative structure-activity relationship (QSAR) models with a machine learning approach. The validated models predicted the affinities of Cannabis Sativa constituents. Four structures of CB2 were acquired from the Protein Data Bank (PDB) and the discriminatory ability of CB2-selective ligands and two sets of decoys were tested. We succeeded in developing the QSAR model by achieving Q2 5-CV > 0.62. The QSAR models helped to identify three prospective CB2-selective molecules that are dissimilar to already tested compounds. In a complementary structure-based virtual screening study that used available PDB structures of CB2, the agonist-bound, Cryogenic Electron Microscopy structure of CB2 showed the best statistical performance in discriminating between CB2-active and non-active ligands. The same structure also performed best in discriminating between CB2-selective ligands from non-selective ligands.

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. While experimentally-resolved structures undoubtedly provide valuable snapshots of specific GPCR conformational states, they give only limited information on their flexibility and dynamics associated with function. 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 requires efficient storage resources and specialized software, hence limiting the dissemination of these data to specialists in the field. Here we present the 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. We demonstrate the power of this resource by performing comparative analyses of multiple GPCR simulations on two mechanisms critical to receptor function: internal water networks and sodium ion interaction.

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.

Due to the involvement of G protein-coupled receptors (GPCRs) in most of the physiological and pathological processes in humans they have been attracting a lot of attention from pharmaceutical industry as well as from scientific community. Therefore, the need for new, high quality structures of GPCRs is enormous. The updated homology modeling service GPCRM (http://gpcrm.biomodellab.eu/) meets those expectations by greatly reducing the execution time of submissions (from days to hours/minutes) with nearly the same average quality of obtained models. Additionally, due to three different scoring functions (Rosetta, Rosetta-MP, BCL::Score) it is possible to select accurate models for the required purposes: the structure of the binding site, the transmembrane domain or the overall shape of the receptor. Currently, no other web service for GPCR modeling provides this possibility. GPCRM is continually upgraded in a semi-automatic way and the number of template structures has increased from 20 in 2013 to over 90 including structures the same receptor with different ligands which can influence the structure not only in the on/off manner. Two types of protein viewers can be used for visual inspection of obtained models. The extended sortable tables with available templates provide links to external databases and display ligand–receptor interactions in visual form.

Cladonia uncialis is a lichen species with confirmed antibacterial activity and whose genome has been recently sequenced, enabling first attempts in its functional characterization. In this work, we investigated activity of the C. uncialis acetone extract (CUE) and usnic acid (UA) enantiomers against ten clinical microbial strains causing skin infections. The results showed that CUE, containing (–)-UA and squamatic acid, assayed at the same concentrations as UA, was noticeably more active than (–)-UA alone, in its pure form. The studied CUE displayed an activity that was comparable to that of (+)-UA observed for Staphylococcus epidermidis and Enterococcus faecium (18-24 mm zone of growth inhibition), but did not display any activity against fungal strains. The CUE demonstrated low cytotoxicity against HaCaT cells, in comparison to UA enantiomers, which is important for its therapeutic use. Results of the antioxidant assay (DPPH) indicated low antioxidant activity (IC50>200 µg/mL) of CUE, while the total phenolic content was 70.36 mg Gallic Acid Equivalent/g of the dry extract.

Vasoactive intestinal peptide receptor 1 (VPAC1) is a member of a secretin-like subfamily of G protein-coupled receptors. Its endogenous neuropeptide (VIP), secreted by neurons and immune cells, modulates various physiological functions such as exocrine and endocrine secretions, immune response, smooth muscles relaxation, vasodilation, and fetal development. As a drug target, VPAC1 has been selected for therapy of inflammatory diseases but drug discovery is still hampered by lack of its crystal structure. In this study we presented the homology model of this receptor constructed with the well-known web service GPCRM. The VPAC1 model is composed of extracellular and transmembrane domains that form a complex with an endogenous hormone VIP. Using the homology model of VPAC1 the mechanism of action of potential drug candidates for VPAC1 was described. Only two series of small-molecule antagonists of confirmed biological activity for VPAC1 have been described thus far. Molecular docking and a series of molecular dynamics simulations were performed to elucidate their binding to VPAC1 and resulting antagonist effect. The presented work provides the basis for the possible binding mode of VPAC1 antagonists and determinants of their molecular recognition in the context of other class B GPCRs. Until the crystal structure of VPAC1 will be released, the presented homology model of VPAC1 can serve as a scaffold for drug discovery studies and is available from the author upon request., info:eu-repo/semantics/published

Predicting the functional preferences of the ligands was always a highly demanding task, much harder that predicting whether a ligand can bind to the receptor. This is because of significant similarities of agonists, antagonists and inverse agonists which are binding usually in the same binding site of the receptor and only small structural changes can push receptor toward a particular activation state. For G protein-coupled receptors, due to a large progress in crystallization techniques and also in receptor thermal stabilization, it was possible to obtain a large number of high-quality structures of complexes of these receptors with agonists and non-agonists. Additionally, the long-time-scale molecular dynamics simulations revealed how the activation processes of GPCRs can take place. Using both theoretical and experimental knowledge it was possible to employ many clever and sophisticated methods which can help to differentiate agonists and non-agonists, so one can interconvert them in search of the optimal drug.

Predicting the functional preferences of the ligands was always a highly demanding task, much harder that predicting whether a ligand can bind to the receptor. This is because of significant similarities of agonists, antagonists and inverse agonists which are binding usually in the same binding site of the receptor and only small structural changes can push receptor toward a particular activation state. For G protein-coupled receptors, due to a large progress in crystallization techniques and also in receptor thermal stabilization, it was possible to obtain a large number of high-quality structures of complexes of these receptors with agonists and non-agonists. Additionally, the long-time-scale molecular dynamics simulations revealed how the activation processes of GPCRs can take place. Using both theoretical and experimental knowledge it was possible to employ many clever and sophisticated methods which can help to differentiate agonists and non-agonists, so one can interconvert them in search of the optimal drug.

Background Membrane proteins are difficult targets for structure prediction due to the limited structural data deposited in Protein Data Bank. Most computational methods for membrane protein structure prediction are based on the comparative modeling. There are only few de novo methods targeting that distinct protein family. In this work an example of such de novo method was used to structurally and functionally characterize two representatives of distinct membrane proteins families of solute carrier transporters and G protein-coupled receptors. The well-known Rosetta program and one of its protocols named Broker was used in two test cases. The first case was de novo structure prediction of three N-terminal transmembrane helices of the human concentrative nucleoside transporter 3 (hCNT3) homotrimer belonging to the solute carrier 28 family of transporters (SLC28). The second case concerned the large scale refinement of transmembrane helices of a homology model of the corticotropin-releasing factor receptor 1 (CRFR1) belonging to the G protein-coupled receptors family. Results The inward-facing model of the hCNT3 homotrimer was used to propose the functional impact of its single nucleotide polymorphisms. Additionally, the 100 ns molecular dynamics simulation of the unliganded hCNT3 model confirmed its validity and revealed mobility of the selected binding site and homotrimer interface residues. The large scale refinement of transmembrane helices of the CRFR1 homology model resulted in the significant improvement of its accuracy with respect to the crystal structure of CRFR1, especially in the binding site area. Consequently, the antagonist CP-376395 could be docked with Autodock VINA to the CRFR1 model without any steric clashes. Conclusions The presented work demonstrated that Rosetta Broker can be a versatile tool for solving various issues referring to protein biology. Two distinct examples of de novo membrane protein structure prediction presented here provided important insights into three major areas of protein biology. Namely, the dynamics of the inward-facing hCNT3 homotrimer system, the structural changes of the CRFR1 receptor upon the antagonist binding and finally, the role of single nucleotide polymorphisms in both, hCNT3 and CRFR1 proteins, were investigated. Electronic supplementary material The online version of this article (doi:10.1186/s12900-017-0078-8) contains supplementary material, which is available to authorized users.

A novel versatile tool named GPCRM has been developed. It targets structure prediction of a distinct protein family of G protein-coupled receptors (GPCRs). In principle, GPCRM builds a GPCR model using a MODELLER-based homology modeling procedure. In addition, that commonly used procedure was improved by using comparison of sequence profiles, multiple template structures and the extensive loop refinement in Rosetta. We applied our method to predict a three dimensional structure of a gastric inhibitory polypeptide receptor (GIPR) from the secretin-like class B of human GPCRs. The GIPR model was also tested in an ensemble docking study in which we investigated plausible interactions of four potential antagonists with that receptor. Out of those four ligands we suggested ChEMBL_1933363 as the most potent antagonist of GIPR based on the Glide docking results.

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.

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.

The recent GPCR Dock 2013 assessment of serotonin receptor 5-HT1B and 5-HT2B, and smoothened receptor SMO targets, exposed the strengths and weaknesses of the currently used computational approaches. The test cases of 5-HT1B and 5-HT2B demonstrated that both the receptor structure and the ligand binding mode can be predicted with the atomic-detail accuracy, as long as the target-template sequence similarity is relatively high. On the other hand, the observation of a low target-template sequence similarity, e.g., between SMO from the frizzled GPCR family and members of the rhodopsin family, hampers the GPCR structure prediction and ligand docking. Indeed, in GPCR Dock 2013, accurate prediction of the SMO target was still beyond the capabilities of most research groups. Another bottleneck in the current GPCR research, as demonstrated by the 5-HT2B target, is the reliable prediction of global conformational changes induced by activation of GPCRs. In this work, we report details of our protocol used during GPCR Dock 2013. Our structure prediction and ligand docking protocol was especially successful in the case of 5-HT1B and 5-HT2B-ergotamine complexes for which we provide one of the most accurate predictions. In addition to a description of the GPCR Dock 2013 results, we propose a novel hybrid computational methodology to improve GPCR structure and function prediction. This computational methodology employs two separate rankings for filtering GPCR models. The first ranking is ligand-based while the second is based on the scoring scheme of the recently published BCL method. In this work, we prove that the use of knowledge-based potentials implemented in BCL is an efficient way to cope with major bottlenecks in the GPCR structure prediction. Thereby, we also demonstrate that the knowledge-based potentials for membrane proteins were significantly improved, because of the recent surge in available experimental structures.

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.

International audience; Despite tremendous successes of GPCR crystallography, the receptors with available structures represent only a small fraction of human GPCRs. An important role of the modeling community is to maximize structural insights for the remaining receptors and complexes. The community-wide GPCR Dock assessment was established to stimulate and monitor the progress in molecular modeling and ligand docking for GPCRs. The four targets in the present third assessment round presented new and diverse challenges for modelers, including prediction of allosteric ligand interaction and activation states in 5-hydroxytryptamine receptors 1B and 2B, and modeling by extremely distant homology for smoothened receptor. Forty-four modeling groups participated in the assessment. State-of-the-art modeling approaches achieved close-to-experimental accuracy for small rigid orthosteric ligands and models built by close homology, and they correctly predicted protein fold for distant homology targets. Predictions of long loops and GPCR activation states remain unsolved problems.

The membrane proteins are still the “Wild West” of structural biology. Although more and more membrane proteins structures are determined, their functioning is still difficult to investigate because they are fully functional only in the membranous environments. Several specific methodologies were developed to investigate various aspects of their cellular life but still they are challenging for computational methods. In this chapter we summarize the efforts made on elucidation the structural and dynamical properties of different types of membrane proteins emphasizing on those computational methods which were designed and employed particularly to study membrane proteins including their interactions in complex membranous systems.

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. Availability GPCRM server and database: http://gpcrm.biomodellab.eu

Vacuolar ATPases are a potential therapeutic target because of their involvement in a variety of severe diseases such as osteoporosis or cancer. Archazolide A (1) and related analogs have been previously identified as selective inhibitors of V-ATPases with potency down to the sub- nanomolar range. Herein we report on the determination of the ligand binding mode by a combination of molecular dock- ing, molecular dynamics simulations, and biochemical experi- ments, resulting in a sound model for the inhibitory mech- anism of this class of putative anticancer agents. The binding site of archazolides was confirmed to be located in the equatorial region of the membrane-embedded VO-rotor, as recently proposed on the basis of site-directed mutagenesis. Quantification of the bioactivity of a series of archazolide derivatives, together with the docking-derived binding mode of archazolides to the V-ATPase, revealed favorable ligand profiles, which can guide the development of a simplified archazolide analog with potential therapeutic relevance.

The Formyl Peptide Receptor 1 (FPR1) is an important chemotaxis receptor involved in various aspects of host defense and inflammatory processes. We constructed a model of FPR1 using as a novel template the chemokine receptor CXCR4 from the same branch of the phylogenetic tree of G-protein-coupled receptors. The previously employed template of rhodopsin contained a bulge at the extracellular part of TM2 which directly influenced binding of ligands. We also conducted molecular dynamics (MD) simulations of FPR1 in the apo form as well as in a form complexed with the agonist fMLF and the antagonist tBocMLF in the model membrane. During all MD simulation of the fMLF-FPR1 complex a water molecule transiently bridged the hydrogen bond between W254(6.48) and N108(3.35) in the middle of the receptor. We also observed a change in the cytoplasmic part of FPR1 of a rotamer of the Y301(7.53) residue (tyrosine rotamer switch). This effect facilitated movement of more water molecules toward the receptor center. Such rotamer of Y301(7.53) was not observed in any crystal structures of GPCRs which can suggest that this state is temporarily formed to pass the water molecules during the activation process. The presence of a distance between agonist and residues R201(5.38) and R205(5.42) on helix TM5 may suggest that the activation of FPR1 is similar to the activation of β-adrenergic receptors since their agonists are separated from serine residues on helix TM5. The removal of water molecules bridging these interactions in FPR1 can result in shrinking of the binding site during activation similarly to the shrinking observed in β-ARs. The number of GPCR crystal structures with agonists is still scarce so the designing of new ligands with agonistic properties is hampered, therefore homology modeling and docking can provide suitable models. Additionally, the MD simulations can be beneficial to outline the mechanisms of receptor activation and the agonist/antagonist sensing.

The community-wide GPCR Dock assessment is conducted to evaluate the status of molecular modeling and ligand docking for human G protein-coupled receptors. The present round of the assessment was based on the recent structures of dopamine D3 and CXCR4 chemokine receptors bound to small molecule antagonists and CXCR4 with a synthetic cyclopeptide. Thirty-five groups submitted their receptor-ligand complex structure predictions prior to the release of the crystallographic coordinates. With closely related homology modeling templates, as for dopamine D3 receptor, and with incorporation of biochemical and QSAR data, modern computational techniques predicted complex details with accuracy approaching experimental. In contrast, CXCR4 complexes that had less-characterized interactions and only distant homology to the known GPCR structures still remained very challenging. The assessment results provide guidance for modeling and crystallographic communities in method development and target selection for further expansion of the structural coverage of the GPCR universe. © 2011 Elsevier Ltd. All rights reserved.

Background Several different methods for contact prediction succeeded within the Sixth Critical Assessment of Techniques for Protein Structure Prediction (CASP6). The most relevant were non-local contact predictions for targets from the most difficult categories: fold recognition-analogy and new fold. Such contacts could provide valuable structural information in case a template structure cannot be found in the PDB. Results We described comprehensive tests of the effectiveness of contact data in various aspects of de novo modeling with CABS, an algorithm which was used successfully in CASP6 by the Kolinski-Bujnicki group. We used the predicted contacts in a simple scoring function for the post-simulation ranking of protein models and as a soft bias in the folding simulations and in the fold-refinement procedure. The latter approach turned out to be the most successful. The CABS force field used in the Replica Exchange Monte Carlo simulations cooperated with the true contacts and discriminated the false ones, which resulted in an improvement of the majority of Kolinski-Bujnicki's protein models. In the modeling we tested different sets of predicted contact data submitted to the CASP6 server. According to our results, the best performing were the contacts with the accuracy balanced with the coverage, obtained either from the best two predictors only or by a consensus from as many predictors as possible. Conclusion Our tests have shown that theoretically predicted contacts can be very beneficial for protein structure prediction. Depending on the protein modeling method, a contact data set applied should be prepared with differently balanced coverage and accuracy of predicted contacts. Namely, high coverage of contact data is important for the model ranking and high accuracy for the folding simulations.

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., Author Summary The activation of G-protein-coupled receptors (GPCRs) depends on small differences in agonist and antagonist structures resulting in specific forces they impose on the helical bundle of the receptor. Having the crystal structures of GPCRs in different stages of activation it is possible to investigate the successive conformational changes leading to full activation. The long molecular dynamics simulations can fill the gap spanning between those structures and provide an overview of the activation processes. The water molecules are recognized to be crucial in the activation process which link shifting of ligand in the binding site, the actions of molecular switches and finally the movements of fragments of TM helices. Here, we present five 700 ns MD simulations of lipid S1P1 receptor, either in Apo form, or bound to antagonist ML056 or natural agonist S1P. The Apo and antagonist-bound receptor structures exhibited similar behavior, with their TM bundles nearly intact, while in the case of the agonist-bound receptor we observed movements of intracellular ends of some of TM helices.

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.