A Surprising Recipe for Designing Biased Ligands.

The determination of the potential value of receptor trafficking at melanocortin receptors has been hampered by the absence of known biased ligands. Heterobivalent MC4R ligands linking agonist to antagonist small peptidic moieties were designed and found to act as Gαs enhancers while minimally activating ß-arrestin recruitment. The strategy invoked offers intriguing promise as a surprising approach that is possibly broadly applicable to the challenge of designing biased ligands at other GPCRs.

Current and Future Challenges in GPCR Drug Discovery.

GPCRs play a pervasive physiological role and, in turn, are the leading target class for pharmaceuticals. Beginning with the determination of the structure of rhodopsin, and dramatically accelerating since the reporting of the first ligand-mediated GPCR X-ray structures, our understanding of the structural and functional characteristics of these proteins has grown dramatically. Deploying this now rapidly emerging information for drug discovery has already been extensively demonstrated through a watershed of studies appearing in numerous scientific reports. Included in these expositions are areas such as sites and characteristics of ligand to GPCR binding, protein activation, effector bias, allosteric mechanisms, dimerization, polypharmacology and others. Computational chemistry studies are demonstrating an increasing role in capitalizing on the structural studies to further advance our understanding of these proteins as well as to drive drug discovery. Such drug discovery activities range from the design of orthosteric site inhibitors through, for example, allosteric modulators, biased ligands, partial agonists and bitopic ligands. Herein, these topics are outlined through specific examples in the hopes of providing a glimpse of the state of the field.

Latch and trigger role for R445 in DAT transport explains molecular basis of DTDS.

A recent study reports on five different mutations as sources of dopamine transporter (DAT) deficiency syndrome (DTDS). One of these mutations, R445C, is believed to be located on the intracellular side of DAT distal to the primary (S1) or secondary (S2) sites to which substrate binding is understood to occur. Thus, the molecular mechanism by which the R445C mutation results in DAT transport deficiency has eluded explanation. However, the recently reported X-ray structures of the endogenous amine transporters for dDAT and hSERT revealed the presence of a putative salt bridge between R445 and E428 suggesting a possible mechanism. To evaluate whether the R445C effect is a result of a salt bridge interaction, the mutants R445E, E428R, and the double mutant E428R/R445E were generated. The single mutants R445E and E428R displayed loss of binding and transport properties of the substrate [3H]DA and inhibitor [3H]CFT at the cell surface while the double mutant E428R/R445E, although nonfunctional, restored [3H]DA and [3H]CFT binding affinity to that of WT. Structure based analyses of these results led to a model wherein R445 plays a dual role in normal DAT function. R445 acts as a component of a latch in its formation of a salt bridge with E428 which holds the primary substrate binding site (S1) in place and helps enforce the inward closed protein state. When this salt bridge is broken, R445 acts as a trigger which disrupts a local polar network and leads to the release of the N-terminus from its position inducing the inward closed state to one allowing the inward open state. In this manner, both the loss of binding and transport properties of the R445C variant are explained.

X-ray structure based evaluation of analogs of citalopram: Compounds with increased affinity and selectivity compared with R-citalopram for the allosteric site (S2) on hSERT.

The recent publication of X-ray structures of SERT includes structures with the potent antidepressant S-Citalopram (S-Cit). Earlier predictions of ligand binding at both a primary (S1) and an allosteric modulator site (S2), were confirmed. We provide herein examples of a series of Citalopram analogs, showing distinct structure-activity relationship (SAR) at both sites that is independent of the SAR at the other site. Analogs with a higher affinity and selectivity than benchmark R-Citalopram (R-Cit) for the S2 versus the S1 site were identified. We deploy structural and computational analyses to explain this SAR and demonstrate the potential utility of the newly emerging X-ray structures within the neurotransmitter:sodium Symporter family for drug design.

Exploration of insights, opportunities and caveats provided by the X-ray structures of hSERT.

The recently reported X-ray structures of the human serotonin (5-HT) transporter SERT with bound inhibitors open new opportunities for drug discovery at SERT, selectivity design with respect to other neurotransmitter sodium transporters, and enhanced understanding of the molecular events involved in SERT action. Through computational and structural analyses, we explore the binding and migration of 5-HT at SERT. Consistent with earlier studies of leucine migration at the bacterial homolog of SERT, LeuT, we find multiple potential 'stopover' sites for 5-HT binding at SERT including the two (transmembrane S1 and extracellular vestibule S2) seen in the binding of the SSRI (S)-citalopram (S-Cit) to SERT, as well as other sites. Docking studies reveal the possibility of both hetero- (S-Cit+5-HT) and homo-dimeric (5-HT+5-HT) co-binding at both these sites which may explain earlier published allosteric activity observations and provide novel design strategies. Comparisons with substrate bound X-ray structures of the dopamine transporter reveal a number of potential sources of selectivity, some of which may be 'artificial' including target based, species related, experimental design related, and ligand dependent examples including substrate versus inhibitor related features.

7TM X-ray structures for class C GPCRs as new drug-discovery tools. 1. mGluR5.

We illustrate, with a focus on mGluR5, how the recently published, first X-ray structures of mGluR 7TM domains, specifically those of mGluR1 and mGluR5 complexed with negative allosteric modulators (NAMs), will begin to influence ligand- (e.g., drug- or sweetener-) discovery efforts involving class C GPCRs. With an extensive docking study allowing full ligand flexibility and full side chain flexibility of all residues in the ligand-binding cavity, we have predicted and analyzed the binding modes of a variety of structurally diverse mGluR5 NAM ligands, showing how the X-ray structures serve to effectively rationalize each ligand's binding characteristics. We demonstrated that the features that are inherent in our earlier overlay model are preserved in the protein structure-based docking models. We identified structurally diverse compounds, which potentially act as mGluR NAMs, and revealed binding-site differences by performing high-throughput docking using a database of approximately six million structures of commercially available compounds and the mGluR1 and mGluR5 X-ray structures. By comparing the 7TM domains of the mGluR5 and mGluR1 X-rays structures, we identified selectivity factors within group I of the mGluRs. Similarly, using homology models that we built for mGluR2 and mGluR4, we have identified the factors leading to the selectivity between group I and groups II and III for ligands occupying the deepest portion of the mGluR5 binding cavity. Finally, we have proposed a structure-based explanation of the pharmacological switching within a set of positive allosteric modulators (PAMs) and their corresponding, very close NAM analogs.

The role of experimental and computational structural approaches in 7TM drug discovery.

INTRODUCTION: Starting with the published X-ray structures of ligand-mediated 7TM proteins in 2007, experimental approaches, led by X-ray structure determinations, and computational approaches, led by docking and molecular dynamics, have converged to elaborate our understanding of this field and demonstrate their effectiveness in drug discovery. AREAS COVERED: The authors review the structural information that has emerged for ligand-mediated 7TM proteins, including the class A, B, C, and F receptors, focusing on the 7TM domains for the multi-domain proteins. The authors describe the key regions associated with ligand binding as well as features responsible for function such as activation versus inhibition and biased signaling. Furthermore, the authors summarize the effectiveness of computational studies to help clarify the structure-function information and their use for drug discovery. EXPERT OPINION: There is now a significant amount of structural information covering a range of 7TM protein classes (A, B, C, and F) and activation states. For these and closely related proteins, structure-based drug discovery has proven to be a powerful tool. More structural information is needed with respect to dimerization, 7TM proteins with ß-arrestin to help in understanding the control of biased signaling, and full-protein structure determinations for non-class A proteins to help in understanding and controlling their functioning. Finally, the use of the existing structural information to target new sites on these proteins needs further exploration.

X-ray structural information of GPCRs in drug design: what are the limitations and where do we go?

INTRODUCTION: In 2007, the X-ray structural determination of non-rhodopsin G-Protein coupled receptors (GPCRs), considered the most extensively targeted protein class for marketed drugs, commenced. With the relatively rapid availability of additional structures, an assessment of the progression made is needed in addition to the assessment of the understandings gleaned, deployment successes and forthcoming prospects. AREAS COVERED: The author reviews the approaches and tools that have made it possible to determine the three dimensional structures of GPCRs using X-ray crystallography. Furthermore, the author describes the methods suited for crystallization of membrane bound GPCR proteins including the lipidic cubic phase and various protein modification approaches. The author also provides highlights, from the literature, of the structures determined to date including targets solved, the nature of the content provided (such as selectivity, activating vs. inactivating determinants) and how these structural features relate to drug design strategies. EXPERT OPINION: The GPCR X-ray structures that have been so far determined have yielded significant information. This has presented dramatic evidence concerning their ability to impact the discovery of compounds through their action as traditional, orthosteric modulators. It is, however, noted that more challenging design strategies, such as identifying biased agonists and the use of sites remote from the orthosteric site for allosteric modulation, are still in their infancy.

Exploration of structure-based drug design opportunities for mGluRs.

The metabotropic glutamate receptors (mGluRs) are a subset of the Class C G-Protein Coupled Receptors (GPCRs). Recently, an emerging strategy for drug-discovery efforts targeting mGluRs has been to develop compounds acting at the so-called allosteric site in the 7-transmembrane (7TM) domain, common to all GPCRs, rather than the extracellular (EC) domain containing the orthosteric glutamate-binding site. We examine herein some of the intrinsic relative merits of targeting these two domains. Comparisons are made among amino-acid sequences in the two domains and among X-ray structures and homology models of the EC domain. We show that there is greater sequence diversity in the EC domains than in the transmembrane (TM) domains. Thus, contrary to generally accepted descriptions of there being greater evolutionary pressure to preserve the EC domain, it is the 7TM domain that is more highly conserved. Within the EC domain, the glutamate-binding site of the Venus flytrap region has hitherto received the most attention as a target site. Analysis of examples of the three-dimensional structures of the EC domains at the glutamate-binding site reveals differences as well, thereby supporting the viability of targeting the EC domain, even at the glutamate-binding site, for drug discovery. To exemplify this strategy, we present examples of active compounds identified via high-throughput docking in the EC region.

Electrostatic chiral distinction: tetrahedral model dimers.

Although chiral distinction plays a pervasive role in chemistry, a complete understanding of how this takes place is still lacking. In this work, we expand the earlier described minimal requirement of so called four-point interactions (vide infra). We focus on chiral point charge model systems as a means to aid in the dissection of the underlying, operative principles. We also construct models with defined symmetry characteristics. By considering extensive constellations of diastereomeric complexes, we are able to identify emerging principles for chiral distinction. As previously postulated, all the diastereomeric complexes, regardless of their nominal contact-points, possess a chiral distinction energy. In the comparison of complexes, we find that, contrary to chemical intuition, the magnitude of chiral distinction does not correlate with the stability of the complexes, i.e., consideration of low energy complexes may not be an effective way to evaluate chiral distinction. Similarly, we do not find a correlation between the number of contact-points and chiral distinction. Moreover, favorable interactions and facile chiral distinction appear to be unrelated. We also see some tendency for greater chiral distinction in less symmetric systems, although this may not be general. These studies can now form the basis to fold in higher levels of complexity into the models so as to gain further insights into the nature of chiral distinction.

X-ray structure breakthroughs in the GPCR transmembrane region.

G-protein-coupled receptor (GPCR) proteins [Lundstrom KH, Chiu ML, editors. G protein-coupled receptors in drug discovery. CRC Press; 2006] are the single largest drug target, representing 25-50% of marketed drugs [Overington JP, Al-Lazikani B, Hopkins AL. How many drug targets are there? Nat Rev Drug Discov 2006;5(12):993-6; Parrill AL. Crystal structures of a second G protein-coupled receptor: triumphs and implications. ChemMedChem 2008;3:1021-3]. While there are six subclasses of GPCR proteins, the hallmark of all GPCR proteins is the transmembrane-spanning region. The general architecture of this transmembrane (TM) region has been known for some time to contain seven alpha-helices. From a drug discovery and design perspective, structural information of the GPCRs has been sought as a tool for structure-based drug design. The advances in the past decade of technologies for structure-based design have proven to be useful in a number of areas. Invoking these approaches for GPCR targets has remained challenging. Until recently, the most closely related structures available for GPCR modeling have been those of bovine rhodopsin. While a representative of class A GPCRs, bovine rhodopsin is not a ligand-activated GPCR and is fairly distant in sequence homology to other class A GPCRs. Thus, there is a variable degree of uncertainty in the use of the rhodopsin X-ray structure as a template for homology modeling of other GPCR targets. Recent publications of X-ray structures of class A GPCRs now offer the opportunity to better understand the molecular mechanism of action at the atomic level, to deploy X-ray structures directly for their use in structure-based design, and to provide more promising templates for many other ligand-mediated GPCRs. We summarize herein some of the recent findings in this area and provide an initial perspective of the emerging opportunities, possible limitations, and remaining questions. Other aspects of the recent X-ray structures are described by Weis and Kobilka [Weis WI, Kobilka BK. Structural insights into G-protein-coupled receptor activation. Curr Opin Struct Biol 2008;18:734-40] and Mustafi and Palczewski [Mustafi D, Palczewski K. Topology of class A G protein-coupled receptors: insights gained from crystal structures of rhodopsins, adrenergic and adenosine receptors. Mol Pharmacol 2009;75:1-12].

Driving forces for ligand migration in the leucine transporter.

The potential mechanisms for driving a ligand to and through a transporter are examined for the leucine transporter using the recently published X-ray structure. Through analyses with computational methods, including investigation of electrostatic properties, site and channel identification, and docking studies, a picture emerges whereby dipolar patches which characterize the electrostatic field serve as the primary driving force to attract the ligand to the protein and begin propagation into it. The electrostatic forces are then augmented by hydrophobic forces in the transport stages, with conformational changes in the protein helping to accommodate the migration. We identify 12 sites that might be involved in ligand recognition and migration. One of these sites corresponds to the tricyclic antidepressant binding site observed in the recently published X-ray structures.

Use of the X-ray structure of the beta2-adrenergic receptor for drug discovery. Part 2: Identification of active compounds.

The recently published X-ray structures of the beta(2)-adrenergic receptor are the first examples of ligand-mediated GPCR crystal structures. We have previously performed computational studies that examine the potential viability of these structures for use in drug design, exploiting known ligand activities. Our previous study and a newly reported beta(2)/Timolol X-ray complex provide validation of the computational approaches. In the present work, we use the X-ray structures to extract, via in silico high-throughput docking, compounds from proprietary and commercial databases and demonstrate the successful identification of active compounds by radioligand binding.

Use of the X-ray structure of the Beta2-adrenergic receptor for drug discovery.

The recently reported X-ray structure of the Beta2-adrenergic receptor, the first reported crystal structure of a ligand-mediated GPCR, is used to explore its utility in computer-aided drug design. Validations were conducted with known beta blockers. This was followed by high-throughput docking studies with proprietary and commercial databases to further validate the X-ray structure's usefulness as a design tool and to explore the potential for discovery of novel chemical classes acting as Beta2 inhibitors. Our results include the finding of ligands with traditional beta-blocker motifs as well as new motifs, thereby serving to both validate the approach and project its usefulness in the finding and design of novel compounds.

Homology modeling of the serotonin transporter: insights into the primary escitalopram-binding site.

The serotonin transporter (SERT) is one of the neurotransmitter transporters that plays a critical role in the regulation of endogenous amine concentrations and therefore is an important target for therapeutic agents affecting the central nervous system. The recently published, high resolution X-ray structure of the closely related amino acid transporter, Aquifex aeolicus leucine transporter (LeuT), provides an opportunity to develop a three-dimensional model of the structure of SERT. We present herein a homology model of SERT using LeuT as the template and containing escitalopram as a bound ligand. Our model explains selectivities known from mutational studies and varying ligand data, which are discussed and illustrated in the paper.

Two contact-point chiral distinction: model CHFClBr dimers.

Dimers of the simple chiral molecule CHFClBr have been studied using a variety of computational approaches, including HF, MP2, and DFT B3LYP and the 6-31G*, 6-31++G**, and 6-311++G** basis sets. Both heterochiral and homochiral dimers were studied to allow analysis of the chiral distinction in these systems. The dimers were arranged in edge-to-edge orientations with assorted combinations of two contact-points ("2:2e") between the dimers. The monomers were constrained to tetrahedral symmetry. We demonstrate that chiral distinction does indeed occur in these two contact-point models. While the stabilization energies are driven by the interactions of the nearest atoms (contacts) in the complexes, the degree of chiral distinction is driven by the profile of changing atoms, which, in the present systems, are often the distal atoms of the complexes. Moreover, the chiral distinction does not correlate with the stabilization energies. The terms contact-points and interactions are defined.

Threading with chemostructural restrictions method for predicting fold and functionally significant residues: application to dipeptidylpeptidase IV (DPP-IV).

We present a new method for more accurate modeling of protein structure, called threading with chemostructural restrictions. This method addresses those cases in which a target sequence has only remote homologues of known structure for which sequence comparison methods cannot provide accurate alignments. Although remote homologues cannot provide an accurate model for the whole chain, they can be used in constructing practically useful models for the most conserved-and often the most interesting-part of the structure. For many proteins of interest, one can suggest certain chemostructural patterns for the native structure based on the available information on the structural superfamily of the protein, the type of activity, the sequence location of the functionally significant residues, and other factors. We use such patterns to restrict (1) a number of possible templates, and (2) a number of allowed chain conformations on a template. The latter restrictions are imposed in the form of additional template potentials (including terms acting as sequence anchors) that act on certain residues. This approach is tested on remote homologues of alpha/beta-hydrolases that have significant structural similarity in the positions of their catalytic triads. The study shows that, in spite of significant deviations between the model and the native structures, the surroundings of the catalytic triad (positions of C(alpha) atoms of 20-30 nearby residues) can be reproduced with accuracy of 2-3 A. We then apply the approach to predict the structure of dipeptidylpeptidase IV (DPP-IV). Using experimentally available data identifying the catalytic triad residues of DPP-IV (David et al., J Biol Chem 1993;268:17247-17252); we predict a model structure of the catalytic domain of DPP-IV based on the 3D fold of prolyl oligopeptidase (Fulop et al., Cell 1998;94:161-170) and use this structure for modeling the interaction of DPP-IV with inhibitor.

Two-point contact chiral distinction--a theoretical appraisal.

Ab initio calculations reveal chiral distinction in two-point contact CHFCIBr dimers, with chiral distinction energy of 1.5 kJ mol-1 between the SR and SS dimers fully optimized at the MP2/6-311++G** level.