New insights from structural biology into the druggability of G protein-coupled receptors

Efficient Refinement of Complex Structures of Flexible Histone Peptides Using Post-Docking Molecular Dynamics Protocols

Histones are keys to many epigenetic events and their complexes have therapeutic and diagnostic importance. The determination of the structures of histone complexes is fundamental in the design of new drugs. Computational molecular docking is widely used for the prediction of target-ligand complexes. Large, linear peptides like the tail regions of histones are challenging ligands for docking due to their large conformational flexibility, extensive hydration, and weak interactions with the shallow binding pockets of their reader proteins. Thus, fast docking methods often fail to produce complex structures of such peptide ligands at a level appropriate for drug design. To address this challenge, and improve the structural quality of the docked complexes, post-docking refinement has been applied using various molecular dynamics (MD) approaches. However, a final consensus has not been reached on the desired MD refinement protocol. In this present study, MD refinement strategies were systematically explored on a set of problematic complexes of histone peptide ligands with relatively large errors in their docked geometries. Six protocols were compared that differ in their MD simulation parameters. In all cases, pre-MD hydration of the complex interface regions was applied to avoid the unwanted presence of empty cavities. The best-performing protocol achieved a median of 32% improvement over the docked structures in terms of the change in root mean squared deviations from the experimental references. The influence of structural factors and explicit hydration on the performance of post-docking MD refinements are also discussed to help with their implementation in future methods and applications.

Free Energy Density of a Fluid and Its Role in Solvation and Binding

The concept that a fluid has a position-dependent free energy density appears in the literature but has not been fully developed or accepted. We set this concept on an unambiguous theoretical footing via the following strategy. First, we set forth four desiderata that should be satisfied by any definition of the position-dependent free energy density, f(R), in a system comprising only a fluid and a rigid solute: its volume integral, plus the fixed internal energy of the solute, should be the system free energy; it deviates from its bulk value, fbulk, near a solute but should asymptotically approach fbulk with increasing distance from the solute; it should go to zero where the solvent density goes to zero; and it should be well-defined in the most general case of a fluid made up of flexible molecules with an arbitrary interaction potential. Second, we use statistical thermodynamics to formulate a definition of the free energy density that satisfies these desiderata. Third, we show how any free energy density satisfying the desiderata may be used to analyze molecular processes in solution. In particular, because the spatial integral of f(R) equals the free energy of the system, it can be used to compute free energy changes that result from the rearrangement of solutes as well as the forces exerted on the solutes by the solvent. This enables the use of a thermodynamic analysis of water in protein binding sites to inform ligand design. Finally, we discuss related literature and address published concerns regarding the thermodynamic plausibility of a position-dependent free energy density. The theory presented here has applications in theoretical and computational chemistry and may be further generalizable beyond fluids, such as to solids and macromolecules.

Structure-based pharmacophore modeling 1. Automated random pharmacophore model generation

Pharmacophores are three-dimensional arrangements of molecular features required for biological activity that are often used in virtual screening efforts to prioritize ligands for experimental testing. G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for ligand discovery and drug development. Ligand-based pharmacophore models can be constructed to identify structural commonalities between known bioactive ligands for targets including GPCR. However, structure-based pharmacophores (which only require an experimentally determined or modeled structure for a protein target) have gained more attention to aid in virtual screening efforts as the number of publicly available experimentally determined GPCR structures have increased (140 unique GPCR represented as of October 24, 2022). Thus, the goal of this study was to develop a method of structure-based pharmacophore model generation applicable to ligand discovery for GPCR that have few known ligands. Pharmacophore models were generated within the active sites of 8 class A GPCR crystal structures via automated annotation of 5 randomly selected functional group fragments to sample diverse combinations of pharmacophore features. Each of the 5000 generated pharmacophores was then used to search a database containing active and decoy/inactive compounds for 30 class A GPCR and scored using enrichment factor and goodness-of-hit metrics to assess performance. Application of this method to the set of 8 class A GPCR produced pharmacophore models possessing the theoretical maximum enrichment factor value in both resolved structures (8 of 8 cases) and homology models (7 of 8 cases), indicating that generated pharmacophore models can prove useful in the context of virtual screening.

Rational design, synthesis, and pharmacological evaluation of a cohort of novel beta-adrenergic receptors ligands enables an assessment of structure-activity relationships

Biomedical applications of molecules that are able to modulate β-adrenergic signaling have become increasingly attractive over the last decade, revealing that β-adrenergic receptors (β-ARs) are key targets for a plethora of therapeutic interventions, including cancer. Despite successes in β-AR drug discovery, identification of β-AR ligands that are useful as selective chemical tools in pharmacological studies of the three β-AR subtypes, or lead compounds for drug development is still a highly challenging task. This is mainly due to the intrinsic plasticity of β-ARs as G protein-coupled receptors in conjunction with the requirement for functional receptor subtype selectivity, tissue specificity and minimal off-target effects. With the aim to provide insight into structure-activity relationships for the three β-AR subtypes, we have synthesized and obtained the pharmacological profile of a series of structurally diverse compounds (named MC) that were designed based on the aryloxy-propanolamine scaffold of SR59230A. Comparative analysis of their predicted binding mode within the active and inactive states of the receptors in combination with their pharmacological profile revealed key structural elements that control their activity as agonists or antagonists, in addition to clues about substituents that mediate selectivity for one receptor subtype over the others. We anticipate that these results will facilitate selective β-AR drug development efforts.

Membrane Lipids Are an Integral Part of Transmembrane Allosteric Sites in GPCRs: A Case Study of Cannabinoid CB1 Receptor Bound to a Negative Allosteric Modulator, ORG27569, and Analogs

A growing number of G-protein-coupled receptor (GPCR) structures reveal novel transmembrane lipid-exposed allosteric sites. Ligands must first partition into the surrounding membrane and take lipid paths to these sites. Remarkably, a significant part of the bound ligands appears exposed to the membrane lipids. The experimental structures do not usually account for the surrounding lipids, and their apparent contribution to ligand access and binding is often overlooked and poorly understood. Using classical and enhanced molecular dynamics simulations, we show that membrane lipids are critical in the access and binding of ORG27569 and its analogs at the transmembrane site of cannabinoid CB1 receptor. The observed differences in the binding affinity and cooperativity arise from the functional groups that interact primarily with lipids. Our results demonstrate the significance of incorporating membrane lipids as an integral component of transmembrane sites for accurate characterization, binding-affinity calculations, and lead optimization in drug discovery.

Membrane-facilitated receptor access and binding mechanisms of long-acting β2-adrenergic receptor (β2-AR) agonists

The drugs salmeterol, formoterol, and salbutamol constitute the frontline treatment for asthma and other chronic pulmonary diseases. These drugs activate the β2-adrenergic receptors (β2-AR), a class A G-protein-coupled receptor (GPCR) and differ significantly in their clinical onset and duration of actions. According to the "microkinetic model," the long duration of action of salmeterol and formoterol compared to salbutamol were attributed, at least in part, to their high lipophilicity and increased local concentrations in the membrane near the receptor. However, the structural and molecular bases of how the lipophilic drugs reach the binding site of the receptor from the surrounding membrane remain unknown. Using a variety of classical and enhanced molecular dynamics simulation techniques, we investigated the membrane partitioning characteristics, binding, and unbinding mechanisms of the ligands. The obtained results offer remarkable insight into the functional role of membrane lipids in the ligand association process. Strikingly, salmeterol entered the binding site from the bilayer through transmembrane helices 1 and 7. The entry was preceded by membrane-facilitated rearrangement and presentation of its phenyl-alkoxy-alkyl tail as a passkey to an access route gated by F193, a residue known critical for salmeterol's affinity. Formoterol's access is through the aqueous path shared by other β2-AR agents. We observed a novel secondary path for salbutamol that is distinct from its primary route. Our study offers a mechanistic description for the membrane-facilitated access and binding of ligands to β2-AR and establishes a groundwork for recognizing membrane lipids as an integral component in the molecular recognition process. Significance Statement The cell membrane's functional role behind the duration of action of long-acting β2-adrenergic receptor (β2-AR) agonists such as salmeterol has been a subject of debate for a long time. We investigated the binding and unbinding mechanisms of the three commonly used β2-AR agonists, salmeterol, formoterol, and salbutamol, using advanced simulation techniques. The obtained results offer unprecedented insights into the active role of membrane lipids in facilitating access and binding of the ligands, affecting the molecular recognition process and their pharmacology.

Discovery of BMS-753426: A Potent Orally Bioavailable Antagonist of CC Chemokine Receptor 2

To improve the metabolic stability profile of BMS-741672 (1a), we undertook a structure-activity relationship study in our trisubstituted cyclohexylamine series. This ultimately led to the identification of 2d (BMS-753426) as a potent and orally bioavailable antagonist of CCR2. Compared to previous clinical candidate 1a, the tert-butyl amine 2d showed significant improvements in pharmacokinetic properties, with lower clearance and higher oral bioavailability. Furthermore, compound 2d exhibited improved affinity for CCR5 and good activity in models of both monocyte migration and multiple sclerosis in the hCCR2 knock-in mouse. The synthesis of 2d was facilitated by the development of a simplified approach to key intermediate (4R)-9b that deployed a stereoselective reductive amination which may prove to be of general interest.

Sweet taste of heavy water

Hydrogen to deuterium isotopic substitution has only a minor effect on physical and chemical properties of water and, as such, is not supposed to influence its neutral taste. Here we conclusively demonstrate that humans are, nevertheless, able to distinguish D₂O from H₂O by taste. Indeed, highly purified heavy water has a distinctly sweeter taste than same-purity normal water and can add to perceived sweetness of sweeteners. In contrast, mice do not prefer D₂O over H₂O, indicating that they are not likely to perceive heavy water as sweet. HEK 293T cells transfected with the TAS1R2/TAS1R3 heterodimer and chimeric G-proteins are activated by D₂O but not by H₂O. Lactisole, which is a known sweetness inhibitor acting via the TAS1R3 monomer of the TAS1R2/TAS1R3, suppresses the sweetness of D₂O in human sensory tests, as well as the calcium release elicited by D₂O in sweet taste receptor-expressing cells. The present multifaceted experimental study, complemented by homology modelling and molecular dynamics simulations, resolves a long-standing controversy about the taste of heavy water, shows that its sweet taste is mediated by the human TAS1R2/TAS1R3 taste receptor, and opens way to future studies of the detailed mechanism of action.

Ben Abu, Mason and colleagues use molecular dynamics, cell-based experiments, mouse models, and human subjects to determine that, unlike ordinary water, heavy water tastes sweet to humans, but not mice. Mechanistically, this effect is mediated by the human TAS1R/TAS1R3 sweet taste receptor.

Sweet taste of heavy water

Hydrogen to deuterium isotopic substitution has only a minor effect on physical and chemical properties of water and, as such, is not supposed to influence its neutral taste. Here we conclusively demonstrate that humans are, nevertheless, able to distinguish D₂O from H₂O by taste. Indeed, highly purified heavy water has a distinctly sweeter taste than same-purity normal water and can add to perceived sweetness of sweeteners. In contrast, mice do not prefer D₂O over H₂O, indicating that they are not likely to perceive heavy water as sweet. HEK 293T cells transfected with the TAS1R2/TAS1R3 heterodimer and chimeric G-proteins are activated by D₂O but not by H₂O. Lactisole, which is a known sweetness inhibitor acting via the TAS1R3 monomer of the TAS1R2/TAS1R3, suppresses the sweetness of D₂O in human sensory tests, as well as the calcium release elicited by D₂O in sweet taste receptor-expressing cells. The present multifaceted experimental study, complemented by homology modelling and molecular dynamics simulations, resolves a long-standing controversy about the taste of heavy water, shows that its sweet taste is mediated by the human TAS1R2/TAS1R3 taste receptor, and opens way to future studies of the detailed mechanism of action.

Development and Application of Subtype-Selective Fluorescent Antagonists for the Study of the Human Adenosine A1 Receptor in Living Cells

The adenosine A₁ receptor (A₁AR) is a G-protein-coupled receptor (GPCR) that provides important therapeutic opportunities for a number of conditions including congestive heart failure, tachycardia, and neuropathic pain. The development of A₁AR-selective fluorescent ligands will enhance our understanding of the subcellular mechanisms underlying A₁AR pharmacology facilitating the development of more efficacious and selective therapies. Herein, we report the design, synthesis, and application of a novel series of A₁AR-selective fluorescent probes based on 8-functionalized bicyclo[2.2.2]octylxanthine and 3-functionalized 8-(adamant-1-yl) xanthine scaffolds. These fluorescent conjugates allowed quantification of kinetic and equilibrium ligand binding parameters using NanoBRET and visualization of specific receptor distribution patterns in living cells by confocal imaging and total internal reflection fluorescence (TIRF) microscopy. As such, the novel A₁AR-selective fluorescent antagonists described herein can be applied in conjunction with a series of fluorescence-based techniques to foster understanding of A₁AR molecular pharmacology and signaling in living cells.

The Effect of Deuteration on the H2 Receptor Histamine Binding Profile: A Computational Insight into Modified Hydrogen Bonding Interactions

We used a range of computational techniques to reveal an increased histamine affinity for its H₂ receptor upon deuteration, which was interpreted through altered hydrogen bonding interactions within the receptor and the aqueous environment preceding the binding. Molecular docking identified the area between third and fifth transmembrane α-helices as the likely binding pocket for several histamine poses, with the most favorable binding energy of −7.4 kcal mol⁻¹ closely matching the experimental value of −5.9 kcal mol⁻¹. The subsequent molecular dynamics simulation and MM-GBSA analysis recognized Asp98 as the most dominant residue, accounting for 40% of the total binding energy, established through a persistent hydrogen bonding with the histamine −NH₃⁺ group, the latter further held in place through the N–H∙∙∙O hydrogen bonding with Tyr250. Unlike earlier literature proposals, the important role of Thr190 is not evident in hydrogen bonds through its −OH group, but rather in the C–H∙∙∙π contacts with the imidazole ring, while its former moiety is constantly engaged in the hydrogen bonding with Asp186. Lastly, quantum-chemical calculations within the receptor cluster model and utilizing the empirical quantization of the ionizable X–H bonds (X = N, O, S), supported the deuteration-induced affinity increase, with the calculated difference in the binding free energy of −0.85 kcal mol⁻¹, being in excellent agreement with an experimental value of −0.75 kcal mol⁻¹, thus confirming the relevance of hydrogen bonding for the H₂ receptor activation.

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation

The computational prediction of relative binding free energies is a crucial goal for drug discovery, and G protein-coupled receptors (GPCRs) are arguably the most important drug target class. However, they present increased complexity to model compared to soluble globular proteins. Despite breakthroughs, experimental X-ray crystal and cryo-EM structures are challenging to attain, meaning computational models of the receptor and ligand binding mode are sometimes necessary. This leads to uncertainty in understanding ligand-protein binding induced changes such as, water positioning and displacement, side chain positioning, hydrogen bond networks, and the overall structure of the hydration shell around the ligand and protein. In other words, the very elements that define structure activity relationships (SARs) and are crucial for accurate binding free energy calculations are typically more uncertain for GPCRs. In this work we use free energy perturbation (FEP) to predict the relative binding free energies for ligands of two different GPCRs. We pinpoint the key aspects for success such as the important role of key water molecules, amino acid ionization states, and the benefit of equilibration with specific ligands. Initial calculations following typical FEP setup and execution protocols delivered no correlation with experiment, but we show how results are improved in a logical and systematic way. This approach gave, in the best cases, a coefficient of determination (R²) compared with experiment in the range of 0.6-0.9 and mean unsigned errors compared to experiment of 0.6-0.7 kcal/mol. We anticipate that our findings will be applicable to other difficult-to-model protein ligand data sets and be of wide interest for the community to continue improving FE binding energy predictions.

In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors

Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.

Insights into adenosine A2A receptor activation through cooperative modulation of agonist and allosteric lipid interactions

The activation process of G protein-coupled receptors (GPCRs) has been extensively studied, both experimentally and computationally. In particular, Molecular Dynamics (MD) simulations have proven useful in exploring GPCR conformational space. The typical behaviour of class A GPCRs, when subjected to unbiased MD simulations from their crystallized inactive state, is to fluctuate between inactive and intermediate(s) conformations, even with bound agonist. Fully active conformation(s) are rarely stabilized unless a G protein is also bound. Despite several crystal structures of the adenosine A2a receptor (A2aR) having been resolved in complex with co-crystallized agonists and Gs protein, its agonist-mediated activation process is still not completely understood. In order to thoroughly examine the conformational landscape of A2aR activation, we performed unbiased microsecond-length MD simulations in quadruplicate, starting from the inactive conformation either in apo or with bound agonists: endogenous adenosine or synthetic NECA, embedded in two homogeneous phospholipid membranes: 1,2-dioleoyl-sn-glycerol-3-phosphoglycerol (DOPG) or 1,2-dioleoyl-sn-glycerol-3-phosphocholine (DOPC). In DOPC with bound adenosine or NECA, we observe transition to an intermediate receptor conformation consistent with the known adenosine-bound crystal state. In apo state in DOPG, two different intermediate conformations are obtained. One is similar to that observed with bound adenosine in DOPC, while the other is closer to the active state but not yet fully active. Exclusively, in DOPG with bound adenosine or NECA, we reproducibly identify receptor conformations with fully active features, which are able to dock Gs protein. These different receptor conformations can be attributed to the action/absence of agonist and phospholipid-mediated allosteric effects on the intracellular side of the receptor.

Comparison of Orexin 1 and Orexin 2 Ligand Binding Modes Using X-ray Crystallography and Computational Analysis

The orexin system, which consists of the two G protein-coupled receptors OX1 and OX2, activated by the neuropeptides OX-A and OX-B, is firmly established as a key regulator of behavioral arousal, sleep, and wakefulness and has been an area of intense research effort over the past two decades. X-ray structures of the receptors in complex with 10 new antagonist ligands from diverse chemotypes are presented, which complement the existing structural information for the system and highlight the critical importance of lipophilic hotspots and water molecules for these peptidergic GPCR targets. Learnings from the structural information regarding the utility of pharmacophore models and how selectivity between OX1 and OX2 can be achieved are discussed.

Water envelope has a critical impact on the design of protein–protein interaction inhibitors

We show that a water envelope network plays a critical role in protein–protein interactions (PPI).

Surface coatings alter transcriptional responses to silver nanoparticles following oral exposure

Silver nanoparticles (AgNPs) are used in food packaging materials, dental care products and other consumer goods and can result in oral exposure. To determine whether AgNP coatings modulate transcriptional responses to AgNP exposure, we exposed mice orally to 20 nm citrate (cit)-coated AgNPs (cit-AgNPs) or polyvinylpyrrolidone (PVP)-coated AgNPs (PVP-AgNPs) at a 4 mg/kg dose for 7 consecutive days and analyzed changes in the expression of protein-coding genes and long noncoding RNAs (lncRNAs), a new class of regulatory RNAs, in the liver. We identified unique and common expression signatures of protein-coding and lncRNA genes, altered biological processes and signaling pathways, and coding-non-coding gene interactions for cit-AgNPs and PVP-AgNPs. Commonly regulated genes comprised only about 10 and 20 percent of all differentially expressed genes in PVP-AgNP and cit-AgNP exposed mice, respectively. Commonly regulated biological processes included glutathione metabolic process and cellular oxidant detoxification. Commonly regulated pathways included Keap-Nrf2, PPAR, MAPK and IL-6 signaling pathways. The coding-non-coding gene co-expression analysis revealed that protein-coding genes were co-expressed with a variable number of lncRNAs ranging from one to twenty three and may share functional roles with the protein-coding genes. PVP-AgNP exposure induced a more robust transcriptional response than cit-AgNP exposure characterized by more than two-fold higher number of differentially expressed both protein- coding and lncRNA genes. Our data demonstrate that the surface coating strongly modulates the spectrum and the number of differentially expressed genes after oral AgNP exposure. On the other hand, our data suggest that AgNP exposure can alter drug and chemical sensitivity, metabolic homeostasis and cancer risk irrespective of the coating type, warranting further investigations.

Massive in Silico Study of Noble Gas Binding to the Structural Proteome

Noble gases are chemically inert, and it was therefore thought they would have little effect on biology. Paradoxically, it was found that they do exhibit a wide range of biological effects, many of which are target-specific and potentially useful and some of which have been demonstrated in vivo. The underlying mechanisms by which useful pharmacology, such as tissue and neuroprotection, anti-addiction effects, and analgesia, is elicited are relatively unexplored. Experiments to probe the interactions of noble gases with specific proteins are more difficult with gases than those with other chemicals. It is clearly impractical to conduct the large number of gas-protein experiments required to gain a complete picture of noble gas biology. Given the simplicity of atoms as ligands, in silico methods provide an opportunity to gain insight into which noble gas-protein interactions are worthy of further experimental or advanced computational investigation. Our previous validation studies showed that in silico methods can accurately predict experimentally determined noble gas binding sites in X-ray structures of proteins. Here, we summarize the largest reported in silico reverse docking study involving 127 854 protein structures and the five nonradioactive noble gases. We describe how these computational screening methods are implemented, summarize the main types of interactions that occur between noble gases and target proteins, describe how the massive data set that this study generated can be analyzed (freely available at group18.csiro.au), and provide the NDMA receptor as an example of how these data can be used to understand the molecular pharmacology underlying the biology of the noble gases. We encourage chemical biologists to access the data and use them to expand the knowledge base of noble gas pharmacology, and to use this information, together with more efficient delivery systems, to develop "atomic drugs" that can fully exploit their considerable and relatively unexplored potential in medicine.

Chemokine Receptor Crystal Structures: What Can Be Learned from Them?

Chemokine receptors belong to the class A of G protein-coupled receptors (GPCRs) and are implicated in a wide variety of physiologic functions, mostly related to the homeostasis of the immune system. Chemokine receptors are also involved in multiple pathologic processes, including immune and autoimmune diseases, as well as cancer. Hence, several members of this GPCR subfamily are considered to be very relevant therapeutic targets. Since drug discovery efforts can be significantly reinforced by the availability of crystal structures, substantial efforts in the area of chemokine receptor structural biology could dramatically increase the outcome of drug discovery campaigns. This short review summarizes the available data on chemokine receptor crystal structures, discusses the numerous applications from chemokine receptor structures that can enhance the daily work of molecular pharmacologists, and describes the challenges and pitfalls to consider when relying on crystal structures for further research applications. SIGNIFICANCE STATEMENT: This short review summarizes the available data on chemokine receptor crystal structures, discusses the numerous applications from chemokine receptor structures that can enhance the daily work of molecular pharmacologists, and describes the challenges and pitfalls to consider when relying on crystal structures for further research applications.

The role of water molecules in phototransduction of retinal proteins and G protein-coupled receptors

G protein coupled receptors (GPCRs) are a key family of membrane proteins in all eukaryotes and also very important drug targets for medical intervention. The extensively studied visual pigment rhodopsin is a prime example of a family A GPCR. Its chromophore ligand retinal is covalently linked to a lysine in helix seven forming a protonated Schiff base. Interestingly, this is the same situation in other-non-GPCR-retinal proteins, like the prototype light-driven microbial proton pump bacteriorhodopsin, albeit there is no (or only a very remote) phylogenetical link. Close to the retinal ligand, several water molecules help to organise a functionally important hydrogen bond network that undergoes significant changes during photo-activation. Such water-mediated networks are also critical for ligand binding of other GPCRs and they are becoming increasingly important in drug discovery. GPCRs also contain a partially conserved water mediated hydrogen bond network that stabilises the ground state of the receptor, and rearrangement of this network leads to the stabilization of the active state. Some water molecules have a specific role in this process to appropriately orient specific residues relative to the Schiff base, and to modulate the fine structure of the transmembrane bundle, particularly near the intracellular G protein binding site. While the atomic details of these mechanisms are still missing, the recent developments in free electron lasers (FELs) are enabling us to begin to observe the changes in waters and relevant side chains shortly after photo activation at an unprecedented level of spatial and temporal resolution.

Understanding Ligand Binding Selectivity in a Prototypical GPCR Family

Adenosine receptors are involved in many pathological conditions and are thus promising drug targets. However, developing drugs that target this GPCR subfamily is a challenging task. A number of drug candidates fail due to lack of selectivity which results in unwanted side effects. The extensive structural similarity of adenosine receptors complicates the design of selective ligands. The problem of selective targeting is a general concern in GPCRs, and in this respect adenosine receptors are a prototypical example. Here we use enhanced sampling simulations to decipher the determinants of selectivity of ligands in A_2a and A₁ adenosine receptors. Our model shows how small differences in the binding pocket and in the water network around the ligand can be leveraged to achieve selectivity.

A piggyBac-based TANGO GFP assay for high throughput screening of GPCR ligands in live cells

BACKGROUND

GPCRs are considered essential for various physiological processes and have been the most productive drug targets. Therefore, development of the methods of GPCR ligands screening is a high priority for pharmaceutical industries and research institutions.

METHODS

We developed a potential method (piggyBac-TANGO) based on the TANGO and PRESTO-TANGO assays. The system was optimized with a piggyBac transposon as a transgene vehicle, and eGFP was used as a reporter instead of luciferase. The assay was validated in the HEK 293T and U87-MG cell lines and antagonist activities of the compounds were assessed. The transgene copy number and long-term stability were evaluated by qPCR. Then, we performed a DRD2-targeted screening for natural products using the piggyBac-TANGO assay.

RESULTS

The validation assay showed that using the piggyBac transposon as a transgene vehicle produced high signal-to-background ratio and stable readout confirmed by investigation of the transgene copy number and long-term stability. Use of eGFP instead of luciferase as a reporter enabled to create a high throughput system suitable for live cells. Moreover, the piggyBac-TANGO assay permitted versatile detection of antagonist activity of compounds and was not limited to a particular cell type. With the use of the piggyBac-TANGO assay, we have successfully identified a novel agonist of DRD2.

CONCLUSION

Thus, the results indicate that the piggyBac-TANGO method is a user-friendly, robust and imaging-based assay that provides a novel approach to high throughput GPCR-targeted ligand screening and drug development.

Structure-Function Relationships of Olfactory and Taste Receptors

The field of chemical senses has made major progress in understanding the cellular mechanisms of olfaction and taste in the past 2 decades. However, the molecular understanding of odor and taste recognition is still lagging far behind and will require solving multiple structures of the relevant full-length receptors in complex with native ligands to achieve this goal. However, the development of multiple complimentary strategies for the structure determination of G protein-coupled receptors (GPCRs) makes this goal realistic. The common conundrum of how multi-specific receptors that recognize a large number of different ligands results in a sensory perception in the brain will only be fully understood by a combination of high-resolution receptor structures and functional studies. This review discusses the first steps on this pathway, including biochemical and physiological assays, forward genetics approaches, molecular modeling, and the first steps towards the structural biology of olfactory and taste receptors.

Structure-based exploration and pharmacological evaluation of N-substituted piperidin-4-yl-methanamine CXCR4 chemokine receptor antagonists

Using the available structural information of the chemokine receptor CXCR4, we present hit finding and hit exploration studies that make use of virtual fragment screening, design, synthesis and structure-activity relationship (SAR) studies. Fragment 2 was identified as virtual screening hit and used as a starting point for the exploration of 31 N-substituted piperidin-4-yl-methanamine derivatives to investigate and improve the interactions with the CXCR4 binding site. Additionally, subtle structural ligand changes lead to distinct interactions with CXCR4 resulting in a full to partial displacement of CXCL12 binding and competitive and/or non-competitive antagonism. Three-dimensional quantitative structure-activity relationship (3D-QSAR) and binding model studies were used to identify important hydrophobic interactions that determine binding affinity and indicate key ligand-receptor interactions.

Could the presence of sodium ion influence the accuracy and precision of the ligand-posing in the human A2A adenosine receptor orthosteric binding site using a molecular docking approach? Insights from Dockbench

The allosteric modulation of G protein-coupled receptors (GPCRs) by sodium ions has received considerable attention as crystal structures of several receptors, in their inactive conformation, show a Na⁺ ion bound to specific residues which, in the human A_2A adenosine receptor (hA_2A AR), are Ser91^3.39, Trp246^6.48, Asn280^7.45, and Asn284^7.49. A cluster of water molecules completes the coordination of the sodium ion in the putative allosteric site. It is absolutely consolidated that the progress made in the field of GPCRs structural determination has increased the adoption of docking-driven approaches for the identification or the optimization of novel potent and selective ligands. Despite the extensive use of docking protocols in virtual screening approaches, to date, almost any of these studies have been carried out without taking into account the presence of the sodium cation and its first solvation shell in the putative allosteric binding site. In this study, we have focused our attention on determining how the presence of sodium ion binding and additionally its first hydration sphere, in hA_2AAR could influence the ligand positioning accuracy during molecular docking simulations for most of the available resting and activated hA_2A AR crystal structures, using DockBench as a comparative benchmarking tool and implementing a new correlation coefficient (EM). This work provides indications on the evidence that the posing performance (accuracy and/or precision) of the docking protocols in reproducing the crystallographic poses of different hA_2A AR antagonists is generally increased in the presence of the sodium cation and its first solvation shell, in agreement with experimental observations. Consequently, the inclusion of sodium ion and its first solvation shell should be considered in order to facilitate the selection of new potential ligands in all molecular docking-based virtual screening protocols that aim to find novel GPCRs antagonists and inverse agonists.

Psychedelics as anti-inflammatory agents

Serotonin (5-hydroxytryptamine, 5-HT)_2A receptor agonists have recently emerged as promising new treatment options for a variety of disorders. The recent success of these agonists, also known as psychedelics, like psilocybin for the treatment of anxiety, depression, obsessive-compulsive disorder (OCD), and addiction, has ushered in a renaissance in the way these compounds are perceived in the medical community and populace at large. One emerging therapeutic area that holds significant promise is their use as anti-inflammatory agents. Activation of 5-HT_2A receptors produces potent anti-inflammatory effects in animal models of human inflammatory disorders at sub-behavioural levels. This review discusses the role of the 5-HT_2A receptor in the inflammatory response, as well as highlight studies using the 5-HT_2A agonist (R)-2,5-dimethoxy-4-iodoamphetamine [(R)-DOI] to treat inflammation in cellular and animal models. It also examines potential mechanisms by which 5-HT_2A agonists produce their therapeutic effects. Overall, psychedelics regulate inflammatory pathways via novel mechanisms, and may represent a new and exciting treatment strategy for several inflammatory disorders.

Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands

Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.

Binding Affinity via Docking: Fact and Fiction

In 1982, Kuntz et al. published an article with the title “A Geometric Approach to Macromolecule-Ligand Interactions”, where they described a method “to explore geometrically feasible alignment of ligands and receptors of known structure”. Since then, small molecule docking has been employed as a fast way to estimate the binding pose of a given compound within a specific target protein and also to predict binding affinity. Remarkably, the first docking method suggested by Kuntz and colleagues aimed to predict binding poses but very little was specified about binding affinity. This raises the question as to whether docking is the right tool to estimate binding affinity. The short answer is no, and this has been concluded in several comprehensive analyses. However, in this opinion paper we discuss several critical aspects that need to be reconsidered before a reliable binding affinity prediction through docking is realistic. These are not the only issues that need to be considered, but they are perhaps the most critical ones. We also consider that in spite of the huge efforts to enhance scoring functions, the accuracy of binding affinity predictions is perhaps only as good as it was 10–20 years ago. There are several underlying reasons for this poor performance and these are analyzed. In particular, we focus on the role of the solvent (water), the poor description of H-bonding and the lack of the systems’ true dynamics. We hope to provide readers with potential insights and tools to overcome the challenging issues related to binding affinity prediction via docking.

Chemical Diversity in the G Protein-Coupled Receptor Superfamily

G protein-coupled receptors (GPCRs) are the largest family of cell signaling transmembrane proteins that can be modulated by a plethora of chemical compounds. Systematic cheminformatics analysis of structurally and pharmacologically characterized GPCR ligands shows that cocrystallized GPCR ligands cover a significant part of chemical ligand space, despite their limited number. Many GPCR ligands and substructures interact with multiple receptors, providing a basis for polypharmacological ligand design. Experimentally determined GPCR structures represent a variety of binding sites and receptor-ligand interactions that can be translated to chemically similar ligands for which structural data are lacking. This integration of structural, pharmacological, and chemical information on GPCR-ligand interactions enables the extension of the structural GPCR-ligand interactome and the structure-based design of novel modulators of GPCR function.

Quantitative metrics for drug-target ligandability

Ligandability is a prerequisite for druggability and is a much easier concept to understand, model and predict because it does not depend on the complex pharmacodynamic and pharmacokinetic mechanisms in the human body. In this review, we consider a metric for quantifying ligandability from experimental data. We discuss ligandability in terms of the balance between effort and reward. The metric is evaluated for a standard set of well-studied drug targets - some traditionally considered to be ligandable and some regarded as difficult. We suggest that this metric should be used to systematically improve computational predictions of ligandability, which can then be applied to novel drug targets to predict their tractability.

Approaches to target tractability assessment - a practical perspective

The assessment of the suitability of novel targets to intervention by different modalities, e.g. small molecules or antibodies, is increasingly seen as important in helping to select the most progressable targets at the outset of a drug discovery project. This perspective considers differing aspects of tractability and how it can be assessed using in silico and experimental approaches. We also share some of our experiences in using these approaches.

Methodologies for the Examination of Water in GPCRs

The following chapter examines some of the current "state-of-the-art" tools for predicting, scoring, and examining explicit water molecules in proteins and protein/ligand complexes, highlighting some of the ways information can be readily examined in a manner that is useful in a drug discovery process.

Development by Genetic Immunization of Monovalent Antibodies Against Human Vasoactive Intestinal Peptide Receptor 1 (VPAC1), New Innovative, and Versatile Tools to Study VPAC1 Receptor Function

Multi-membrane spanning proteins, such as G protein-coupled receptors (GPCRs) and ion channels, are extremely difficult to purify as native proteins. Consequently, the generation of antibodies that recognize the native conformation can be challenging. By combining genetic immunization, phage display, and biopanning, we identified a panel of monovalent antibodies (nanobodies) targeting the vasoactive intestinal peptide receptor 1 (VPAC1) receptor. The nine unique nanobodies that were classified into four different families based on their CDR3 amino acid sequence and length, were highly specific for the human receptor and bind VPAC1 with moderate affinity. They all recognize a similar epitope localized in the extracellular N-terminal domain of the receptor and distinct from the orthosteric binding site. In agreement with binding studies, which showed that the nanobodies did not interfere with VIP binding, all nanobodies were devoid of any functional properties. However, we observed that the binding of two nanobodies was slightly increased in the presence of VPAC1 agonists [vasoactive intestinal polypeptide (VIP) and pituitary adenylate cyclase-activating polypeptide-27 (PACAP-27)], but decreased in the presence of VPAC1 antagonist. As no evidence of allosteric activity was seen in VIP binding studies nor in functional assays, it is, therefore, possible that the two nanobodies may behave as very weak allosteric modulators of VPAC1, detectable only in some sensitive settings, but not in others. We demonstrated that the fluorescently labeled nanobodies detect VPAC1 on the surface of human leukocytes as efficiently as a reference mouse monoclonal antibody. We also developed a protocol allowing efficient detection of VPAC1 by immunohistochemistry in paraffin-embedded human gastrointestinal tissue sections. Thus, these nanobodies constitute new original tools to further investigate the role of VPAC1 in physiological and pathological conditions.

Methods for Virtual Screening of GPCR Targets: Approaches and Challenges

Virtual screening (VS) has become an integral part of the drug discovery process and is a valuable tool for finding novel chemical starting points for GPCR targets. Ligand-based VS makes use of biochemical data for known, active compounds and has been applied successfully to many diverse GPCRs. Recent progress in GPCR X-ray crystallography has made it possible to incorporate detailed structural information into the VS process. This chapter outlines the latest VS techniques along with examples that highlight successful applications of these methods. Best practices for increasing the likelihood of VS success, as well as ongoing challenges, are also discussed.

Molecular modeling study on resistance of WT/D473H SMO to antagonists LDE-225 and LEQ-506

The smoothened (SMO) receptor, an essential signal transducer in the Hedgehog pathway, was targeted with antagonists to suppress the tumor. It is interesting that SMO D473H mutation confers resistance on inhibitor LDE-225 rather than LEQ-506. In this paper, the binding modes of them against the wild type and mutant SMO receptors were identified to gain insights into the resistant and non-resistant factors, based on a comprehensive protocol involving molecular docking, molecular dynamic simulations, free energy calculation and decomposition. A comparison of resistant LDE-225 and non-resistant LEQ-506 indicates that the volume of the binding cavity decreases seriously in the mutant complex with resistant LDE-225. In addition, the D473H mutation disrupts the hydrogen bond network with residues R400 and Q477, which results in the TM6 conformation inward. Owing to the absence of the hydrogen bond, residues R400 and Q477 make weak contributions to LDE-225. However, the D473H mutation along with TM6 conformational change has no effect on non-resistant LEQ-506. Finally, the resistance ascribes to adverse interaction between the greater polarity of mutant residue H473 and the nonpolar phenmethyl of LDE-225. The elaborate insights into structural and energetic mechanism of drug resistance provide an effective strategy to design rationally non-resistant antagonists.

Enhancing Structure Prediction and Design of Soluble and Membrane Proteins with Explicit Solvent-Protein Interactions

Solvent molecules interact intimately with proteins and can profoundly regulate their structure and function. However, accurately and efficiently modeling protein solvation effects at the molecular level has been challenging. Here, we present a method that improves the atomic-level modeling of soluble and membrane protein structures and binding by efficiently predicting de novo protein-solvent molecule interactions. The method predicted with unprecedented accuracy buried water molecule positions, solvated protein conformations, and challenging mutational effects on protein binding. When applied to homology modeling, solvent-bound membrane protein structures, pockets, and cavities were recapitulated with near-atomic precision even from distant homologs. Blindly refined atomic-level structures of evolutionary distant G protein-coupled receptors imply strikingly different functional roles of buried solvent between receptor classes. The method should prove useful for refining low-resolution protein structures, accurately modeling drug-binding sites in structurally uncharacterized receptors, and designing solvent-mediated protein catalysis, recognition, ligand binding, and membrane protein signaling.

Discovery of new GPCR ligands to illuminate new biology

Although a plurality of drugs target G-protein-coupled receptors (GPCRs), most have emerged from classical medicinal chemistry and pharmacology programs and resemble one another structurally and functionally. Though effective, these drugs are often promiscuous. With the realization that GPCRs signal via multiple pathways, and with the emergence of crystal structures for this family of proteins, there is an opportunity to target GPCRs with new chemotypes and confer new signaling modalities. We consider structure-based and physical screening methods that have led to the discovery of new reagents, focusing particularly on the former. We illustrate their use against previously untargeted or orphan GPCRs, against allosteric sites, and against classical orthosteric sites that selectively activate one downstream pathway over others. The ligands that emerge are often chemically novel, which can lead to new biological effects.

How Ligands Illuminate GPCR Molecular Pharmacology

G protein-coupled receptors (GPCRs), which are modulated by a variety of endogenous and synthetic ligands, represent the largest family of druggable targets in the human genome. Recent structural and molecular studies have both transformed and expanded classical concepts of receptor pharmacology and have begun to illuminate the distinct mechanisms by which structurally, chemically, and functionally diverse ligands modulate GPCR function. These molecular insights into ligand engagement and action have enabled new computational methods and accelerated the discovery of novel ligands and tool compounds, especially for understudied and orphan GPCRs. These advances promise to streamline the development of GPCR-targeted medications.

Signalling assemblies: the odds of symmetry

The assembly of proteins into complexes is fundamental to nearly all biological signalling processes. Symmetry is a dominant feature of the structures of experimentally determined protein complexes, observed in the vast majority of homomers and many heteromers. However, some asymmetric structures exist, and asymmetry also often forms transiently, intractable to traditional structure determination methods. Here, we explore the role of protein complex symmetry and asymmetry in cellular signalling, focusing on receptors, transcription factors and transmembrane channels, among other signalling assemblies. We highlight a recurrent tendency for asymmetry to be crucial for signalling function, often being associated with activated states. We conclude with a discussion of how consideration of protein complex symmetry and asymmetry has significant potential implications and applications for pharmacology and human disease.

Paralog-divergent Features May Help Reduce Off-target Effects of Drugs: Hints from Glucagon Subfamily Analysis

Side effects from targeted drugs remain a serious concern. One reason is the nonselective binding of a drug to unintended proteins such as its paralogs, which are highly homologous in sequences and have similar structures and drug-binding pockets. To identify targetable differences between paralogs, we analyzed two types (type-I and type-II) of functional divergence between two paralogs in the known target protein receptor family G-protein coupled receptors (GPCRs) at the amino acid level. Paralogous protein receptors in glucagon-like subfamily, glucagon receptor (GCGR) and glucagon-like peptide-1 receptor (GLP-1R), exhibit divergence in ligands and are clinically validated drug targets for type 2 diabetes. Our data showed that type-II amino acids were significantly enriched in the binding sites of antagonist MK-0893 to GCGR, which had a radical shift in physicochemical properties between GCGR and GLP-1R. We also examined the role of type-I amino acids between GCGR and GLP-1R. The divergent features between GCGR and GLP-1R paralogs may be helpful in their discrimination, thus enabling the identification of binding sites to reduce undesirable side effects and increase the target specificity of drugs.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Regulation of G Protein-Coupled Receptors by Ubiquitination

G protein-coupled receptors (GPCRs) comprise the largest family of membrane receptors that control many cellular processes and consequently often serve as drug targets. These receptors undergo a strict regulation by mechanisms such as internalization and desensitization, which are strongly influenced by posttranslational modifications. Ubiquitination is a posttranslational modification with a broad range of functions that is currently gaining increased appreciation as a regulator of GPCR activity. The role of ubiquitination in directing GPCRs for lysosomal degradation has already been well-established. Furthermore, this modification can also play a role in targeting membrane and endoplasmic reticulum-associated receptors to the proteasome. Most recently, ubiquitination was also shown to be involved in GPCR signaling. In this review, we present current knowledge on the molecular basis of GPCR regulation by ubiquitination, and highlight the importance of E3 ubiquitin ligases, deubiquitinating enzymes and β-arrestins. Finally, we discuss classical and newly-discovered functions of ubiquitination in controlling GPCR activity.

Structure-Based Discovery of New Antagonist and Biased Agonist Chemotypes for the Kappa Opioid Receptor

The ongoing epidemics of opioid overdose raises an urgent need for effective antiaddiction therapies and addiction-free painkillers. The κ-opioid receptor (KOR) has emerged as a promising target for both indications, raising demand for new chemotypes of KOR antagonists as well as G-protein-biased agonists. We employed the crystal structure of the KOR-JDTic complex and ligand-optimized structural templates to perform virtual screening of available compound libraries for new KOR ligands. The prospective virtual screening campaign yielded a high 32% hit rate, identifying novel fragment-like and lead-like chemotypes of KOR ligands. A round of optimization resulted in 11 new submicromolar KOR binders (best K_i = 90 nM). Functional assessment confirmed at least two compounds as potent KOR antagonists, while compound 81 was identified as a potent G_i biased agonist for KOR with minimal β-arrestin recruitment. These results support virtual screening as an effective tool for discovery of new lead chemotypes with therapeutically relevant functional profiles.

Rearrangement of a polar core provides a conserved mechanism for constitutive activation of class B G protein-coupled receptors

The glucagon receptor (GCGR) belongs to the secretin-like (class B) family of G protein-coupled receptors (GPCRs) and is activated by the peptide hormone glucagon. The structures of an activated class B GPCR have remained unsolved, preventing a mechanistic understanding of how these receptors are activated. Using a combination of structural modeling and mutagenesis studies, we present here two modes of ligand-independent activation of GCGR. First, we identified a GCGR-specific hydrophobic lock comprising Met-338 and Phe-345 within the IC3 loop and transmembrane helix 6 (TM6) and found that this lock stabilizes the TM6 helix in the inactive conformation. Disruption of this hydrophobic lock led to constitutive G protein and arrestin signaling. Second, we discovered a polar core comprising conserved residues in TM2, TM3, TM6, and TM7, and mutations that disrupt this polar core led to constitutive GCGR activity. On the basis of these results, we propose a mechanistic model of GCGR activation in which TM6 is held in an inactive conformation by the conserved polar core and the hydrophobic lock. Mutations that disrupt these inhibitory elements allow TM6 to swing outward to adopt an active TM6 conformation similar to that of the canonical β₂-adrenergic receptor complexed with G protein and to that of rhodopsin complexed with arrestin. Importantly, mutations in the corresponding polar core of several other members of class B GPCRs, including PTH1R, PAC1R, VIP1R, and CRFR1, also induce constitutive G protein signaling, suggesting that the rearrangement of the polar core is a conserved mechanism for class B GPCR activation.

In silico design of novel probes for the atypical opioid receptor MRGPRX2

The primate-exclusive MRGPRX2 G protein-coupled receptor (GPCR) has been suggested to modulate pain and itch. Despite putative peptide and small molecule MRGPRX2 agonists, selective nanomolar potency probes have not yet been reported. To identify a MRGPRX2 probe, we first screened 5,695 small molecules and found many opioid compounds activated MRGPRX2, including (−)- and (+)-morphine, hydrocodone, sinomenine, dextromethorphan and the prodynorphin-derived peptides, dynorphin A, dynorphin B, and α- and β-neoendorphin. We used these to select for mutagenesis-validated homology models and docked almost 4 million small molecules. From this docking, we predicted ZINC-3573, which represents a potent MRGPRX2-selective agonist, showing little activity against 315 other GPCRs and 97 representative kinases, and an essentially inactive enantiomer. ZINC-3573 activates endogenous MRGPRX2 in a human mast cell line inducing degranulation and calcium release. MRGPRX2 is a unique atypical opioid-like receptor important for modulating mast cell degranulation, which can now be specifically modulated with ZINC-3573.

Structural Analysis of Chemokine Receptor-Ligand Interactions

This review focuses on the construction and application of structural chemokine receptor models for the elucidation of molecular determinants of chemokine receptor modulation and the structure-based discovery and design of chemokine receptor ligands. A comparative analysis of ligand binding pockets in chemokine receptors is presented, including a detailed description of the CXCR4, CCR2, CCR5, CCR9, and US28 X-ray structures, and their implication for modeling molecular interactions of chemokine receptors with small-molecule ligands, peptide ligands, and large antibodies and chemokines. These studies demonstrate how the integration of new structural information on chemokine receptors with extensive structure–activity relationship and site-directed mutagenesis data facilitates the prediction of the structure of chemokine receptor–ligand complexes that have not been crystallized. Finally, a review of structure-based ligand discovery and design studies based on chemokine receptor crystal structures and homology models illustrates the possibilities and challenges to find novel ligands for chemokine receptors.

Recent Advances and Applications of Molecular Docking to G Protein-Coupled Receptors

The growing number of studies on G protein-coupled receptors (GPCRs) family are a source of noticeable improvement in our understanding of the functioning of these proteins. GPCRs are responsible for a vast part of signaling in vertebrates and, as such, invariably remain in the spotlight of medicinal chemistry. A deeper insight into the underlying mechanisms of interesting phenomena observed in GPCRs, such as biased signaling or allosteric modulation, can be gained with experimental and computational studies. The latter play an important role in this process, since they allow for observations on scales inaccessible for most other methods. One of the key steps in such studies is proper computational reconstruction of actual ligand-receptor or protein-protein interactions, a process called molecular docking. A number of improvements and innovative applications of this method were documented recently. In this review, we focus particularly on innovations in docking to GPCRs.