Modulating G-protein coupled receptor/G-protein signal transduction by small molecules suggested by virtual screening

The Retinoid and Non-Retinoid Ligands of the Rod Visual G Protein-Coupled Receptor

G protein-coupled receptors (GPCRs) play a predominant role in the drug discovery effort. These cell surface receptors are activated by a variety of specific ligands that bind to the orthosteric binding pocket located in the extracellular part of the receptor. In addition, the potential binding sites located on the surface of the receptor enable their allosteric modulation with critical consequences for their function and pharmacology. For decades, drug discovery focused on targeting the GPCR orthosteric binding sites. However, finding that GPCRs can be modulated allosterically opened a new venue for developing novel pharmacological modulators with higher specificity. Alternatively, focus on discovering of non-retinoid small molecules beneficial in retinopathies associated with mutations in rhodopsin is currently a fast-growing pharmacological field. In this review, we summarize the accumulated knowledge on retinoid ligands and non-retinoid modulators of the light-sensing GPCR, rhodopsin and their potential in combating the specific vision-related pathologies. Also, recent findings reporting the potential of biologically active compounds derived from natural products as potent rod opsin modulators with beneficial effects against degenerative diseases related to this receptor are highlighted here.

Tauroursodeoxycholic acid binds to the G-protein site on light activated rhodopsin

The heterotrimeric G-protein binding site on G-protein coupled receptors remains relatively unexplored regarding its potential as a new target of therapeutic intervention or as a secondary site of action by the existing drugs. Tauroursodeoxycholic acid bears structural resemblance to several compounds that were previously identified to specifically bind to the light-activated form of the visual receptor rhodopsin and to inhibit its activation of transducin. We show that TUDCA stabilizes the active form of rhodopsin, metarhodopsin II, and does not display the detergent-like effects of common amphiphilic compounds that share the cholesterol scaffold structure, such as deoxycholic acid. Computer docking of TUDCA to the model of light-activated rhodopsin revealed that it interacts using similar mode of binding to the C-terminal domain of transducin alpha subunit. The ring regions of TUDCA made hydrophobic contacts with loop 3 region of rhodopsin, while the tail of TUDCA is exposed to solvent. The results show that TUDCA interacts specifically with rhodopsin, which may contribute to its wide-ranging effects on retina physiology and as a potential therapeutic compound for retina degenerative diseases.

Intracellular Receptor Modulation: Novel Approach to Target GPCRs

Recent crystal structures of multiple G protein-coupled receptors (GPCRs) have revealed a highly conserved intracellular pocket that can be used to modulate these receptors from the inside. This novel intracellular site partially overlaps with the G protein and β-arrestin binding site, providing a new manner of pharmacological intervention. Here we provide an update of the architecture and function of the intracellular region of GPCRs, until now portrayed as the signaling domain. We review the available evidence on the presence of intracellular binding sites among chemokine receptors and other class A GPCRs, as well as different strategies to target it, including small molecules, pepducins, and nanobodies. Finally, the potential advantages of intracellular (allosteric) ligands over orthosteric ligands are also discussed.

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.

PubChem structure-activity relationship (SAR) clusters

BACKGROUND

Developing structure-activity relationships (SARs) of molecules is an important approach in facilitating hit exploration in the early stage of drug discovery. Although information on millions of compounds and their bioactivities is freely available to the public, it is very challenging to infer a meaningful and novel SAR from that information.

RESULTS

Research discussed in the present paper employed a bioactivity-centered clustering approach to group 843,845 non-inactive compounds stored in PubChem according to both structural similarity and bioactivity similarity, with the aim of mining bioactivity data in PubChem for useful SAR information. The compounds were clustered in three bioactivity similarity contexts: (1) non-inactive in a given bioassay, (2) non-inactive against a given protein, and (3) non-inactive against proteins involved in a given pathway. In each context, these small molecules were clustered according to their two-dimensional (2-D) and three-dimensional (3-D) structural similarities. The resulting 18 million clusters, named "PubChem SAR clusters", were delivered in such a way that each cluster contains a group of small molecules similar to each other in both structure and bioactivity.

CONCLUSIONS

The PubChem SAR clusters, pre-computed using publicly available bioactivity information, make it possible to quickly navigate and narrow down the compounds of interest. Each SAR cluster can be a useful resource in developing a meaningful SAR or enable one to design or expand compound libraries from the cluster. It can also help to predict the potential therapeutic effects and pharmacological actions of less-known compounds from those of well-known compounds (i.e., drugs) in the same cluster.

The receptor concept in 3D: from hypothesis and metaphor to GPCR-ligand structures

The first mentioning of the word "receptor" for the structure with which a bioactive compound should react for obtaining its specific influence on a physiological system goes back to the years around 1900. The receptor concept was adapted from the lock and key theory for the enzyme substrate and blockers interactions. Through the years the concept, in the beginning rather being a metaphor, not a model, was refined and became reality in recent years. Not only the structures of receptors were elucidated, also the receptor machineries were unraveled. Following a brief historical review we will describe how the recent breakthroughs in the experimental determination of G protein-coupled receptor (GPCR) crystal structures can be complemented by computational modeling, medicinal chemistry, biochemical, and molecular pharmacological studies to obtain new insights into the molecular determinants of GPCR-ligand binding and activation. We will furthermore discuss how this information can be used for structure-based discovery of novel GPCR ligands that bind specific (allosteric) binding sites with desired effects on GPCR functional activity.

Inhibition of bacterial virulence: drug-like molecules targeting the Salmonella enterica PhoP response regulator

Two-component signal transduction (TCST) is the predominant signaling scheme used in bacteria to sense and respond to environmental changes in order to survive and thrive. A typical TCST system consists of a sensor histidine kinase to detect external signals and an effector response regulator to respond to external changes. In the signaling scheme, the histidine kinase phosphorylates and activates the response regulator, which functions as a transcription factor to modulate gene expression. One promising strategy toward antibacterial development is to target TCST regulatory systems, specifically the response regulators to disrupt the expression of genes important for virulence. In Salmonella enterica, the PhoQ/PhoP signal transduction system is used to sense and respond to low magnesium levels and regulates the expression for over 40 genes necessary for growth under these conditions, and more interestingly, genes that are important for virulence. In this study, a hybrid approach coupling computational and experimental methods was applied to identify drug-like compounds to target the PhoP response regulator. A computational approach of structure-based virtual screening combined with a series of biochemical and biophysical assays was used to test the predictability of the computational strategy and to characterize the mode of action of the compounds. Eight compounds from virtual screening inhibit the formation of the PhoP-DNA complex necessary for virulence gene regulation. This investigation served as an initial case study for targeting TCST response regulators to modulate the gene expression of a signal transduction pathway important for bacterial virulence. With the increasing resistance of pathogenic bacteria to current antibiotics, targeting TCST response regulators that control virulence is a viable strategy for the development of antimicrobial therapeutics with novel modes of action.

The cytoplasmic rhodopsin-protein interface: potential for drug discovery

The mammalian dim-light photoreceptor rhodopsin is a prototypic G protein coupled receptor (GPCR), interacting with the G protein, transducin, rhodopsin kinase, and arrestin. All of these proteins interact with rhodopsin at its cytoplasmic surface. Structural and modeling studies have provided in-depth descriptions of the respective interfaces. Overlap and thus competition for binding surfaces is a major regulatory mechanism for signal processing. Recently, it was found that the same surface is also targeted by small molecules. These ligands can directly interfere with the binding and activation of the proteins of the signal transduction cascade, but they can also allosterically modulate the retinal ligand binding pocket. Because the pocket that is targeted contains residues that are highly conserved across Class A GPCRs, these findings imply that it may be possible to target multiple GPCRs with the same ligand(s). This is desirable for example in complex diseases such as cancer where multiple GPCRs participate in the disease networks.

Allosteric Modulation of G Protein Coupled Receptors by Cytoplasmic, Transmembrane and Extracellular Ligands

G protein coupled receptors (GPCRs) bind diverse classes of ligands, and depending on the receptor, these may bind in their transmembrane or the extracellular domains, demonstrating the principal ability of GPCRs to bind ligand in either domains. Most recently, it was also observed that small molecule ligands can bind in the cytoplasmic domain, and modulate binding and response to extracellular or transmembrane ligands. Thus, all three domains in GPCRs are potential sites for allosteric ligands, and whether a ligand is allosteric or orthosteric depends on the receptor. Here, we will review the evidence supporting the presence of putative binding pockets in all three domains of GPCRs and discuss possible pathways of communication between these pockets.

Combining 3-D quantitative structure-activity relationship with ligand based and structure based alignment procedures for in silico screening of new hepatitis C virus NS5B polymerase inhibitors

The viral NS5B RNA-dependent RNA-polymerase (RdRp) is one of the best-studied and promising targets for the development of novel therapeutics against hepatitis C virus (HCV). Allosteric inhibition of this enzyme has emerged as a viable strategy toward blocking replication of viral RNA in cell based systems. Herein, we describe how the combination of a complete computational procedure together with biological studies led to the identification of novel molecular scaffolds, hitherto untested toward NS5B polymerase. Structure based 3-D quantitative structure-activity relationship (QSAR) models were generated employing NS5B non-nucleoside inhibitors (NNIs), whose bound conformations were readily available from the protein database (PDB). These were grouped into two training sets of structurally diverse NS5B NNIs, based on their binding to the enzyme thumb (15 NNIs) or palm (10 NNIs) domains. Ligand based (LB) and structure based (SB) alignments were rigorously investigated to assess the reliability on the correct molecular alignment for unknown binding mode modeled compounds. Both Surflex and Autodock programs were able to reproduce with minimal errors the experimental binding conformations of 24 experimental NS5B allosteric inhibitors. Eighty-one (thumb) and 223 (palm) modeled compounds taken from literature were LB and SB aligned and used as external validation sets for the development of 3-D QSAR models. Low error of prediction proved the 3-D QSARs to be useful scoring functions for the in silico screening procedure. Finally, the virtual screening of the NCI Diversity Set led to the selection for enzymatic assays of 20 top-scoring molecules for each final model. Among the 40 selected molecules, preliminary data yielded four derivatives exhibiting IC(50) values ranging between 45 and 75 microM. Binding mode analysis of hit compounds within the NS5B polymerase thumb domain showed that one of them, NSC 123526, exhibited a docked conformation which was in good agreement with the thumb training set most active compound (6).

The membrane complex between transducin and dark-state rhodopsin exhibits large-amplitude interface dynamics on the sub-microsecond timescale: insights from all-atom MD simulations

Rhodopsin, the prototype class A G-protein-coupled receptor, is a very important model system for all seven-transmembrane domain proteins. Characterization of the interactions between rhodopsin and transducin, its intracellular G-protein counterpart, and the fluctuations in these interactions due to thermal motions is required for an understanding of early events in the mechanism of signal transduction. In this study, we used all-atom molecular dynamics simulations of a transmembrane protein complex between rhodopsin and the heterotrimeric transducin (G alpha beta gamma) in an all-atom DOPC (1,2-dioleoylsn-glycero-3-phosphocholine) membrane-water environment. Based on the analysis of a microsecond-timescale simulation trajectory, we characterized the dynamics of the system and its effects in the structural features of the protein subunits. Our simulations describe a highly dynamic interaction interface where the system is alternating between distinct domain orientations at the 10- to 100-ns timescale that can be further classified into interaction modes involving contacts between distinct structural features on the protein subunits. We related our results with experimental measurements from a variety of studies and high-resolution models of activated rhodopsin. Monitoring key structural features that are involved in the activation process along our simulation trajectory indicates the presence of extensive dynamics in the dark-adapted state, including a motion of Y223 from helix 3 toward the "ionic-lock" interactions of the conserved ERY motif. The dynamic picture shown here is consistent with a framework in which the dark-state fluctuations sample conformations consistent with the activated state. These results provide an atomic-level description of the dynamics of the full complex and further suggest novel mutagenesis experiments that can be used to investigate the stability and dynamics of this model membrane protein receptor system.

Using ligand-based virtual screening to allosterically stabilize the activated state of a GPCR

G-protein coupled receptors play an essential role in many biological processes. Despite an increase in the number of solved X-ray crystal structures of G-protein coupled receptors, capturing a G-protein coupled receptor in its activated state for structural analysis has proven to be difficult. An unexplored paradigm is stabilization of one or more conformational states of a G-protein coupled receptor via binding a small molecule to the intracellular loops. A short tetrazole peptidomimetic based on the photoactivated state of rhodopsin-bound structure of Gt(alpha)(340-350) was previously designed and shown to stabilize the photoactivated state of rhodopsin, the G-protein coupled receptor involved in vision. A pharmacophore model derived from the designed tetrazole tetrapeptide was used for ligand-based virtual screening to enhance the possible discovery of novel scaffolds. Maybridge Hitfinder and National Cancer Institute diversity libraries were screened for compounds containing the pharmacophore. Forty-seven compounds resulted from virtually screening the Maybridge library, whereas no hits resulted with the National Cancer Institute library. Three of the 47 Maybridge compounds were found to stabilize the MII state. As these compounds did not inhibit binding of transducin to photoactivated state of rhodopsin, they were assumed to be allosteric ligands. These compounds are potentially useful for crystallographic studies where complexes with these compounds might capture rhodopsin in its activated conformational state.

Recent advances in structure-based virtual screening of G-protein coupled receptors

In addition to the rhodopsin crystal structure, high-resolution crystal structures of ligand-mediated G-protein-coupled receptors (GPCRs) have recently become available, and these have become attractive templates for developing homology models of several GPCRs of therapeutic interest. These crystal structures and the homology models derived from them have provided significant insights into ligand-receptor interactions. Moreover, several studies have demonstrated that the structural models are indeed suitable for virtual screening of compound databases to identify new ligands for various GPCRs. Recent examples of such virtual screening against GPCRs are discussed in this review.