In this work, a high-throughput, intracellular calcium fluorescence assay for 384-well plates to screen small molecule libraries on recombinant G protein-coupled receptors (GPCRs) is described. The target, the kinin receptor from the cattle fever tick, Rhipicephalus microplus, is expressed in CHO-K1 cells. This assay identifies agonists and antagonists using the same cells in one “dual-addition” assay.
G protein-coupled receptors (GPCRs) represent the largest superfamily of receptors and are the targets of numerous human drugs. High-throughput screening (HTS) of random small molecule libraries against GPCRs is used by the pharmaceutical industry for target-specific drug discovery. In this study, an HTS was employed to identify novel small-molecule ligands of invertebrate-specific neuropeptide GPCRs as probes for physiological studies of vectors of deadly human and veterinary pathogens.
The invertebrate-specific kinin receptor was chosen as a target because it regulates many important physiological processes in invertebrates, including diuresis, feeding, and digestion. Furthermore, the pharmacology of many invertebrate GPCRs is poorly characterized or not characterized at all; therefore, the differential pharmacology of these groups of receptors with respect to the related GPCRs in other metazoans, especially humans, adds knowledge to the structure-activity relationships of GPCRs as a superfamily. An HTS assay was developed for cells in 384-well plates for the discovery of ligands of the kinin receptor from the cattle fever tick, or southern cattle tick, Rhipicephalus microplus. The tick kinin receptor was stably expressed in CHO-K1 cells.
The kinin receptor, when activated by endogenous kinin neuropeptides or other small molecule agonists, triggers Ca2+ release from calcium stores into the cytoplasm. This calcium fluorescence assay combined with a "dual-addition" approach can detect functional agonist and antagonist "hit" molecules in the same assay plate. Each assay was conducted using drug plates carrying an array of 320 random small molecules. A reliable Z' factor of 0.7 was obtained, and three agonist and two antagonist hit molecules were identified when the HTS was at a 2 µM final concentration. The calcium fluorescence assay reported here can be adapted to screen other GPCRs that activate the Ca2+ signaling cascade.
G protein-coupled receptors (GPCRs), which are present from yeast to humans, represent the largest superfamily of receptors in many organisms1. They play critical roles in regulating nearly all biological processes in animals. There are 50-200 GPCRs in the genome of arthropods, meaning they represent the largest membrane receptor superfamily2. They are classified into six major classes, A-F, based on their sequence similarity and functions3. GPCRs transduce various extracellular signals, such as those of hormones, neuropeptides, biogenic amines, glutamate, proton, lipoglycoproteins, and photons4. GPCRs couple to heterotrimer G proteins (Gα, Gβ, and Gγ) to transmit downstream signals. GPCRs coupled to Gαs or Gαi/o proteins increase or decrease, respectively, the intracellular 3', 5'-cyclic adenosine monophosphate (cAMP) levels by activating or inhibiting adenylyl cyclase. GPCRs coupled to Gαq/11 induce calcium release from the endoplasmic reticulum calcium stores by activating the phospholipase C (PLC)-inositol-1,4,5-triphosphate (IP3) pathway. GPCRs coupled to Gα12/13 activate RhoGTPase nucleotide exchange factors5,6. GPCRs are the target of more than 50% of human drugs and an acaricide, amitraz4. As GPCRs transduce such diverse signals, they are promising targets for developing novel pesticides that disrupt invertebrate-specific physiological functions.
The goal of HTS is to identify hit molecules that can modulate receptor functions. HTS involves assay development, miniaturization, and automation7. Arthropod neuropeptide GPCRs are involved in most physiological functions, such as development, molting and ecdysis, excretion, energy mobilization, and reproduction4. Most of the neuropeptide GPCRs of arthropods and metazoans signal through the calcium signaling cascade2,6,8,9,10, such as in the myoinhibitory peptide and SIFamide receptors of the blacklegged tick Ixodes scapularis; their ligands are antagonistic in hindgut motility assays, with SIF eliciting contraction and MIP inhibiting it11,12. An NPY-like receptor of the yellow fever mosquito, Aedes aegypti, regulates female host seeking13. Compared to other alternative calcium mobilization assays such as the aequorin calcium bioluminescence assay14, the calcium fluorescence assay is easy to perform, does not require the transfection of other recombinant calcium detecting proteins, and is cost-effective. The calcium fluorescence assay produces a prolonged signal compared to the fast kinetic signal obtained in the aequorin calcium bioluminescence assay14,15.
In the example here, the kinin receptor from the cattle fever tick, Rhipicephalus microplus, was recombinantly expressed in the CHO-K1 cell line and used for the calcium fluorescence assay. There is only one kinin receptor gene found in R. microplus; the receptor signals through a Gq protein-dependent signaling pathway and triggers the efflux of Ca2+ from calcium stores into the intracellular space16. This process can be detected and quantified by a fluorophore, which elicits a fluorescence signal when binding calcium ions (Figure 1).
The kinin receptor is an invertebrate-specific GPCR, which belongs to the Class A Rhodopsin-like receptors. Kinin is an ancient signaling neuropeptide that is present in Mollusca, Crustacea, Insecta, and Acari4,17,18. Coleopterans (beetles) lack the kinin signaling system; in the mosquito Aedes aegypti, there is only one kinin receptor that binds three aedeskinins, while Drosophila melanogaster has one kinin receptor with drosokinin as a unique ligand19,20,21. There are no homologous kinins or kinin receptors in vertebrates. Although the exact function of kinin is unknown in ticks, the kinin receptor RNAi-silenced females of R. microplus show significantly reduced reproductive fitness22. Kinins are pleotropic peptides in insects. In Drosophila melanogaster, they are involved in both the central and peripheral nervous regulatory systems23, pre-ecdysis24, feeding25, metabolism26, and sleep activity patterns26,27, as well as larval locomotion28. Kinins regulate hindgut contraction, diuresis, and feeding in the mosquito A. aegypti29,30,31. The kinin peptides have a conserved C-terminal pentapeptide Phe-X1-X2-Trp-Gly-NH2, which is the minimum required sequence for biological activity32. The arthropod specificity, the small size of the endogenous ligand, which makes them amenable to small-molecule interference, and the pleiotropic functions in insects make the kinin receptor a promising target for pest control4.
The "dual-addition" assay (Figure 2) allows the identification of agonists or antagonists in the same HTS assay15. It is adapted from a "dual-addition" assay that is commonly used in the pharmaceutical industry for drug discovery33. In brief, the first addition of drugs into the cell plate allows the identification of potential agonists in the chemical library when a higher fluorescence signal is detected compared to the application of the solvent control. After 5 min of incubation with these small molecules, a known agonist (kinin peptide) is applied to all the wells. Those wells that randomly received an antagonist from the drug plate display a lower fluorescence signal upon agonist addition compared to the control wells that received the solvent in the first addition. This assay then allows the identification of potential agonists and antagonists with the same cells. In a standard HTS project, these hit molecules would be further validated through dose-response assays and by additional biological activity assays, which are not shown here.
Figure 1: Illustration of the calcium fluorescence assay mechanism. The Gq protein triggers the intracellular calcium signaling pathway. The kinin receptor (G protein-coupled receptor) was recombinantly expressed in CHO-K1 cells. When the agonist ligand binds to the receptor, the Gq protein associated with the kinin receptor activates PLC, which catalyzes the conversion of a PIP2 molecule into IP3 and DAG. IP3 then binds to the IP3R on the surface of the endoplasmic reticulum, leading to the release of Ca2+ into the cytoplasm, where Ca2+ ions bind to the fluorophores and elicit a fluorescence signal. The fluorescence signal can be obtained by excitation at 490 nm and detected at 514 nm. Abbreviations: GPCR = G protein-coupled receptor; PLC = phospholipase C; PIP2 = phosphatidylinositol 4,5-bisphosphate; IP3 = inositol trisphosphate; DAG = diacylglycerol; IP3R = IP3 receptor. Created with BioRender.com. Please click here to view a larger version of this figure.
Figure 2: The workflow for the high-throughput screening of small molecules on a G protein-coupled receptor expressed in CHO-K1 cells. (A) Recombinant CHO-K1 cells stably expressing the kinin receptor were added to the 384-well plate (10,000 cells/well) using a liquid handling system (25 µL/well) and incubated in a humidified CO2 incubator for 12-16 h. (B) The assay buffer containing the fluorescent dye (25 µL/well) was added into the cell plate using a liquid handling system. The plate was incubated for 30 min at 37 °C for 30 min and equilibrated at RT for another 30 min. (C) The background fluorescence signal of the cells in each well was measured with a plate reader. (D) Drug solutions from a 384-well library plate and blank solvent (all at 0.5 µL/well) were added into the cellular assay plate using a liquid handling system. (E) Cellular calcium fluorescence responses were measured with the plate reader immediately after the addition of the drug solutions; compound(s) eliciting higher than average fluorescence signals were picked out as agonist hit(s). Antagonist hits that block the GPCR (icon below) were revealed after the addition of the peptide agonist during step G. (F) In the same assay plate, after 5 min of incubation of the cells with screening compounds, an endogenous agonist peptide Rhimi-K-1 (QFSPWGamide) of the tick kinin receptor was added to each well (1 µM). (G) Cellular fluorescence responses after the agonist peptide addition were measured by the plate reader immediately. Compound(s) inhibiting the fluorescence signal were selected as antagonist hit(s). Abbreviations: GPCR = G protein-coupled receptor; RT = room temperature; RFU = relative fluorescence units. Created with BioRender.com. Please click here to view a larger version of this figure.
1. Cell maintenance
NOTE: A CHO-K1 cell line that stably expresses the kinin receptor from R. microplus, named BMLK3, was developed by Holmes et al.16. The details of the development of the cell line are presented elsewhere14. All the following steps are performed under sterile conditions in a class II biosafety cabinet.
2. Calcium fluorescence assay
3. Data analysis
An in-house drug plate (SAC2-34-6170) composed of 320 random small molecules was used for demonstrating this HTS assay as an example. The HTS had excellent assay quality with a Z' factor of 0.7 (Table 1). This Z' factor reflects the assay quality independent of the tested compounds34. A Z′-factor of 0.5 or greater indicates a good assay signal dynamic range between the RFUs of the positive controls and the negative controls. Assays with a Z' factor less than 0.5 have a suboptimal signal dynamic range, are usually not informative for hit selection, and are, therefore, discarded. The response (measured in relative fluorescence units, RFU) from positive controls (average RFU of 64 positive controls responses, µc+ = 130,139) was on average 19 times higher than that of the negative controls (average RFU of 64 negative controls responses, µc- = 6,675). The cutoff fluorescence for selecting "agonist hits" here was set to a normalized percent activation (NPA) max ratio of >50%. The cutoff fluorescence for selecting "antagonist hits" was set to an observed inhibitory activity (Io) of >50%. Three agonist hit molecules and two antagonist hits were selected: SACC-0090237 (NPA = 153%), SACC-0036982 (NPA = 118%), and SACC-0006532 (NPA = 152%) as agonists, and SACC-0103392 (Io = 53%) and SACC-0041280 (Io = 56%) as antagonists (the raw data are shown in Supplemental Table S2). Examples of calculations of NPA and Io for agonists (green rows) and antagonists (orange rows), respectively, are shown in Supplemental Table S2.
Compound addition | µs1 | σs | µc- | σc- | Z’ value |
(n = 320) | (n = 320) | (n = 64) | (n = 64) | ||
9,212 | 18,130 | 6,675 | 1,895 | 0.7 | |
Agonist peptide addition | µs | σs | µc+ | σc+ | |
1,22,322 | 15,987 | 1,30,139 | 11,062 |
Table 1: Summary of the normalized cell responses from the whole plate in the HTS assay. The mean (µ) and standard deviation (σ) are obtained from the CDD vault analysis. Abbreviations: HTS = high-throughput screening; RFUs = relative fluorescence units; 1s (samples) = 320 tested compounds; µs = average value of cellular responses (RFUs) from 320 wells; µc- = average value of RFUs from negative controls; µc+ = average value of RFUs from positive controls; σs = standard deviation of cellular responses (RFUs) from 320 wells; σc- = standard deviation of RFUs from negative controls; σc+ = standard deviation of RFUs from positive controls.
Supplemental Table S1: Custom programs of the liquid handling system for the automatic pipetting steps. This table was modified from Xiong et al.31. Please click here to download this File.
Supplemental Table S2: The raw data of the representative results. The rows highlighted in green represent the agonist hits, and the rows highlighted in orange represent the antagonist hits. Please click here to download this File.
The goal of HTS is to identify hit molecules through screening massive numbers of small molecules. Therefore, the results from this example only represent a small part of a conventional HTS experiment. Furthermore, the hit molecules identified need to be validated in downstream assays such as a dose-dependent assay on the same recombinant cell line and on a CHO-K1 cell line carrying only the empty vector, which can be performed simultaneously to save small molecules. Cytotoxicity assays will help demonstrate that a lack of response in the second addition is not due to cell mortality. Other bioassays should be performed to validate activity in the isolated tissues in vitro, or by applying or feeding the small molecules to live arthropods. The complete pipeline of small-molecule ligand identification of the tick kinin receptor can be found in a previous publication31. Through this pipeline, a total of 36 antagonist hits of the tick kinin receptor were identified. Three of these antagonists were further tested in a tissue contraction-inhibition assay and showed anti-myotropic activity on the mosquito hindgut, where orthologous mosquito kinin receptors are expressed31.
Each GPCR-expressing cell line should be validated for suitable experimental conditions before HTS is performed. These conditions include the cell density (10,000 cells per well), the concentration of the screened compounds (2 µM), the final DMSO concentration (0.2%) in the assay, the concentration of the agonist ligand used for the second addition (1 µM), the pipetting speed and depth, and the reading time post the addition of the agonist peptide (5 min). It is critical to adjust the depth and speed of the liquid handling system to ensure that the dispensing steps do not disturb the cells if the cells are adherent. The validation of these conditions before embarking on HTS will aid in troubleshooting if the assay quality is poor. It is important to have the cells in optimal conditions: the cells should reach 90%-100% confluency but not be overconfluent in the cell plate when used for HTS. If the background signals in certain wells are abnormally low or high, after confirming under the microscope for obviously low or high cell density, the outliers should be removed manually when analyzing the data. Before adapting the calcium fluorescence assay to HTS, the kinetic curve of the positive peptide must be known, as the kinetic response of the cells to the agonist peptide should last for 1-2 min for it to be adapted into an HTS assay and last long enough to read the whole plate. To achieve the prolonged signal, a relatively high concentration of the agonist peptide should be applied, or a higher cell density should be tested. If the kinetic response of the certain receptor cells is too short (lasts less than 1 min), it is recommended to read at a faster speed in the end-point readout or use only a portion of the plate for HTS. Certainly, the quality of the recombinant cell line with respect to the level of receptor expression is crucial for success. Different clonal cell lines will likely vary, as the insertion of the plasmid is random in the genome, and this could certainly influence receptor expression. It would be ideal to incorporate an early calcium fluorescence test of the transfected pool and during the process of selecting clonal cell lines.
The limitation of the calcium fluorescence assay is that it measures the activity of the secondary messengers instead of the direct binding of the ligand to the receptor. Although it is possible to exclude the off-target agonist hits by screening the hit molecules on the cells that do not express the target receptor, it is impossible to completely exclude the off-target antagonist hit molecules using this assay because it is possible that some molecules may inhibit the cellular calcium release pathways downstream of the GPCR signaling. Therefore, further validation is needed using either a ligand binding assay or a bioassay for validation at the organism level, such as in a tick tissue assay with a known agonist, as done previously with mosquitos31.
Most functional assays that detect GPCR activation primarily measure intracellular events in two major signaling pathways: G protein-dependent and G protein-independent pathways. G protein-dependent assays measure Ca2+, inositol triphosphate (IP3), or cAMP concentrations, while G protein-independent assays detect primarily receptor trafficking or β-arrestin recruitment7. Different assays utilize different technologies for readouts, such as Ca2+-sensitive dyes (fluorescence assays), photoproteins such as aequorins (Ca2+ bioluminescence assays)14, tagged fluorescence proteins, fluorescence resonance energy transfer (FRET)-based assays, bioluminescence resonance energy transfer (BRET)-based assays, or label-free optical biosensor-based cellular assays. These assays and the corresponding available readout technologies were reviewed in detail by Zhang and Xie35. Some assays are applicable to any GPCRs, and some are specific to certain GPCRs (e.g., the Ca2+ assay for Gq-coupled receptors, cAMP for Gi/o– or Gs– coupled receptors). However, by promiscuously co-expressing certain G proteins with the receptor in the recombinant cells, the signaling pathway can be altered. For example, Gq15/16 proteins can be co-expressed with non-Gq-coupled receptors for PLC-IP3-Ca2+ signaling detection7,36. In contrast to the assays that measure a single molecular event, there are high-content analysis assays, which collect a massive amount of biological data to improve the ligand specificity37. This type of assay is also more complex and expensive to perform7.
In invertebrate research, the most widely used functional assays detect GPCR activation by analyzing the fluctuations in the concentrations of second messengers, such as Ca2+ or cAMP. Due to the differences in the signaling pathways between mammals and invertebrates, most of the other abovementioned readout technologies have not been yet broadly adopted for invertebrate GPCR studies. However, G-proteins are relatively conserved across the animal kingdom. Recently, Lismont et al.38 demonstrated that some human BRET-based G-protein biosensors can also be adapted to detect the activation of an insect GPCR when co-expressed in the receptor recombinant cells. By co-expressing each of the different G protein-based biosensors (Gs, Gi/o, Gq/11, or G12/13) with specific GPCRs, the G protein-coupling mechanism can be elucidated. In the study by Lismont et al.38, the Schgr-CRF-DH could dose-dependently activate the Gi/o and Gq/11 biosensors and did not activate the Gs and G12/13 biosensors.
Overall, the calcium fluorescence assay performed in a 384-well plate in combination with the "dual-addition" approach described here allows the screening of tens of thousands of molecules in a relatively short time. This is a relatively time- and cost-efficient approach to discover hit molecules for downstream validations and applications.
The authors have nothing to disclose.
This work was supported by the USDA-NIFA-AFRI Animal Health and Well-Being Award (Award number 2022-67015-36336, PVP [Project Director]) and from competitive funds from the Texas A&M AgriLife Research Insect Vector Diseases Grant Program (FY'22-23) to P.V.P. The A.W.E.S.O.M.E. faculty group of the College of Agriculture and Life Sciences, TAMU, is acknowledged for help editing the manuscript. Supplemental Table S2 contains data from an in-house, random, small-molecule library obtained from Dr. James Sacchettini's laboratory at Texas A&M University and Texas A&M AgriLife Research.
0.25% trypsin-EDTA | Gibco Invitrogen | 15050-065 | with phenol red |
0.4% trypan blue | MilliporeSigma | T8154 | liquid, sterile |
1.5 mL microcentrifuge tubes | Thermo Fisher | AM12400 | RNase-free Microfuge Tubes |
5 mL serological pipette | Corning | 29443-045 | Corning Costar Stripette individually wrapped |
10 mL serological pipette | Corning | 29443-047 | Corning Costar Stripette individually wrapped |
15 mL conical tubes | Falcon | 352196 | sterile |
20 µL filter tips | USA Scientifc Inc. | P1121 | sterile, barrier |
25 mL serological pipette | Corning | 29443-049 | Corning Costar Stripette individually wrapped |
50 mL conical tubes | Corning | 430828 | graduated, sterile |
150 mL auto-friendly reservior | Integra Bioscience | 6317 | sterile, individually wrapped for cell seeding in day 1 |
150 mL auto-friendly reservior | Integra Bioscience | 6318 | sterile, stacked, for loading dye in day 2 |
384/ 12.5 µL low retention tips | Integra Bioscience | 6405 | long, sterile filter |
384/ 12.5 µL tips | Integra Bioscience | 6404 | long, sterile filter |
384-well plate | Greiner | 781091 | CELLSTAR, clear polystyrene, µClear, Black/Flat |
Aluminum plate seals | Axygen Scientific | PCR-AS-200 | polyester-based |
Aluminum foil wrap | Walmart | ||
Biosafty cabinet II | NuAire | NU-540-300 | |
Cell counter | Nexcelom | AutoT4 | |
cell counting slides | Nexcelom | SD-100 | 20 µL chamber |
CO2 humidified incubator | Thermo Fisher | Forma Series II | |
Desk Lamp | SunvaleeyTEK | RS1000B | |
Dimethyl sulfoxide | MilliporeSigma | 276855 | anhydous, >99.9% |
Drug plate | Corning | 3680 | |
Dulbecco's phosphate-buffered saline | Corning | 21-031-CV | DPBS, 1x without calcium amd magnesium |
Ethanol | Koptec | 2000 | |
F-12K Nutrient Mixture | Corning | 45000-354 | (Kaighn's Mod.) with L-glutamine |
Fetal bovine serum | Equitech-Bio | SFBU30 | |
Fluorescent calcium assay kit | ENZO Lifescience | ENZ-51017 | 10×96 tests |
G418 sulfate | Gibco Invitrogen | 10131-027 | Geneticin selective antibiotic 50 mg/mL |
Hank's buffer | MilliporeSigma | 55037C | HBSS modified, with calcium, with magnesium, without phenol read |
HEPES buffer | Gibco Invitrogen | 15630-080 | 1 Molar |
HTS data storage plateform | CDD vault | https://www.collaborativedrug.com/ | |
Liquid handling system | Integra Bioscience | Viaflo | 384/12.5 µL |
Plate centrifuge | Thermo Fisher | Sorvall ST8 | |
Plate reader | BMG technology | Clariostar | |
Poly-D-lysine | MilliporeSigma | P6407 | |
Rhimi-K-1 agonist peptide | Genscript | custom order | QFSPWGamide |
T-75 flask | Falcon | 353136 |