Discovery of agonist-antagonist pairs for the modulation of Ca [2]+ and voltage-gated K+ channels of large conductance that contain beta1 subunits.

Large conductance, calcium/voltage-gated potassium channels (BK) regulate critical body processes, including neuronal, secretory and smooth muscle (SM) function. While BK-forming alpha subunits are ubiquitous, accessory beta1 subunits are highly expressed in SM. This makes beta1 an attractive target for pharmaceutical development to treat SM disorders, such as hypertension or cerebrovascular spasm. Compounds activating BK via beta1 have been identified, yet they exhibit low potency and off-target effects while antagonists that limit agonist activity via beta 1 remain unexplored. Beta1-dependent BK ligand-based pharmacophore modeling and ZINC database searches identified 15 commercially available hits. Concentration-response curves on BK alpha + beta1 subunit-mediated currents were obtained in CHO cells. One potent (EC50 = 20 nM) and highly efficacious activator (maximal activation = ×10.3 of control) was identified along with a potent antagonist (KB = 3.02 nM), both of which were dependent on beta1. Our study provides the first proof-of-principle that an agonist/antagonist pair can be used to control beta1-containing BK activity.

Optical control of sphingosine-1-phosphate formation and function.

Sphingosine-1-phosphate (S1P) plays important roles as a signaling lipid in a variety of physiological and pathophysiological processes. S1P signals via a family of G-protein-coupled receptors (GPCRs) (S1P1-5) and intracellular targets. Here, we report on photoswitchable analogs of S1P and its precursor sphingosine, respectively termed PhotoS1P and PhotoSph. PhotoS1P enables optical control of S1P1-3, shown through electrophysiology and Ca2+ mobilization assays. We evaluated PhotoS1P in vivo, where it reversibly controlled S1P3-dependent pain hypersensitivity in mice. The hypersensitivity induced by PhotoS1P is comparable to that induced by S1P. PhotoS1P is uniquely suited for the study of S1P biology in cultured cells and in vivo because it exhibits prolonged metabolic stability compared to the rapidly metabolized S1P. Using lipid mass spectrometry analysis, we constructed a metabolic map of PhotoS1P and PhotoSph. The formation of these photoswitchable lipids was found to be light dependent, providing a novel approach to optically probe sphingolipid biology.

Optical Control of Lysophosphatidic Acid Signaling.

Lysophosphatidic acid (LPA) is a phospholipid that acts as an extracellular signaling molecule and activates the family of lysophosphatidic acid receptors (LPA1-6). These G protein-coupled receptors (GPCRs) are broadly expressed and are particularly important in development as well as in the nervous, cardiovascular, reproductive, gastrointestinal, and pulmonary systems. Here, we report on a photoswitchable analogue of LPA, termed AzoLPA, which contains an azobenzene photoswitch embedded in the acyl chain. AzoLPA enables optical control of LPA receptor activation, shown through its ability to rapidly control LPA-evoked increases in intracellular Ca2+ levels. AzoLPA shows greater activation of LPA receptors in its light-induced cis-form than its dark-adapted (or 460 nm light-induced) trans-form. AzoLPA enabled the optical control of neurite retraction through its activation of the LPA2 receptor.

Benchmarking GPCR homology model template selection in combination with de novo loop generation.

G protein-coupled receptors (GPCR) comprise the largest family of membrane proteins and are of considerable interest as targets for drug development. However, many GPCR structures remain unsolved. To address the structural ambiguity of these receptors, computational tools such as homology modeling and loop modeling are often employed to generate predictive receptor structures. Here we combined both methods to benchmark a protocol incorporating homology modeling based on a locally selected template and extracellular loop modeling that additionally evaluates the presence of template ligands during these modeling steps. Ligands were also docked using three docking methods and two pose selection methods to elucidate an optimal ligand pose selection method. Results suggest that local template-based homology models followed by loop modeling produce more accurate and predictive receptor models than models produced without loop modeling, with decreases in average receptor and ligand RMSD of 0.54 Å and 2.91 Å, respectively. Ligand docking results showcased the ability of MOE induced fit docking to produce ligand poses with atom root-mean-square deviation (RMSD) values at least 0.20 Å lower (on average) than the other two methods benchmarked in this study. In addition, pose selection methods (software-based scoring, ligand complementation) selected lower RMSD poses with MOE induced fit docking than either of the other methods (averaging at least 1.57 Å lower), indicating that MOE induced fit docking is most suited for docking into GPCR homology models in our hands. In addition, target receptor models produced with a template ligand present throughout the modeling process most often produced target ligand poses with RMSD values ≤ 4.5 Å and Tanimoto coefficients > 0.6 after selection based on ligand complementation than target receptor models produced in the absence of template ligands. Overall, the findings produced by this study support the use of local template homology modeling in combination with de novo ECL2 modeling in the presence of a ligand from the template crystal structure to generate GPCR models intended to study ligand binding interactions.

Molecular modelling guided design, synthesis and QSAR analysis of new small molecule non-lipid autotaxin inhibitors.

The lysophospholipase D autotaxin (ATX) generates lysophosphatidic acid (LPA) that activates six cognate G-protein coupled receptors (GPCR) in cancerous cells, promoting their motility and invasion. Four novel compounds were generated aided by molecular docking guided design and synthesis techniques to obtain new dual inhibitors of ATX and the lysophosphatidic acid receptor subtype 1 (LPAR1). Biological evaluation of these compounds revealed two compounds, 10 and 11, as new ATX enzyme inhibitors with potencies in the range of 218-220 nM and water solubility (>100 µg/mL), but with no LPAR1 inhibitory activity. A QSAR model was generated that included four newly designed compounds and twenty-one additional compounds that we have reported previously. The QSAR model provided excellent predictability of the pharmacological activity and potency among structurally related drug candidates. This model will be highly useful in guiding the synthesis of new ATX inhibitors in the future.

A benchmark study of loop modeling methods applied to G protein-coupled receptors.

G protein-coupled receptors (GPCR) are important drug discovery targets. Despite progress, many GPCR structures have not yet been solved. For these targets, comparative modeling is used in virtual ligand screening to prioritize experimental efforts. However, the structure of extracellular loop 2 (ECL2) is often poorly predicted. This is significant due to involvement of ECL2 in ligand binding for many Class A GPCR. Here we examine the performance of loop modeling protocols available in the Rosetta (cyclic coordinate descent [CCD], KIC with fragments [KICF] and next generation KIC [NGK]) and Molecular Operating Environment (MOE) software suites (de novo search). ECL2 from GPCR crystal structures served as the structure prediction targets and were divided into four sets depending on loop length. Results suggest that KICF and NGK sampled and scored more loop models with sub-angstrom and near-atomic accuracy than CCD or de novo search for loops of 24 or fewer residues. None of the methods were able to sample loop conformations with near-atomic accuracy for the longest targets ranging from 25 to 32 residues based on 1000 models generated. For these long loop targets, increased conformational sampling is necessary. The strongly conserved disulfide bond between Cys3.25 and Cys45.50 in ECL2 proved an effective filter. Setting an upper limit of 5.1 Å on the S-S distance improved the lowest RMSD model included in the top 10 scored structures in Groups 1-4 on average between 0.33 and 1.27 Å. Disulfide bond formation and geometry optimization of ECL2 provided an additional incremental benefit in structure quality.

Phospholipase D2-dependent inhibition of the nuclear hormone receptor PPARgamma by cyclic phosphatidic acid.

Cyclic phosphatidic acid (1-acyl-2,3-cyclic-glycerophosphate, CPA), one of nature's simplest phospholipids, is found in cells from slime mold to humans and has a largely unknown function. We find here that CPA is generated in mammalian cells in a stimulus-coupled manner by phospholipase D2 (PLD2) and binds to and inhibits the nuclear hormone receptor PPARgamma with nanomolar affinity and high specificity through stabilizing its interaction with the corepressor SMRT. CPA production inhibits the PPARgamma target-gene transcription that normally drives adipocytic differentiation of 3T3-L1 cells, lipid accumulation in RAW264.7 cells and primary mouse macrophages, and arterial wall remodeling in a rat model in vivo. Inhibition of PLD2 by shRNA, a dominant-negative mutant, or a small molecule inhibitor blocks CPA production and relieves PPARgamma inhibition. We conclude that CPA is a second messenger and a physiological inhibitor of PPARgamma, revealing that PPARgamma is regulated by endogenous agonists as well as by antagonists.

Discovery and synthetic optimization of a novel scaffold for hydrophobic tunnel-targeted autotaxin inhibition.

Autotaxin (ATX) is a ubiquitous ectoenzyme that hydrolyzes lysophosphatidylcholine (LPC) to form the bioactive lipid mediator lysophosphatidic acid (LPA). LPA activates specific G-protein coupled receptors to elicit downstream effects leading to cellular motility, survival, and invasion. Through these pathways, upregulation of ATX is linked to diseases such as cancer and cardiovascular disease. Recent crystal structures confirm that the catalytic domain of ATX contains multiple binding regions including a polar active site, hydrophobic tunnel, and a hydrophobic pocket. This finding is consistent with the promiscuous nature of ATX hydrolysis of multiple and diverse substrates and prior investigations of inhibitor impacts on ATX enzyme kinetics. The current study used virtual screening methods to guide experimental identification and characterization of inhibitors targeting the hydrophobic region of ATX. An initially discovered inhibitor, GRI392104 (IC50 4µM) was used as a lead for synthetic optimization. In total twelve newly synthesized inhibitors of ATX were more potent than GRI392104 and were selective for ATX as they had no effect on other LPC-specific NPP family members or on LPA1-5 GPCR.

Activation of calcium- and voltage-gated potassium channels of large conductance by leukotriene B4.

Calcium/voltage-gated, large conductance potassium (BK) channels control numerous physiological processes, including myogenic tone. BK channel regulation by direct interaction between lipid and channel protein sites has received increasing attention. Leukotrienes (LTA4, LTB4, LTC4, LTD4, and LTE4) are inflammatory lipid mediators. We performed patch clamp studies in Xenopus oocytes that co-expressed BK channel-forming (cbv1) and accessory ß1 subunits cloned from rat cerebral artery myocytes. Leukotrienes were applied at 0.1 nm-10 µm to either leaflet of cell-free membranes at a wide range of [Ca(2+)]i and voltages. Only LTB4 reversibly increased BK steady-state activity (EC50 = 1 nm; Emax reached at 10 nm), with physiological [Ca(2+)]i and voltages favoring this activation. Homomeric cbv1 or cbv1-ß2 channels were LTB4-resistant. Computational modeling predicted that LTB4 docked onto the cholane steroid-sensing site in the BK ß1 transmembrane domain 2 (TM2). Co-application of LTB4 and cholane steroid did not further increase LTB4-induced activation. LTB4 failed to activate ß1 subunit-containing channels when ß1 carried T169A, A176S, or K179I within the docking site. Co-application of LTB4 with LTA4, LTC4, LTD4, or LTE4 suppressed LTB4-induced activation. Inactive leukotrienes docked onto a portion of the site, probably preventing tight docking of LTB4. In summary, we document the ability of two endogenous lipids from different chemical families to share their site of action on a channel accessory subunit. Thus, cross-talk between leukotrienes and cholane steroids might converge on regulation of smooth muscle contractility via BK ß1. Moreover, the identification of LTB4 as a highly potent ligand for BK channels is critical for the future development of ß1-specific BK channel activators.

Integrating the puzzle pieces: the current atomistic picture of phospholipid-G protein coupled receptor interactions.

A compelling question of how phospholipids interact with their target receptors has been of interest since the first receptor-mediated effects were reported. The recent report of a crystal structure for the S1P(1) receptor in complex with an antagonist phospholipid provides interesting perspective on the insights that had previously been gained through structure-activity studies of the phospholipids, as well as modeling and mutagenesis studies of the receptors. This review integrates these varied lines of investigation in the context of their various contributions to our current understanding of phospholipid-receptor interactions. This article is part of a Special Issue entitled Advances in Lysophospholipid Research.

The steroid interaction site in transmembrane domain 2 of the large conductance, voltage- and calcium-gated potassium (BK) channel accessory ß1 subunit.

Large conductance, voltage- and calcium-gated potassium (BK) channels regulate several physiological processes, including myogenic tone and thus, artery diameter. Nongenomic modulation of BK activity by steroids is increasingly recognized, but the precise location of steroid action remains unknown. We have shown that artery dilation by lithocholate (LC) and related cholane steroids is caused by a 2× increase in vascular myocyte BK activity (EC(50) = 45 µM), an action that requires ß1 but not other (ß2-ß4) BK accessory subunits. Combining mutagenesis and patch-clamping under physiological conditions of calcium and voltage on BK α- (cbv1) and ß1 subunits from rat cerebral artery myocytes, we identify the steroid interaction site from two regions in BK ß1 transmembrane domain 2 proposed by computational dynamics: the outer site includes L157, L158, and T165, whereas the inner site includes T169, L172, and L173. As expected from computational modeling, cbv1+rß1T165A,T169A channels were LC-unresponsive. However, cbv1 + rß1T165A and cbv1 + rß1T165A,L157A,L158A were fully sensitive to LC. Data indicate that the transmembrane domain 2 outer site does not contribute to steroid action. Cbv1 + rß1T169A was LC-insensitive, with rß1T169S being unable to rescue responsiveness to LC. Moreover, cbv1 + rß1L172A, and cbv1 + rß1L173A channels were LC-insensitive. These data and computational modeling indicate that tight hydrogen bonding between T169 and the steroid α-hydroxyl, and hydrophobic interactions between L172,L173 and the steroid rings are both necessary for LC action. Therefore, ß1 TM2 T169,L172,L173 provides the interaction area for cholane steroid activation of BK channels. Because this amino acid triplet is unique to BK ß1, our study provides a structural basis for advancing ß1 subunit-specific pharmacology of BK channels.

Cerebrovascular dilation via selective targeting of the cholane steroid-recognition site in the BK channel ß1-subunit by a novel nonsteroidal agent.

The Ca(2+)/voltage-gated K(+) large conductance (BK) channel ß1 subunit is particularly abundant in vascular smooth muscle. By determining their phenotype, BK ß1 allows the BK channels to reduce myogenic tone, facilitating vasodilation. The endogenous steroid lithocholic acid (LCA) dilates cerebral arteries via BK channel activation, which requires recognition by a BK ß1 site that includes Thr169. Whether exogenous nonsteroidal agents can access this site to selectively activate ß1-containing BK channels and evoke vasodilation remain unknown. We performed a chemical structure database similarity search using LCA as a template, along with a two-step reaction to generate sodium 3-hydroxyolean-12-en-30-oate (HENA). HENA activated the BK (cbv1 + ß1) channels cloned from rat cerebral artery myocytes with a potency (EC50 = 53 µM) similar to and an efficacy (×2.5 potentiation) significantly greater than that of LCA. This HENA action was replicated on native channels in rat cerebral artery myocytes. HENA failed to activate the channels made of cbv1 + ß2, ß3, ß4, or ß1T169A, indicating that this drug selectively targets ß1-containing BK channels via the BK ß1 steroid-sensing site. HENA (3-45 µM) dilated the rat and C57BL/6 mouse pressurized cerebral arteries. Consistent with the electrophysiologic results, this effect was larger than that of LCA. HENA failed to dilate the arteries from the KCNMB1 knockout mouse, underscoring BK ß1's role in HENA action. Finally, carotid artery-infusion of HENA (45 µM) dilated the pial cerebral arterioles via selective BK-channel targeting. In conclusion, we have identified for the first time a nonsteroidal agent that selectively activates ß1-containing BK channels by targeting the steroid-sensing site in BK ß1, rendering vasodilation.

Hits of a high-throughput screen identify the hydrophobic pocket of autotaxin/lysophospholipase D as an inhibitory surface.

Autotaxin (ATX), a lysophospholipase D, plays an important role in cancer invasion, metastasis, tumor progression, tumorigenesis, neuropathic pain, fibrotic diseases, cholestatic pruritus, lymphocyte homing, and thrombotic diseases by producing the lipid mediator lysophosphatidic acid (LPA). A high-throughput screen of ATX inhibition using the lysophosphatidylcholine-like substrate fluorogenic substrate 3 (FS-3) and ∼10,000 compounds from the University of Cincinnati Drug Discovery Center identified several small-molecule inhibitors with IC50 vales ranging from nanomolar to low micromolar. The pharmacology of the three most potent compounds: 918013 (1; 2,4-dichloro-N-(3-fluorophenyl)-5-(4-morpholinylsulfonyl) benzamide), 931126 (2; 4-oxo-4-{2-[(5-phenoxy-1H-indol-2-yl)carbonyl]hydrazino}-N-(4-phenylbutan-2-yl)butanamide), and 966791 (3; N-(2,6-dimethylphenyl)-2-[N-(2-furylmethyl)(4-(1,2,3,4-tetraazolyl)phenyl)carbonylamino]-2-(4-hydroxy-3-methoxyphenyl) acetamide), were further characterized in enzyme, cellular, and whole animal models. Compounds 1 and 2 were competitive inhibitors of ATX-mediated hydrolysis of the lysophospholipase substrate FS-3. In contrast, compound 3 was a competitive inhibitor of both FS-3 and the phosphodiesterase substrate p-nitrophenyl thymidine 5'-monophosphate. Computational docking and mutagenesis suggested that compounds 1 and 2 target the hydrophobic pocket, thereby blocking access to the active site of ATX. The potencies of compounds 1-3 were comparable to each other in each of the assays. All of these compounds significantly reduced invasion of A2058 human melanoma cells in vitro and the colonization of lung metastases by B16-F10 murine melanoma cells in C57BL/6 mice. The compounds had no agonist or antagonist effects on select LPA or sphingosine 1-phosphate receptors, nor did they inhibit nucleotide pyrophosphatase/phosphodiesterase (NPP) enzymes NPP6 and NPP7. These results identify the molecular surface of the hydrophobic pocket of ATX as a target-binding site for inhibitors of enzymatic activity.

Multiple cholesterol recognition/interaction amino acid consensus (CRAC) motifs in cytosolic C tail of Slo1 subunit determine cholesterol sensitivity of Ca2+- and voltage-gated K+ (BK) channels.

Large conductance, Ca(2+)- and voltage-gated K(+) (BK) channel proteins are ubiquitously expressed in cell membranes and control a wide variety of biological processes. Membrane cholesterol regulates the activity of membrane-associated proteins, including BK channels. Cholesterol modulation of BK channels alters action potential firing, colonic ion transport, smooth muscle contractility, endothelial function, and the channel alcohol response. The structural bases underlying cholesterol-BK channel interaction are unknown. Such interaction is determined by strict chemical requirements for the sterol molecule, suggesting cholesterol recognition by a protein surface. Here, we demonstrate that cholesterol action on BK channel-forming Cbv1 proteins is mediated by their cytosolic C tail domain, where we identified seven cholesterol recognition/interaction amino acid consensus motifs (CRAC4 to 10), a distinct feature of BK proteins. Cholesterol sensitivity is provided by the membrane-adjacent CRAC4, where Val-444, Tyr-450, and Lys-453 are required for cholesterol sensing, with hydrogen bonding and hydrophobic interactions participating in cholesterol location and recognition. However, cumulative truncations or Tyr-to-Phe substitutions in CRAC5 to 10 progressively blunt cholesterol sensitivity, documenting involvement of multiple CRACs in cholesterol-BK channel interaction. In conclusion, our study provides for the first time the structural bases of BK channel cholesterol sensitivity; the presence of membrane-adjacent CRAC4 and the long cytosolic C tail domain with several other CRAC motifs, which are not found in other members of the TM6 superfamily of ion channels, very likely explains the unique cholesterol sensitivity of BK channels.

Gold nanorods carrying paclitaxel for photothermal-chemotherapy of cancer.

Nanotechnology-based photothermal therapy has emerged as a promising treatment for cancer during the past decade. However, heterogeneous laser heating and limited light penetration can lead to incomplete tumor cell eradication. Here, we developed a method to overcome these limitations by combining chemotherapy with photothermal therapy using paclitaxel-loaded gold nanorods. Paclitaxel was loaded to gold nanorods with high density (2.0 × 10(4) paclitaxel per gold nanorod) via nonspecific adsorption, followed by stabilization with poly(ethylene glycol) linked with 11-mercaptoundecanoic acid. Paclitaxel was entrapped in the hydrophobic pocket of the polymeric monolayer on the surface of gold nanorods, which allows direct cellular delivery of the hydrophobic drugs via the lipophilic plasma membrane. Highly efficient drug release was demonstrated in a cell membrane mimicking two-phase solution. Combined photothermal therapy and chemotherapy with the paclitaxel-loaded gold nanorods was shown to be highly effective in killing head and neck cancer cells and lung cancer cells, superior to photothermal therapy or chemotherapy alone due to a synergistic effect. The paclitaxel-gold nanorod enabled photothermal chemotherapy has the potential of preventing tumor reoccurrence and metastasis and may have an important impact on the treatment of head and neck cancer and other malignancies in the clinic.

Autotaxin inhibition: development and application of computational tools to identify site-selective lead compounds.

Autotaxin (ATX) catalyzes the conversion of lysophosphatidyl choline (LPC) to lysophosphatidic acid (LPA). Both ATX and LPA have been linked to pathophysiologies ranging from cancer to neuropathic pain. Inhibition of LPA production by ATX is therefore of therapeutic interest. Here we report the application of previously-developed, subsite-targeted pharmacophore models in a screening workflow that involves either docking or binary QSAR as secondary filters to identify ATX inhibitors from previously unreported structural types, four of which have sub-micromolar inhibition constants. Cell-based assays demonstrate that ATX inhibition and cytotoxicity structure-activity-relationships (SAR) exhibit selectivity cliffs, characterized by structurally similar compounds exhibiting similar biological activities with respect to ATX inhibition but very different biological activities with respect to cytotoxicity. Thus, general cytotoxicity should not be used as an early filter to eliminate candidate ATX inhibitor scaffolds from further SAR studies. Assays using two substrates of vastly different sizes demonstrate that the tools developed to identify compounds binding outside the central core of the active site did identify compounds acting at an allosteric site. In contrast, tools developed to identify active-site directed compounds did not identify active-site directed compounds. The stronger volume overlap imposed when selecting screening candidates expected to bind outside the active site is likely responsible for the stronger match between intended and actual target site.

Structural Characterization of an LPA1 Second Extracellular Loop Mimetic with a Self-Assembling Coiled-Coil Folding Constraint.

G protein-coupled receptor (GPCR) structures are of interest as a means to understand biological signal transduction and as tools for therapeutic discovery. The growing number of GPCR crystal structures demonstrates that the extracellular loops (EL) connecting the membrane-spanning helices show tremendous structural variability relative to the more structurally-conserved seven transmembrane α-helical domains. The EL of the LPA(1) receptor have not yet been conclusively resolved, and bear limited sequence identity to known structures. This study involved development of a peptide to characterize the intrinsic structure of the LPA(1) GPCR second EL. The loop was embedded between two helices that assemble into a coiled-coil, which served as a receptor-mimetic folding constraint (LPA(1)-CC-EL2 peptide). The ensemble of structures from multi-dimensional NMR experiments demonstrated that a robust coiled-coil formed without noticeable deformation due to the EL2 sequence. In contrast, the EL2 sequence showed well-defined structure only near its C-terminal residues. The NMR ensemble was combined with a computational model of the LPA(1) receptor that had previously been validated. The resulting hybrid models were evaluated using docking. Nine different hybrid models interacted with LPA 18:1 as expected, based on prior mutagenesis studies, and one was additionally consistent with antagonist affinity trends.

Application of computational methods for class A GPCR Ligand discovery.

G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for drug development due to their role in transmitting cellular signals in a multitude of biological processes. Of the six classes categorizing GPCR (A, B, C, D, E, and F), class A contains the largest number of therapeutically relevant GPCR. Despite their importance as drug targets, many challenges exist for the discovery of novel class A GPCR ligands serving as drug precursors. Though knowledge of the structural and functional characteristics of GPCR has grown significantly over the past 20 years, a large portion of GPCR lack reported, experimentally determined structures. Furthermore, many GPCR have no known endogenous and/or synthetic ligands, limiting further exploration of their biochemical, cellular, and physiological roles. While many successes in GPCR ligand discovery have resulted from experimental high-throughput screening, computational methods have played an increasingly important role in GPCR ligand identification in the past decade. Here we discuss computational techniques applied to GPCR ligand discovery. This review summarizes class A GPCR structure/function and provides an overview of many obstacles currently faced in GPCR ligand discovery. Furthermore, we discuss applications and recent successes of computational techniques used to predict GPCR structure as well as present a summary of ligand- and structure-based methods used to identify potential GPCR ligands. Finally, we discuss computational hit list generation and refinement and provide comprehensive workflows for GPCR ligand identification.

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.

Structure-based pharmacophore modeling 2. Developing a novel framework for structure-based pharmacophore model generation and selection.

Pharmacophore models are three-dimensional arrangements of molecular features required for biological activity that are used in ligand identification efforts for many biological targets, including G protein-coupled receptors (GPCR). Though GPCR are integral membrane proteins of considerable interest as targets for drug development, many of these receptors lack known ligands or experimentally determined structures necessary for ligand- or structure-based pharmacophore model generation, respectively. Thus, we here present a structure-based pharmacophore modeling approach that uses fragments placed with Multiple Copy Simultaneous Search (MCSS) to generate high-performing pharmacophore models in the context of experimentally determined, as well as modeled GPCR structures. Moreover, we have addressed the oft-neglected topic of pharmacophore model selection via development of a cluster-then-predict machine learning workflow. Herein score-based pharmacophore models were generated in experimentally determined and modeled structures of 13 class A GPCR and resulted in pharmacophore models exhibiting high enrichment factors when used to search a database containing 569 class A GPCR ligands. In addition, classification of pharmacophore models with the best performing cluster-then-predict logistic regression classifier resulted in positive predictive values (PPV) of 0.88 and 0.76 for selecting high enrichment pharmacophore models from among those generated in experimentally determined and modeled structures, respectively.