the royal society academia - industry international conference 2014 focussed on the topic of  gpcr structure , function , drug discovery and crystallography and was held on september 12 in chicheley hall , uk .
this conference brought together 20 renowned experts in gpcr research and drug discovery spanning europe , australia and north america .
approximately half of the attendees were from academia and half from industry ( see fig . 
1attendees of the conference ( listed from left to right ) : front row  isabel moraes , fiona marshall , gebhard schertler and patrick sexton ; second row 
shweta singh , irina tikhonova , tom ceska , roland seifert , simon hodgson , daniel fourmy and alexander heifetz ; back row 
ian storer , mike bodkin , vadim cherezov , christofer tautermann , christopher tate , vsevolod gurevich and peter hunt .
chicheley hall itself can be seen in the background ( the photograph was taken by richard law ) attendees of the conference ( listed from left to right ) : front row  isabel moraes , fiona marshall , gebhard schertler and patrick sexton ; second row 
shweta singh , irina tikhonova , tom ceska , roland seifert , simon hodgson , daniel fourmy and alexander heifetz ; back row 
ian storer , mike bodkin , vadim cherezov , christofer tautermann , christopher tate , vsevolod gurevich and peter hunt .
chicheley hall itself can be seen in the background ( the photograph was taken by richard law ) gpcrs are a large family of integral membrane proteins that have enormous physiological and biomedical importance .
since gpcrs are involved in mediating cell signalling processes , they are implicated in many diseases and are the targets of numerous therapeutic drugs .
this is borne out by the fact that 60  % of all prescription drugs today target gpcrs ( schneberg et al .
2004 ) , developed for just 50 established gpcr targets out of the 800 known members of the gene family .
the importance of gpcrs was recently highlighted with the nobel prize for chemistry 2012 being awarded to two eminent gpcr researchers , prof .
there remains an ongoing need to better understand the interplay between structure and function of these receptors to advance our scientific knowledge and capacity to more effectively harness therapeutic capabilities . as a result
however , there is still a sizeable gap between ongoing academic research and the needs of the pharmaceutical industry .
a major reason for this is that academic and industrial scientists have too few productive opportunities to meet and interact , particularly to establish cross - discipline collaborations .
the intimate atmosphere of such a small conference provided a unique opportunity to stimulate the generation of new networks and partnerships between academia and industry , and to promote current gpcr research and its applications to drug discovery .
the invitation of scientists representing structural biology , protein engineering , pharmacology as well as computational and medicinal chemistry provided an interdisciplinary core to enable fruitful discussion and debate .
schertler pointed out that during 18  years at the mrc lab in cambridge , he had a continuous string of collaborations with industry partners from small and large pharma that created a very valuable network of complementary expertise .
this collaborative network led to breakthroughs in the expression and purification of difficult membrane protein targets , and later this network was an important ingredient in the formation of the drug discovery spin out company heptares therapeutics .
schertler s earlier industry partners became the drivers of this gpcr - oriented company with their expertise on target selection , business models and intellectual property ( ip ) complemented perfectly the expertise of the mrc academic partners .
the spin out company has in the meantime grown from three post docs to about 80 people , and it is able to tackle the most difficult gpcr targets with resources that would otherwise not have been accessible in europe from any funding agency .
the basis of good industry academia collaboration is a clear agreement about the goals of the collaboration . very often selecting a pre - competitive goal allows the academic partner to fully publish the results and gives the industry partner a significant edge in accessing emerging technologies .
the example of heptares therapeutics illustrates that a spin out company can generate resources that would neither be available inside a company or from public funding agencies .
if the spin out is able to reach milestones and refinance , then it can become a powerful tool to drive the application of new technology and lead to a change in the research culture in industry .
individuals have to commit to a longer - term perspective , which is aimed at changing the scientific landscape .
this scientific environment is essential for large companies that exist and are established to be able to recruit an excellent workforce and for small companies to pick up competitive innovation .
in multinational companies , the academic and technology environment can start to dictate the location of departments inside a large organisation . from this
follows that for a stable business development , the academic environment is as essential to the company as the in - house research activities .
one of the most significant advances in previous years has been the structural advances , both in terms of stabilising protein conformational states to make them more amenable to x - ray crystallography , but also in new technology such as x - ray free electron lasers which have the potential to accelerate structural biology not just for gpcrs but for many integral membrane protein drug targets . a significant portion of the meeting discussed these advances and also how they had been used in recent drug - discovery programs within both industry and academia .
another focus in recent times has been the concept of  ligand bias. it has become clear that individual gpcrs can exist in multiple receptor conformations and can elicit numerous functional responses , both g protein- and non - g protein - mediated .
this has led to the discovery that different ligands can stabilise distinct subsets of receptor conformations that can  traffic stimulus to diverse functional outputs with varying prominence , a concept referred to as biased signalling ( also known as functional selectivity , stimulus bias or ligand - directed signalling ) . in principle ,
biased signalling can result in the development of more efficacious and safer drugs , but there are some unresolved questions regarding the best system in which to assess these aspects .
the structural information alongside the realisation of biased signalling has also been explored with in silico modelling , and it was demonstrated at the meeting that this can give very useful insight into underlying properties of signalling control .
the final section of the meeting focussed on how to best resolve problems in the future , including modelling processes at a higher level . in the following sections , we expand on these discussions in more detail .
it is apparent from numerous studies that the stability of the gpcr - ligand complex in detergent solution is an important parameter that will dictate the success of any crystallisation trials ( tate 2012 ) .
although high thermostability alone does not guarantee the formation of diffraction - quality crystals , if the gpcr - ligand complex is too unstable , then crystals may not form or , if they do , they may diffract only poorly . the majority of gpcr structures have been determined from crystals of the receptor bound to a high - affinity antagonist , which usually binds with kd or ki values in the range of 10  pm to 10  nm . however ,
if a ligand binds to a receptor only with low affinity and/or the receptor is unstable in detergent , then it may still be possible to obtain a structure if the receptor is thermostabilised . a method developed to thermostabilise gpcrs
uses systematic scanning mutagenesis coupled to a thermostability assay performed on the detergent - solubilised mutant receptors to identify specific thermostabilising mutations ( tate 2012 ) .
each mutation usually imparts 13  c improvement in thermostability to the receptor , although the most highly stabilising mutation found improved thermostability of the agonist - bound conformation of the adenosine a2a receptor by 14  c .
once the single thermostabilising mutations have been identified , then the best thermostabilising mutations can be combined to make an optimally stable receptor ( shibata et al .
the methodology has been applied to many different gpcrs , in either an antagonist - bound conformation or an agonist - bound conformation , which have been subsequently crystallised and their structures determined .
2014 ) was a 2.1- resolution structure of the 1-adrenergic receptor ( 1ar ) bound to cyanopindolol ( see fig .  2 ) and
the crystals grown in lipidic cubic phase , although without requiring fusions of the receptor to either t4-lysozyme ( t4l ) or apocytochrome b562ril ( bril ) . the structure showed the presence of an intramembrane sodium ion that was in the identical position to the intramembrane sodium ion in the adenosine a2a receptor that is known to act as an allosteric antagonist .
in contrast , the na ion in 1ar does not appear to affect the transition between the inactive and active states of the receptor.fig .
2crystal structures of 1ar bound to novel chemotypes developed by fragment screening and hit optimisation .
structures of the ligand binding pocket are depicted with only portions of selected transmembrane helices ( h3 , h5 , h6 and h7 ) shown with the side chains ( green ) and ligands ( yellow ) depicted in sticks : a cyanopindolol , pdb code 2vt4 ; b compound 19 , pdb code 3zpq ; c compound 20 , pdb code 3zpr  figure adapted from the publication ( christopher et al .
2013 ) crystal structures of 1ar bound to novel chemotypes developed by fragment screening and hit optimisation .
structures of the ligand binding pocket are depicted with only portions of selected transmembrane helices ( h3 , h5 , h6 and h7 ) shown with the side chains ( green ) and ligands ( yellow ) depicted in sticks : a cyanopindolol , pdb code 2vt4 ; b compound 19 , pdb code 3zpq ; c compound 20 , pdb code 3zpr  figure adapted from the publication ( christopher et al .
2013 ) a significant advantage of the thermostabilised receptors is their use in drug discovery , which for the first time opens up the opportunities for structure - based drug design ( congreve et al .
the thermostabilised receptors are readily purified and crystallised on a routine basis , which facilitates co - crystallisation with fragments and lead compounds .
this was recently demonstrated for the 1ar where a fragment screen performed by surface plasmon resonance ( spr ) , followed by minimal hit optimisation , produced nm - affinity antagonists with novel scaffolds that were readily co - crystallised with the receptor ( christopher et al .
the thermostabilisation of gpcrs is a central platform in heptares therapeutics , resulting in numerous crystal structures .
another structural approach gaining momentum is the application of x - ray free - electron lasers ( xfel ) to gpcrs .
structural studies of gpcrs , and other biomedically relevant membrane proteins and complexes , are hampered by challenges related to growing sufficiently large crystals capable of withstanding radiation damage and yielding high - resolution data at synchrotron sources .
the recent introduction of a new generation of x - ray sources , x - ray free electron lasers ( xfels ) , producing ultra - bright pulses of coherent x - rays with an ultrashort duration , holds the promise to advance our understanding of structure and function of these challenging targets .
a novel approach using a membrane mimetic gel - like matrix known as lipidic cubic phase ( lcp ) for growth and delivery of membrane protein microcrystals for data collection by serial femtosecond crystallography ( sfx ) at xfels ( liu et al .
microcrystals are delivered to the intersection point with an xfel beam in random orientations using a specially designed lcp injector ( weierstall et al .
the injector allows adjusting lcp flowrate in a wide range to match the xfel pulse repetition rate and , thus , minimises crystal consumption .
lcp - sfx uses highly intense sub-50-fs xfel pulses to overcome radiation damage and collect room temperature high - resolution data from sub-10-m crystals .
protein consumption is reduced by two to three orders of magnitude compared to previously used liquid injectors , making the lcp - sfx method attractive for structural studies of challenging membrane and soluble proteins , and their complexes ( see fig . 
3schematic summary illustrating serial femtosecond crystallography of gpcrs with using lipidic cubic phase for microcrystal growth and delivery schematic summary illustrating serial femtosecond crystallography of gpcrs with using lipidic cubic phase for microcrystal growth and delivery results demonstrating the great utility of this approach were highlighted at the meeting .
they included the structure of the human delta - opioid receptor bound to a bi - functional peptide ligand ( fenalti et al .
2014 ) , the structure of a major gpcr signalling complex ( unpublished data ) and the first novel gpcr structure solved entirely by the lcp - sfx approach ( unpublished data ) . in the future
, this method could lead to the development of an efficient gpcr structure - based drug design pipeline by removing the major obstacles , such as the difficulties in preparation of large amounts of homogeneous and stable protein and growing sufficiently large crystals for a large number of different protein - ligand complexes .
another advantage of xfels is the ability to freeze protein motion in time and obtain structures of unstable intermediates , advancing our knowledge about the signal transduction mechanism in gpcrs .
it was also noted that there are many signalling complexes in the cell that have component parts that have some flexibility as part of the structure necessary for function .
quite often crystallography can stumble because only very small crystals are attainable , with very weak diffraction , and this can be due , in part , to some partial disorder in a biologically relevant part of the molecule . the crystallographic solution up to now has been to delete or modify the flexible regions of the protein in order to create constructs that are more ordered so that crystals of significant size could be grown . with diffraction from microcrystals now possible , the structures of these more challenging biological assemblies are within reach .
in terms of function , the human 2-adrenergic receptor ( 2ar ) is probably the best - studied gpcr at the molecular , cellular and physiological level ( seifert 2013 ) .
the 2ar was of critical importance for the development of current models of receptor activation including biased signalling .
recent research has also shown for this system that ligand bias can depend on the system studied ( native versus recombinant ) ( seifert 2013 ) and also on the ( patho)physiological state ( healthy versus diseased ) , a feature sometimes also referred to as dynamic bias ( michel et al .
accordingly , the analysis of ligands at receptors such as the 2ar has become much more complex , requiring multidimensional approaches .
the 2ar constitutes an important drug target ; agonists for this receptor are used for treatment of bronchial asthma and chronic - obstructive lung disease
most strikingly , the use of 2ar agonists alone in asthma is associated with increased mortality .
receptor desensitisation and activation of deleterious non - canonical signalling pathways , i.e. pathways different from the canonical gs pathway , could contribute to this situation .
moreover , specific 2ar polymorphisms may be associated with decreased responsiveness to certain ligands , and the use of racemic 2ar ligands may be problematic .
specifically , the distomeric ( not therapeutically active ) ligands may contribute to drug toxicity ( seifert and dove 2009 ) .
based on these concerns , non - canonical signalling pathways , receptor polymorphisms and pure 2ar stereoisomers were examined in a pluridimensional signal transduction matrix .
this matrix included studies with recombinant 2ar and native 2ar expressed in human neutrophils ( see fig .  4 ) .
human neutrophils constitute a relevant cell type for inflammation in bronchial asthma that can be readily isolated in substantial numbers . in neutrophils ,
the 2ar exerts anti - inflammatory effects by inhibiting chemoattractant - stimulated superoxide radical formation .
steroisomers of fenoterol were used as model ligands because these ligands have already been shown to exhibit functional selectivity ( seifert and dove 2009).fig .
the concept of functional selectivity requires that a given gpcr is analysed in multiple different assays .
gpcr - ga fusion proteins ensure a defined 1:1 stoichiometry of receptor and g - protein and allow analysis of gdp / gtp turnover , measured in the gtpase activity assay , and effector system activation , measured in the adenylyl cyclase ( ac ) assay , with high sensitivity .
the measurement of -arrestin recruitment is accomplished in hek cells stably transfected with fusion proteins of the 2ar linked to a luciferase fragment and -arrestin linked to the complementary luciferase fragment . upon binding of an agonist ,
human neutrophilic granulocytes constitute a physiologically relevant model system for the 2ar . in these cells ,
the 2ar couples to ac ( isoform 9 ) , resulting in an increase in the second messenger camp .
the 2ar inhibits formyl peptide receptor ( fpr)-mediated activation of the neutrophilic nadph oxidase ( nox ) that generates reactive oxygen species ( ros ) .
it is generally assumed that this inhibition is mediated via camp , but an increasing number of studies indicate that the inhibition is actually camp independent .
a major goal of current pharmacological research is the development of individualised pharmacotherapy that takes into account individual genetic polymorphisms . at the level of gpcrs ,
therefore , as a model receptor , the 2ar was analysed because for this gpcr several polymorphisms are already known , but assignment of specific polymorphisms to defined disease entities is controversial . after obtaining consent from volunteers and completing a questionnaire ,
a small sample of blood ( 48  ml ) was drawn from healthy male and female subjects .
a fraction of the blood was used to sequence the 2ar gene to identify known ( and unknown ) 2ar polymorphisms .
the remainder of the blood was used to isolate human neutrophils and assess the pharmacological profile of the 2ar with several standard ligands according to the signalling paradigm shown in ( a ) .
additionally , the impact of sex , age , smoking and allergy history on 2ar pharmacology is evaluated .
the study fills a gap in the field because it provides data on the pharmacological properties of gpcr polymorphisms in a physiologically relevant context .
the results of the study will be submitted for publication to a peer - reviewed journal in spring 2015 
 a multidimensional analysis of the 2ar in native and recombinant systems .
the concept of functional selectivity requires that a given gpcr is analysed in multiple different assays .
gpcr - ga fusion proteins ensure a defined 1:1 stoichiometry of receptor and g - protein and allow analysis of gdp / gtp turnover , measured in the gtpase activity assay , and effector system activation , measured in the adenylyl cyclase ( ac ) assay , with high sensitivity .
the measurement of -arrestin recruitment is accomplished in hek cells stably transfected with fusion proteins of the 2ar linked to a luciferase fragment and -arrestin linked to the complementary luciferase fragment . upon binding of an agonist ,
human neutrophilic granulocytes constitute a physiologically relevant model system for the 2ar . in these cells ,
the 2ar couples to ac ( isoform 9 ) , resulting in an increase in the second messenger camp .
the 2ar inhibits formyl peptide receptor ( fpr)-mediated activation of the neutrophilic nadph oxidase ( nox ) that generates reactive oxygen species ( ros ) .
it is generally assumed that this inhibition is mediated via camp , but an increasing number of studies indicate that the inhibition is actually camp independent .
a major goal of current pharmacological research is the development of individualised pharmacotherapy that takes into account individual genetic polymorphisms . at the level of gpcrs
therefore , as a model receptor , the 2ar was analysed because for this gpcr several polymorphisms are already known , but assignment of specific polymorphisms to defined disease entities is controversial . after obtaining consent from volunteers and completing a questionnaire , a small sample of blood ( 48  ml )
a fraction of the blood was used to sequence the 2ar gene to identify known ( and unknown ) 2ar polymorphisms .
the remainder of the blood was used to isolate human neutrophils and assess the pharmacological profile of the 2ar with several standard ligands according to the signalling paradigm shown in ( a ) .
additionally , the impact of sex , age , smoking and allergy history on 2ar pharmacology is evaluated .
the study fills a gap in the field because it provides data on the pharmacological properties of gpcr polymorphisms in a physiologically relevant context .
the results of the study will be submitted for publication to a peer - reviewed journal in spring 2015 in aggregate , these studies revealed that most reported ligands are biased towards canonical gs signalling ( reinartz et al .
this is particularly evident for ligands with large n - alkyl substituents , suggesting that these ligand domains , through constrained mobility of transmembrane helices , impede with coupling to gi proteins and -arrestin . in principle
however , the seifert group did not identify any ligand with bias towards gi or -arrestin . with regard to polymorphisms ,
the pharmacological profile of the 2ar with respect to inhibition of superoxide radical formation was assessed .
overall , substantial variability in the pharmacological profile of the 2ar in neutrophils was noted , but no association of the pharmacological 2ar profile with a specific polymorphism emerged .
thus , at the time being , there is no evidence for the notion that 2ar polymorphisms can be used to optimise asthma therapy .
to summarise this aspect , our increase in knowledge of the 2ar has resulted in a situation that renders future research more complicated .
most importantly , it is not anymore sufficient to determine a single parameter for a receptor such as gs - mediated adenylyl cyclase activation .
rather , multiple parameters have to be determined including non - canonical gi- and -arrestin signalling .
it is important to analyse the 2ar not only in recombinant but also in native systems .
it will also be very informative to resolve crystal structures of the 2ar in complex with various ligands and coupling partners to understand the molecular basis of functional selectivity ( seifert and dove 2009 ) .
the glp-1 receptor represents a good model system for studying class b receptor function ( koole et al .
glp-1 is a key incretin peptide that promotes insulin secretion in response to nutrient ingestion , but also has a range of other actions including preservation of b - cell mass , reduction in gastric emptying and reduction in appetite that make it a desirable target for treatment of type ii diabetes .
class b secretin - like receptors , like many other gpcrs , are pleiotropically coupled to a spectrum of both g - protein - dependent and -independent signalling pathways , and while camp production is the best characterised signalling endpoint for these receptors , physiological and therapeutic responses are the product of the integrated signalling response from all activated pathways
. moreover , the different contacts that are made between distinct ligands and their respective receptor can engender unique receptor conformations that give rise to distinct signalling profiles .
this behaviour can be observed through differences in activation of second messengers , but also through changes to how receptors are desensitised and down - regulated .
biased signalling is further complicated when allosteric drugs are considered , as conformational preferences of the receptor when allosteric and orthosteric ( endogenous ligand ) sites are co - occupied may be different than when either site is individually occupied .
biased signalling is particularly relevant to receptor systems that have multiple endogenous ligands , and where exogenous mimetics are used clinically , as is often seen with class b receptors .
clear ( ligand - directed ) bias for both peptides and small molecule agonists of the glp-1 was demonstrated ( koole et al .
2013a , b ) , and this is consistent with earlier studies of pituitary adenylyl cyclase - activating polypeptide type 1 receptor ( pacapr ) receptors and unclassified class b members like hctrs , suggesting that this is likely to be a common feature of agonist activation of class b receptors .
transmembrane helical packing and conformational transition involved in receptor activation are assumed to involve key hydrogen bond networks formed around polar residues in the transmembrane helices .
the importance of conserved polar residues in class b receptors at the glp-1 receptor has been recently studied ( koole et al .
this work has revealed networks of interaction that differentially contribute to global receptor activation and biased signalling . in particular , there appear to be two key interaction networks , one at the base of the receptor that may serve a similar role to the d[e]ry motif in class a receptors to maintain an inactive state ; the second is located in the mid - region of the helical core and plays a critical role in pathway specific signalling , in a ligand - dependent manner .
mutation of residues within the central network has identified differences in how ligands propagate activation transition for individual signalling pathways and that distinct ligands utilise only subsets of the network for signal propagation providing initial insight into molecular mechanisms for biased signalling .
the concept of biased signalling has also been explored in the cholecystokinin receptor-2 ( cck2r , which also binds the digestive hormone gastrin ) ( magnan et al .
2011 , 2013 ) and is a gpcr for which pharmaceutical companies and academic laboratories have successfully developed non - peptide ligands , mostly antagonists . since cck2r is a potential target in several pathologies of the central nervous system ( anxiety , panic attacks ) , of the gut ( peptic ulcer disease ) and of neuroendocrine cancers , the effects of a series of such ligands on stimulation of phospholipase - c and as well as on recruitment of non - visual arrestins and stimulation of receptor internalisation have been studied ( see fig . 
5schematic representation of the cck2r which can adopt two distinct conformational states upon cck activation .
the cck2r state couples to phospholipase - c activation and the cck2r state recruits -arrestins .
this figure also shows that gv150013x , a competitive antagonist on cck2r , is inefficient to inhibit recruitment of -arrestins by cck2r because of steric hindrance at the orthosteric binding site schematic representation of the cck2r which can adopt two distinct conformational states upon cck activation .
the cck2r state couples to phospholipase - c activation and the cck2r state recruits -arrestins .
this figure also shows that gv150013x , a competitive antagonist on cck2r , is inefficient to inhibit recruitment of -arrestins by cck2r because of steric hindrance at the orthosteric binding site several ligands stimulating phospholipase - c with efficacies reaching up to 50  % of that of cck but which were inefficient to stimulate -arrestin1/2 recruitment and receptor internalisation have been identified ( magnan et al .
these ligands , which were initially considered as full antagonists at the cck2r , are therefore more accurately described as antagonists on the -arrestin - dependent signalling pathway of this receptor , but are partial agonists of the g - protein - dependent signalling pathway .
extending these studies to ligands that are antagonists of g - protein - dependent pathways , it was discovered that the antagonist termed gv150013x ( n-(+)-[1-(adamant-1-ylmethyl)-2,4-dioxo-5-phenyl-2,3,4,5-tetrahydro-1h-1,5-benzodiazepin-3-yl]-n-phenylurea ) ( see fig . 
5 ) could not inhibit -arrestin1/2 recruitment and cck - induced cck2r internalisation ( magnan et al .
schild plot analysis of antagonist activity of gv150013x on cck - induced phospholipase - c activation indicated that this molecule competitively inhibited the cck effect .
this information led the team to dock gv150013x in the orthosteric binding site of the modelled cck2r , which had been previously validated by site - directed mutagenesis .
the docking study predicted that the absence of effect of gv150013x on cck2r recruiting -arrestins was due to the presence of a bulky adamantane moiety in the ligand .
gv150013x competitively inhibited cck - induced g - protein - dependent pathway whereas it was inefficient on the -arrestin - dependent pathway .
these data suggested that the cck2r could adopt two distinct conformations upon cck activation and these two conformations are distinguishable at the binding site level .
fourmy s lab therefore performed a site - directed mutagenesis study of the cck2r orthosteric binding site with a focus on amino acids presumably in contact with the adamantane moiety of gv150013x . in parallel , the team synthesised an analogue of gv150013x in which the adamantane moiety was substituted by a methyl group : 1-(2,3,4,5-tetrahydro-1-methyl-2,4-dioxo-5-phenyl-1h - benzo[b]diazepin-3-yl)-3-phenylurea , termed gv - ch3
. pharmacological studies with cck2r mutants and with gv - ch3 consistently demonstrated that the prediction of modelling and docking study was correct : the absence of effect of gv150013x on recruitment of -arrestins to the cck2r was due to a steric hindrance within the binding site which impedes binding of the antagonist ( magnan et al . 2013 ) .
overall , this area of research shows how pharmacological analysis of gpcr signalling combined with molecular modelling of gpcrs and chemistry of ligands can be used to analyse the origin of biased signalling .
such a strategy together with forthcoming determination of gpcr structures in complex with various signalling proteins ( g proteins , -arrestins ) opens the possibility of rational drug design of biased ligands .
one of the biased - agonism pathways involves arrestins , and they themselves present an alternative line of investigation .
arrestins specifically bind active phosphorylated gpcrs , precluding further g protein activation and channelling the signalling to g - protein - independent pathways ( gurevich and gurevich 2006 ) .
based on the elucidation of an arrestin structure and key functional elements , special arrestins to channel cell signalling in a desired direction were constructed ( gurevich and gurevich 2012 ) .
enhanced phosphorylation - independent arrestin mutants were designed by disrupting key stabilising intra - molecular interactions that hold arrestins in a basal conformation ( gurevich and gurevich 2012 ) .
enhanced arrestin-1 was shown to compensate for the defects of rhodopsin phosphorylation in vivo , prolonging the survival of mutant rod photoreceptors , improving their light sensitivity and speeding up photoresponse recovery ( song et al . 2009 ) . while the compensation with the first - generation enhanced mutant was only partial , new more powerful phosphorylation - independent forms of arrestin-1
rhodopsin - specific arrestin-1 is the prevalent form , so it is clear that one needs to target this subtype to compensate for disease - causing defects of rhodopsin phosphorylation .
however , activating mutations in many gpcrs underlie various human disorders ( schneberg et al .
since the two non - visual arrestins are fairly promiscuous , interacting with hundreds of gpcr subtypes , and most cells express 525 different gpcrs , only one of which is a mutant , to use a compensational approach , one needs receptor - specific non - visual arrestins . to this end ,
the elements of non - visual arrestins that determine their receptor preference were identified ( vishnivetskiy et al . 2011 ) , and on the backbone of the most promiscuous non - visual subtype , arrestin-3 , mutants with high ( > 50-fold ) receptor specificity were created ( gimenez et al .
this finding showed that targeting individual receptors with engineered non - visual arrestins is feasible .
arrestins interact with numerous partners , organising multi - protein complexes and recruiting them to particular sub - cellular compartments ( gurevich and gurevich 2014a , b ) .
this creates the potential of constructing signalling - biased arrestins that activate or inhibit certain pathways without affecting others .
recently , gurevich and colleagues designed an arrestin-3 mutant that acts as a silent scaffold : it binds all kinases in the c - jun n - terminal kinase ( jnk ) activation cascade , but does not promote its phosphorylation ( breitman et al .
this mutant was shown to act in a dominant - negative fashion , suppressing jnk activation in the cell ( breitman et al .
multi - functionality of arrestins suggests that parts acting on particular pathways can be separated and used to modify cell signalling .
indeed , a small element of arrestin-3 that acts as a mini - scaffold , promoting jnk activation in vitro and in cells , has been identified ( zhan et al .
since arrestins play a role in numerous signalling pathways , targeted mutagenesis can yield arrestin - based molecular tools to tell the cell what to do in a language it can not disobey ( gurevich and gurevich 2014a , b ) .
despite many examples of successful marketed drugs that modulate the function of gpcrs , there remain a large number of potential therapeutically relevant gpcrs that are regarded as difficult to drug effectively .
methods used to analyse the gpcr s ligand binding sites with a view to designing ligands were reviewed and included a discussion on the importance of water molecules for mediating interactions between ligands and receptor .
heptares use their star technology to generate thermostabilised gpcrs which can be used for biophysical studies , fragment screening and determination of x - ray structures .
examples were presented from different gpcrs showing that the most potent ligands act to displace high energy or 
an overview of class b and class c structures recently solved at heptares ( see fig .  6 ) was also presented .
the corticotropin releasing factor receptor 1 ( crf1 ) x - ray structure identified a novel allosteric binding pocket deep within the transmembrane domain and illustrated why finding ligands for the orthosteric pocket has been challenging .
the structure of class c metabotropic glutamate 5 receptor ( mglu5 ) was shown and provided an explanation for the tight structure - activity relationships ( sars ) and pharmacology mode switching which have been observed for this receptor .
the conclusions were that structure - based design could in theory now be applied more broadly across the gpcrome.fig .
2014 ) complex with cp-376395 from the class b receptor crf1 and overlays of a selection of ligands from class a receptor structures present in the pdb .
the observed ligand binding positions demonstrate the spectrum of binding modes across gpcr classes ranging from extracellular orthosteric to deeper intracellular allosteric sites crystal structure of the class c mglu5-mavoglurant receptor ( bennett et al .
2014 ) complex with cp-376395 from the class b receptor crf1 and overlays of a selection of ligands from class a receptor structures present in the pdb .
the observed ligand binding positions demonstrate the spectrum of binding modes across gpcr classes ranging from extracellular orthosteric to deeper intracellular allosteric sites antihistamines are one of the most well - studied drugs ( simons and simons 2011 ) .
indeed , it is now more than a century since the discovery of histamine as an important biogenic substance and more than 50  years since the production of the first antihistamine drugs ( e.g. chlorpheniramine , diphenhydramine , hydroxyzine ) .
h1r antagonists ( antihistamines ) are widely used in the treatment of a broad range of allergic diseases like rhinitis , conjunctivitis , urticaria and non - allergic disorders like pruritus and insomnia ( simons and simons 2011 ) .
the first generation of h1r antagonists that had been introduced in the period of 19421980 demonstrated considerable side effects .
poor selectivity for the h1r and the ability to cross the blood  brain barrier ( bbb ) interfering with the histaminergic transmission ( sedation ) were among the most unwanted side effects .
recent studies have shown an increase in the number of allergic diseases , currently affecting more than 30  % of the world population ( qin 2007 ) .
hence , there is an urgent need for more effective and safe anti - allergic drugs .
the second - generation h1r antagonists , introduced in the early 1980s , had notable advantages such as being significantly more selective and non - sedating due to the lack of the ability to cross the bbb .
these second - generation , non - sedating h1r antagonists have been widely used in the treatment of allergic conditions but still demonstrated some cardiotoxic side effects , e.g. induction of torsades - de - pointes arrhythmias .
this has recently led to the development of third - generation h1r antagonists where both sedative and cardiovascular side effects were addressed ( oppenheimer and casale 2002 ; canonica and blaiss 2011 ) .
the crystal structure of the first - generation h1r antagonist doxepine bound to h1r was solved in 2011 ( shimamura et al . 2011 ) .
doxepine has been associated with a large number of different side effects that can be rationalised by its lack of h1r selectivity and being a potent binder of h2r , some members of muscarinic , serotonic and -adrenergic gpcr subfamilies and also of some protein kinases . due to the lack of doxepine selectivity , the need for additional crystal structures of h1r bound to second and third generation of antihistamines that will rationalise the selectivity
dr . moraes reported for the first time the solution of two additional holo h1r crystal structures bound to the highly selective second- and third - generation h1r antagonists : cetirizine ( gillard et al . 2002 ) and fexofenadine ( sharma et al . 2014 ) , respectively ( unpublished data ; see fig . 
cetirizine and fexofenadine are about 600-fold more selective for h1r compared with a wide panel of gpcrs and ion channels ( gillard et al .
this work has resulted from collaboration between evotec ltd and the membrane protein laboratory ( mpl  imperial college london).fig .
7crystal structure of h1r bound to highly selective second- and third - generation antihistamines : cetirizine ( left ) and fexofenadine ( right ) crystal structure of h1r bound to highly selective second- and third - generation antihistamines : cetirizine ( left ) and fexofenadine ( right ) these new structures provide insights into potency and selectivity , the key challenges for the design of new generation of h1r antagonists and provide a significant contribution to potentially aid computational guided structure - based drug discovery of new antihistamine drugs targeting h2 , h3 and h4 receptors where crystal structures are still absent .
structural information in conjunction with careful modelling can provide greater insight into the location and functional relevance of druggable binding locations including both orthosteric and allosteric sites .
several recent structural publications have provided greater clarity on the binding modes and kinetics of existing drugs in both orthosteric sites , such as bronchodilator tiotropium binding to muscarinic m3 receptor ( tautermann et al .
2013 ) , and allosteric binding sites , such as anti - viral maraviroc which acts as a negative allosteric modulator of chemokine receptor ccr5 ( kruse et al .
provision of knowledge of this type should begin to assist in design of more subtype selective ligands , especially when combined with leading edge computational techniques such as homology modelling ( storer et al .
an area where increased structural knowledge could be especially impactful is in the improved design and optimisation of selective ligands .
selectivity between subfamily members of gpcrs has proved essential yet challenging in some instances , a notable example being the design of highly selective serotonin ( 5-ht ) 5-ht2c - receptor agonists ( see fig . 
2014 ) receptor by pyrimido[4,5-d]azepines modelling of binding and activation of 5-ht2c ( storer et al .
2014 ) receptor by pyrimido[4,5-d]azepines agonism of the 5-ht2c receptor has therapeutic implications in a number of important disease areas including cns disorders and obesity .
however , the design of ligands that are selective over agonism of the closely related 5-ht2b receptor , which has been linked with irreversible heart valvulopathy , has proved highly challenging ( storer et al .
furthermore , ligand - protein interaction features that lead to receptor agonism versus antagonism versus inverse agonism are not well understood despite numerous examples of subtle ligand structural changes driving pronounced differences in functional efficacy of gpcrs ( storer et al .
this is further complicated in instances where the same ligand exerts a different functional effect in different isoforms of the same receptor .
this has been demonstrated by differential effects in the human isoform compared to the equivalent preclinical species receptors , hampering both in vivo efficacy and safety studies .
an illustrative example is the histamine h4 receptor antagonist program where an inverse agonist of the human h4-receptor was a partial agonist of rat and an antagonist of dog h4 receptor , complicating the interpretation of preclinical in vivo studies ( mowbray et al .
therefore , a greater appreciation for the structural basis for these effects could ultimately assist in prediction and systematic avoidance of similar issues in future programs .
biased signalling is of current interest to both academia and the pharmaceutical industry based on the hypothesis that it could provide improved therapeutic benefit whilst avoiding undesirable activities that unbiased signalling of certain receptors has historically caused ( tautermann 2014 ) .
this has provided motivation to revisit some gpcr targets that were previously either poorly drugged or deemed undruggable due to lack of therapeutic index over adverse events ( correll and mckittrick 2014 ) . however , despite advances in biology and chemistry providing assay methods to measure bias and clear examples of biased ligands emerging , a deeper structural understanding of protein conformational changes and interactions that lead to differential receptor signalling is still in its infancy but clearly of keen interest to the medicinal chemistry community ( violin et al .
ccr4 appears to have three different binding sites for peptide agonist and small molecule ligands ( see fig . 
complexity was further demonstrated in that two classes of chemically distinct small molecule ligands bind to different sites , both of which are allosteric modulators and also display different functional signalling ( procopiou et al .
9homology model of the ccr4 receptor , with putative multiple binding sites for peptide agonists and small molecule allosteric ligands homology model of the ccr4 receptor , with putative multiple binding sites for peptide agonists and small molecule allosteric ligands this complex picture of ccr4 resulted from extensive chemistry and pharmacological studies , and in hindsight it is clear why 
traditional screening approaches such as binding or whole cell studies did not clearly identify the range of different ligands and sites .
in one such screening approach , it was rationalised that the intracellular ligand binding site ( helix 8) in ccr4 was likely modified owing the proximal chemical modification required in the commonly used flipr format .
whilst the individual targets are well tractable individually , the design of single molecules that  fit and potently antagonise both receptors , and have added properties of broader selectivity and intranasal or oral drug properties , is extremely time consuming and challenging ( procopiou et al . 2011 ) ( see fig .  10).fig .
10one of the low - energy conformations of the extended dual h1h3r ligands docked into the homology model of h1 one of the low - energy conformations of the extended dual h1h3r ligands docked into the homology model of h1 thus , structurally enabled approaches to 
difficult gprc targets should improve the tractability for chemistry by the direct identification of novel multiple ligand binding sites and subsequent chemical optimisation of small molecule ligands to drugs .
this is an opportunity for drug discovery in a class of gpcr targets not easily accessible by previous methods . whilst many gpcr targets are chemically tractable resulting in multiple ligands and successful drugs , there are still many gpcr targets which are not amenable to  traditional approaches , such as focussed or high throughput screening , or where the natural ligands are unsuitable starting points for the design of oral drugs owing to their properties ( e.g. large , lipophilic or peptidic ) .
today , gpcr modelling is widely used in the structure - based drug discovery process .
the availability of structural information on the binding site of a targeted gpcr plays a key role in rationalisation , efficiency and cost - effectiveness of drug design .
x - ray crystallography , a traditional source of structural information , is not currently feasible for every gpcr or gpcr - ligand complex .
this situation significantly limits the ability of crystallography to impact drug discovery for gpcr targets in  real - time and hence there is an urgent need for other practical alternatives .
gpcr modelling is widely used as a practical alternative in the absence of crystallographic data , but can also provide much more useful information . today , it can address such key issues like gpcr flexibility and dynamics , ligand kinetics ( kon / koff rates ) , prediction of water positions and their role in ligand binding , calculation of the free energy of binding ( affinity ) and prediction of the effects of mutations on ligand binding , etc . some of these modelling approaches were reviewed in this meeting . the role of flexibility and dynamics in the development of antipsychotic drugs and
recent pharmacological studies revealed that clinically effective antipsychotic agents act by binding to several bioamine receptors ( roth et al .
in particular , the interaction with the serotonin ( 5-ht2a and 5-ht6 ) and dopamine ( d2 and d3 ) receptors ( group 1 ) induces cognition - enhancing effects , while the histamine ( h1 ) , 5-ht2c and 5-ht2b receptors ( group 2 ) modulation causes unwanted side effects ( selvam et al .
the attempts to develop drugs based on the one - target - one - disease paradigm have been limited ( allen and roth 2011 ) . as a result , there remains an urgent need for innovative approaches to develop new effective multi - target agents that improve patients health while reducing care costs .
ideal candidates should selectively target disease - active members of the family ( group 1 ) while not binding to members responsible for undesired side effects ( group 2 ) .
this is a very challenging goal to achieve , as the residues forming the orthosteric binding pocket , i.e. the binding site of endogenous ligands , are conserved within the receptor family , thus causing the recognition of a drug by many members .
a strategy to overcome this issue is to design allosteric drugs targeting a less conserved allosteric site , which modulates the orthosteric site , or to design bitopic drugs , which bind to both allosteric and orthosteric sites .
the recently released crystal structures of several bioamine receptors in complex with orthosteric and allosteric ligands enable the selectivity issue of antipsychotic drugs to be addressed at the molecular level .
thus ,  structural snapshots provided by crystallography can be used to explore receptor motions using computer simulations .
it is conceivable that allosteric and bitopic modulators interact with binding pockets that exist only in a subset of the receptor conformational space .
computer modelling can contribute to their identification by providing detailed insights into motions and interactions in the entire protein family and subsequently unravelling complex relationships in generated data within a reasonable time and at low cost . in academia , this approach has been undertaken with some promising results .
for example , in the tikhonova group , a computational protocol combining concepts from statistical mechanics and chemoinformatics have been developed to explore the flexibility of the bioamine receptors and identify geometrical and physicochemical properties that characterised the conformational space of the bioamine receptor family ( selvam et al .
2013 ) . figure  11 illustrates the molecular modelling steps undertaken to identify the unique pharmacophoric features of disease - active receptors.fig .
11structure - based computational protocol for selective polypharmacology  figure adapted from a recent publication ( selvam et al .
2013 ) structure - based computational protocol for selective polypharmacology  figure adapted from a recent publication ( selvam et al .
2013 ) the bioamine receptors of groups 1 and 2 have been subjected to molecular dynamics simulations in a realistic environment .
the results of simulations show substantial flexibility and its variability across the members of the receptor family . using molecular probe mapping technique combined with the volume - based tanimoto similarity measurements ,
similar and different geometrical and physicochemical properties were shown across the conformational space of the receptor family and the unique pharmacophoric features of disease - active receptors ( group 1 ) were highlighted .
the unique features are then linked with mutational and ligand structure - activity relationship data and tested in retrospective screening .
the combination of techniques used gives an efficient method to identify unique properties of the disease - related proteins on the reduced diverse conformational space and represents a novel application of existing computational methods for the investigation of structural reasons for selective polypharmacology ( selvam et al .
this protocol can be now exploited by industry for other protein families , involving in cancer and infectious diseases , which require a multi - target approach . in an industry setting , evotec ltd uses a hierarchical gpcr modelling protocol ( hgmp ) that has been developed in conjunction with the university of oxford to support structure - based drug discovery programs ( see fig . 
the hgmp generates a 3d model of gpcr structures and its complexes with small molecules by applying a set of computational methods .
hgmp involves homology modelling , followed by md simulation and flexible ensemble docking , to predict binding poses and function of ligands bound to gpcrs .
the hgmp includes a large set of unique plugins to refine the gpcr models and exclusive scoring functions like the gpcr - likeness assessment score ( glas ) to evaluate model quality ( heifetz et al .
armed with a pairwise protein comparison method ( pros ) used to cluster the structural data generated by the hgmp and to distinguish between different activation sub - states .
recently , the capabilities of hgmp have been extended by the addition of gpcr biased ligand tools .
the optimisation of hgmp has been performed by evotec ltd in real drug discovery projects.fig .
12 
 a hgmp workflow and b a model of 5-ht2c ( in red ) produced by the hgmp workflow .
the water molecules and ions are omitted from the figure for clarity 
 a hgmp workflow and b a model of 5-ht2c ( in red ) produced by the hgmp workflow .
the water molecules and ions are omitted from the figure for clarity the performance of hgmp in gpcr drug discovery projects such as mch-1r for obesity treatment ( heifetz et al .
2012 , 2013b ) and the 5-ht2c for the treatment of metabolic disorders ( tye et al . 2011
additionally , the hgmp technology was used in the solving of the two h1r crystals structures bound to the second- and third - generation antihistamines : cetirizine ( gillard et al .
2002 ) and fexofenadine ( sharma et al . 2014 ) , respectively , as described by dr .
another area that has received increasing attention over recent years has been the role of water networks and their elucidation by md simulations ( tautermann et al .
2014 ) or the a2ar ( liu et al . 2012 ) with resolutions of 1.8   ] , it has become clear that gpcrs often exhibit conserved water networks , which extend from the extracellular side to the intracellular surface .
this conserved solvent network has been implied to be crucial for signalling ( nygaard et al .
2009 ) and a highly solvated conserved allosteric sodium binding site close to the conserved npxxy motifs of class a gpcrs is postulated to be involved in -arrestin signalling ( fenalti et al .
, the consideration of water molecules in gpcr ligand design has been shown to be crucial ( bortolato et al . 2013 ) because several crystal structures show water - mediated ligand - protein interactions ( congreve et al .
therefore , a ligand always has to gain more free energy from binding to the receptor than the removal of water actually costs .
several methods are available which enable the very crude and quick estimation of the energy penalty for water displacement , and most efficiently they are used for growing a bound ligand ( bortolato et al .
they can help medicinal chemists to decide if a certain sub - pocket of the receptor can be explored by hydrophobic moieties or if a displaced water has to be substituted by an entity which has to mimic the hydrogen bond network .
obviously , these quick methods are not thorough in a sense that the estimates of the free energy are physically sound , but often they are sufficiently good for a go / no - go decision . in order to get a ( formally ) correct estimate for the change in the free energy of binding upon ligand modification , methods like thermodynamic integration have to be applied ( christ and fox 2014 ) .
these more accurate methods all rely on long md simulations of the receptor , and therefore they can also capture reorganisations in the solvent network and induced effects in the binding pocket .
one step further beyond assessing the free energy of ligand binding is the estimation of binding kinetics of ligands .
water networks also play a crucial role here as well , as the breaking of ligand - receptor hydrogen bonds usually ends up with solvated hydrogen bond acceptor and donor functional groups .
as the energy barrier of the reaction determines the rate constant , the direct interaction of water molecules with the ligand - protein hydrogen bonds is decisive for the dissociation rate .
if the hydrogen bond is solvent accessible , the breaking of the bond does not require high energy because the re - hydration occurs simultaneously .
if the hydrogen bond is buried , a breaking leads to a high energy penalty because acceptor and donor do not find immediate new interaction partners and are in unfavourable solvation states ( schmidtke et al .
currently , the in silico prediction of association / dissociation rates is only possible for very small fragments rather than drug - like molecules .
therefore , the state - of - the - art is still the explanation of experimental off - rates of drug - like molecules rather than the prediction of them .
recently , the duration of action of tiotropium on the human muscarinic m3 receptor ( hm3r ) was studied ( tautermann et al .
13 ) by methyl does not reduce the pki strongly , but it has a major effect on the receptor half - life of the molecule .
long md simulations ( > 2  s ) were used to investigate the differences in the bound species , especially focussing on the water network .
tiotropium forms a double hydrogen bond to n6.52 in hm3r ( fig .  13 , left panel ) , and during the simulation no water comes close ( blue solid surface ) . when substituting the hydroxy - group by methyl (  methyl ligand ) , water density is observed directly above the ( single ) hydrogen bond ( green mesh ) , corresponding to a large number of md snapshots where n6.52 is ( partly ) solvated
this comes along with a significantly widened exit channel of hm3r , and some md snapshots even show water inserting in the hydrogen bond , as displayed in fig . 
thus , the explanation for the unexpected change in off - rates is the fact that the shielded hydrogen bond in bound tiotropium becomes solvent exposed when modifying the ligand .
this observation was only possible through long md simulations , and an explanation by the static x - ray structure would not have been possible .
to summarise , the deep understanding of water networks within gpcrs has proven to be essential if one wants to understand gpcr signalling as well as gpcr ligand binding and dissociation .
md simulations are a very valuable tool to derive physically meaningful parameters such as free energy differences or solvent maps because they capture the dynamics of the systems , which is crucial especially for flexible proteins such as gpcrs.fig . 13 
 left 
comparison of the water ( oxygen ) densities in the tiotropium and the methyl ligand .
the water densities of the tiotropium md are displayed as blue solid surfaces ; densities in methyl - ligand mds are shown as green mesh .
the significant extra density in the methyl - ligand md is marked by a green ellipse .
top right  chemical structure of tiotropium ( r = oh ) and the methyl - ligand ( r = ch3 ) .
bottom right  a snapshot of the md simulation with the methyl ligand , where water inserts into the ligand - protein hydrogen bond ( in contrast to the tiotropium md , where such water - mediated hydrogen bonds are never observed)figure adapted from the recent publication ( tautermann et al .
2015 ) 
 left  comparison of the water ( oxygen ) densities in the tiotropium and the methyl ligand .
the water densities of the tiotropium md are displayed as blue solid surfaces ; densities in methyl - ligand mds are shown as green mesh .
the significant extra density in the methyl - ligand md is marked by a green ellipse .
top right  chemical structure of tiotropium ( r = oh ) and the methyl - ligand ( r = ch3 ) .
bottom right  a snapshot of the md simulation with the methyl ligand , where water inserts into the ligand - protein hydrogen bond ( in contrast to the tiotropium md , where such water - mediated hydrogen bonds are never observed)figure adapted from the recent publication ( tautermann et al .
2015 ) an essential part of computational modelling is to re - evaluate predictions often in light of new structural data .
researchers at novartis horsham had performed a large amount of work exploring the binding sites of allosteric inhibitors of the chemokine receptor cxcr2 , and several interesting things emerged after re - analysing this work in the light of crystallographic data ( salchow et al .
it had been determined at novartis horsham that antagonists were acting at an allosteric binding site , and the patent literature around the time of the work suggested that the cxc chemokines could have a binding site close to the intracellular face of the receptor .
one of the key residues in this binding site was proposed to be k320 as the cxcr2 antagonists commonly had an acidic functionality of some description within them and this residue was asparagine in the related cxcr1 receptor over which many of the cxcr2 antagonists had selectivity .
the mutagenesis experiments that were undertaken did show that residues at the intracellular end of the tm domains had an influence on the binding and/or potency of the antagonists whilst one proposed as critical within the tm domain did not show any influence .
however , the expected overlay of the antagonists based upon their ligand - only overlays was not reproduced in the effects seen against the various mutant receptors .
the residues proposed for mutagenesis , and the interpretation of their effects , were very dependent upon the model used to create the homology model for cxcr2 at the time , and so a review of how the latest gpcr crystal structures could have influenced this project was presented.fig .
14 
 a the proposed binding modes for the quinoxaline ( yellow carbons ) and urea ( cyan carbons ) series with the influential mutated residues in cxcr2 shown in cpk
. the consistent influence of k320 locates the acidic functionality in the antagonists , yet the varying effects of d143 on representatives of these two series suggest that the hydrophobic groups are located differently .
b overlay of the cxcr1 nmr structures from the pdb ( code 2lnl ) n320 is shown in green cpk and influential mutants from our experiments in purple cpk .
the protein backbones , in ribbon representation , are coloured from n ( blue ) to c ( red ) termini 
 a the proposed binding modes for the quinoxaline ( yellow carbons ) and urea ( cyan carbons ) series with the influential mutated residues in cxcr2 shown in cpk
. the consistent influence of k320 locates the acidic functionality in the antagonists , yet the varying effects of d143 on representatives of these two series suggest that the hydrophobic groups are located differently .
b overlay of the cxcr1 nmr structures from the pdb ( code 2lnl ) n320 is shown in green cpk and influential mutants from our experiments in purple cpk .
the protein backbones , in ribbon representation , are coloured from n ( blue ) to c ( red ) termini the nmr structures for the related cxcr1 receptor ( as reproduced in fig . 
14b ) consistently show that the critical n320 residue is not close to the other influential residues , and it would have been hard to rationalise the mutagenesis from this template .
the ccr5 crystal structure shows distortion in the tm7/helix 8 region and the critical k320 residue was mutated to a glu , but overall the helical alignment compared to rhodopsin was very similar and so would not have provided any benefit compared to the rhodopsin template used originally .
the 2ar - gs - protein structure was a marvellous achievement ( rasmussen et al .
2011 ) , but for the purposes of acting as a template for the intracellular binding area of the antagonists , it is likely to be too different to cxcr2 to be useful .
the most relevant structure could have been the cxcr4 , but here again the tm7/helix 8 region is distorted in the crystal and a cxcr2 homology model was created by using a chimeric template of this receptor and the tm7/helix 8 region from an earlier 2ar structure ( rasmussen et al .
interestingly , this model had a related cxcr2 antagonist bound in the tm domain and the major learning point from the current gpcr symposium was just how much influence changes within the tm domain affected the intracellular coupling and vice versa .
a conclusion drawn from the symposium with respect to the work completed at nibr horsham would be that more mutagenesis information from across the whole of the cxcr2 protein would have been needed in order to fully understand the binding area of these allosteric antagonists . finally ,
network - based approaches , evolutionary algorithms and predictive modelling , all areas that will become more prevalent in the future , were discussed .
omics age has brought with it huge quantities of data around diseases , targets , compounds and their effects .
interactomes can be built with the ultimate goal to understand disease networks and how they are influenced by the changes in small molecules and their properties .
the polypharmacology associated with current typical and atypical anti - psychotics is complex , and as an example , the question of how do we go about designing a novel anti - psychotic given the tools and data we have access to today was raised .
the opensource chembl space polypharmacology network viewer ( fechner et al . 2013 ) was introduced as an interactive way to review the rich pharmacology accessible in the chembl database and identify some good starting points for drug design .
the experimental polypharmacology associated with the hits can be complemented using target prediction ligand similarity - based approaches such as the similarity ensemble approach ( sea ; keiser et al .
2009 ) or broad panel - based predictive modelling approaches ( ghosh and jones 2014 ) .
predictive modelling approaches were also used to build protein target qsars that in combination with pharmacophore triplet compound similarity were used to develop a multi - objective scoring function . given a small organic fragment , an automated evolutionary design algorithm using reaction vectors
was used to grow a molecule by simultaneously optimising the multi - parameters required for the targeted phenotype polypharmacology ( patel et al .
the reaction vector design approach was extended to whole reaction sequences and ultimately reaction networks .
a gpr38 reaction network was built which exemplified that the chemistry phase space around a hit could be readily expanded to that of closely accessible molecules .
15automated multi - objective compound design using reaction vectors ( 26k reaction db and 93k reagents ) starting from piperidine and using four objectives : similarity to haloperidol and ziprasidone pharmacophores , dopamine d2 , 1b adrenergic and histamine qsar models .
the tri - cyclics generated appeared similar to known anti - pyschotics , chlorpromazine and fluphenazine automated multi - objective compound design using reaction vectors ( 26k reaction db and 93k reagents ) starting from piperidine and using four objectives : similarity to haloperidol and ziprasidone pharmacophores , dopamine d2 , 1b adrenergic and histamine qsar models .
the tri - cyclics generated appeared similar to known anti - pyschotics , chlorpromazine and fluphenazine it was also shown how gpcr targets and their interacting partners could be identified from genome wide association studies ( gwas ) using network analysis followed by the analysis of compound gene expression data to complement target disease gene expression as a strategy for network - based drug design .
multiple algorithms have been developed that in combination with omics data and gpcr structure - based design make a powerful arsenal for today s drug designer .
all participants agreed that gpcr research and drug discovery can benefit greatly from the collaboration between academia and industry .
the effectiveness of such collaborations for gpcr research and drug discovery was widely described and exemplified during this conference : the  non - profit gpcr consortium ( http://gpcrconsortium.org ) described by vadim cherezov ; the collaboration between mrc laboratory and heptares mentioned by gebhard schertler and chris tate ; the consortium between evotec ltd , oxford university uk and the membrane protein laboratory at diamond light source , uk described by alexander heifetz ; and many other successful networks like the adhesion - gpcr consortium ( agc ) or glisten ( gutierrez - de - teran 2014 ) ( gpcr - ligand interactions , structures , and transmembrane signalling : a european research network ) .
there was broad agreement with schertler s and ceska s comments that most ventures between academia and industry are dependent on track record and trust and that individuals have to commit to a longer - term perspective which is aimed at changing the scientific landscape .
there was also agreement that meetings between academia and industry such as this conference are very useful to learn about developments within each other s areas of expertise and to share the challenges whilst forging new links and networks . for future meetings of this type
, it was proposed to also include a broader mix of younger - generation scientists early in their career in order to maximise benefit to a broader community .
 g - protein coupled receptors ( gpcrs ) are the targets of over half of all prescribed drugs today . 
 the uniprot database has records for about 800 proteins classified as gpcrs , but drugs have only been developed against 50 of these . 
 thus , there is huge potential in terms of the number of targets for new therapies to be designed . 
 several breakthroughs in gpcrs biased pharmacology , structural biology , modelling and scoring have resulted in a resurgence of interest in gpcrs as drug targets . 
 therefore , an international conference , sponsored by the royal society , with world - renowned researchers from industry and academia was recently held to discuss recent progress and highlight key areas of future research needed to accelerate gpcr drug discovery . 
 several key points emerged . 
 firstly , structures for all three major classes of gpcrs have now been solved and there is increasing coverage across the gpcr phylogenetic tree . 
 this is likely to be substantially enhanced with data from x - ray free electron sources as they move beyond proof of concept . 
 secondly , the concept of biased signalling or functional selectivity is likely to be prevalent in many gpcrs , and this presents exciting new opportunities for selectivity and the control of side effects , especially when combined with increasing data regarding allosteric modulation . 
 thirdly , there will almost certainly be some gpcrs that will remain difficult targets because they exhibit complex ligand dependencies and have many metastable states rendering them difficult to resolve by crystallographic methods . 
 subtle effects within the packing of the transmembrane helices are likely to mask and contribute to this aspect , which may play a role in species dependent behaviour . 
 this is particularly important because it has ramifications for how we interpret pre - clinical data . in summary , 
 collaborative efforts between industry and academia have delivered significant progress in terms of structure and understanding of gpcrs and will be essential for resolving problems associated with the more difficult targets in the future . 
Introduction
Developments in GPCR crystallography
Improvements in the understanding of GPCR function
Challenges and solutions for GPCR drug discovery
In silico driven GPCR drug discovery
Conclusions