Aspects of the invention related to identifying G-coupled receptors in insects such as mosquitoes and the design and use of screens that employ these proteins to identify efficacious insecticides.
Many neglected tropical infectious diseases affecting humans are transmitted by arthropods such as mosquitoes and ticks. New mode-of-action chemistries are urgently sought to enhance vector management practices in countries where arthropod-borne diseases are endemic, especially where vector populations have acquired widespread resistance to insecticides.
Some aspects of the invention include a method of controlling a population of invertebrates. The method comprises contacting an invertebrate with a compound that includes at least one of the group consisting of asenapine, amperozide, cis-(Z)-flupenthixol, benztropine, methiothepin, loxapine, mianserin, clomipramine, amperozide, butaclamol, clozapine, doxepin, and SCH23390 amitriptyline, chlorpromazine chlorprothixene. In some more particular embodiments, the invertebrate is an insect. In some more particular embodiments, the compound includes at least one of the group consisting of asenapine, amperozide, and cis-(z)-flupenthixol. In some more particular embodiments, the compound includes at least one of the group consisting of benztropine, methiothepin, loxapine, chlorprothixene, mianserin, and clomipramine. In some more particular embodiments, the insect is a mosquito. In some more particular embodiments, the compound includes at least one of the group consisting of asenapine, amperozide, cis-(Z)-flupenthixol, chlorpormazine, amitriptyline, and doxepin. In some more particular embodiments, the insect is a termite. In some more particular embodiments, the compound includes at least one of the group consisting of asenapine, cis-(Z)-flupenthixol, amperozide, amitriptyline, and chlorpromazine. In some more particular embodiments, the insect is a cockroach. in some more particular embodiments, the compound includes at least one of the group consisting of amitriptyline, chlorpormazine, cis-(Z)-flupenthixol, asenapine, and amperozide. In some more particular embodiments, the invertebrate is an arthropod. In some more particular embodiments, the invertebrate is a tick. In some more particular embodiments, the compound includes at least one of asenapine, chlorpromazine, and amperozide.
Some aspects of the invention include a method of controlling a population of invertebrates. The method comprises contacting an invertebrate with a compound that includes at least one of the group consisting of amitriptyline hydrochloride, (±)-butaclamol hydrochloride, clozapine, doxepin hydrochloride, cis-(Z)-flupenthixol dihydrochloride, methiothepin maleate, mianserin hydrochloride, niclosamide, piceatannol, and resveratrol. In some more particular embodiments, the invertebrate is any one of an arthropod, an insect, a mosquito, a termite, a cockroach, and a tick.
Some aspects of the invention include a method of controlling a population of invertebrates. The method comprises contacting an invertebrate with a compound that includes at least one of the group of chemical scaffolds consisting of dibenzocycloheptane derivatives, phenothiazine derivatives, thioxanthene derivatives, butyrophenone derivatives, diphenyl amine-containing compounds, quinazoline derivatives, benzodiazoxide derivatives, indole derivatives, piperazinylpyrazolopyrimidines, aryldiaminopryimidine, and hexahydrothienopyridines. In some more particular embodiments, the invertebrate is any one of an arthropod, an insect, a mosquito, a termite, a cockroach, and a tick.
Some aspects of the invention include a method of controlling a population of invertebrates. The method comprises contacting an invertebrate with a compound that includes at least 3 conjugated or non-conjugated rings, wherein said rings are independently selected from the group consisting of aromatic, non-aromatic, heterocyclic and non-heterocyclic rings, and wherein heterocyclic rings may independently includes at least one of the following atoms: C, O, N, or S and wherein each ring may be independently, substituted or unsubstituted with at least one of the following groups or atoms including H, amines, imines, ketones, aldehydes, alcohols, thiols, aromatic rings, alkanes, alkenes, alkynes, halogens and the like. In some more particular embodiments, the invertebrate is any one of an arthropod, an insect, a mosquito, a termite, a cockroach, and a tick.
Some aspects of the invention include a method of controlling a population of insects. The method includes contacting an insect with at least polynucleotide that interferes with the expression of a gene having at least 90 percent homology to SEQ ID. NO. 1 or SEQ ID. NO. 3. In some more particular embodiments, the polynucleotide interferes with the expression of at gene that includes or that hybridizes under stringent conditions to SEQ ID. NO. 1 or SEQ ID. NO. 3.
Some aspects of the invention include a method of screening. The method includes expressing a polynucleotide having at least 90 or 95 percent homology to SEQ ID. NO. 1 or SEQ ID. NO. 3, and contacting cells that express said polynucleotide with at least one exogenous compound, wherein said polynucleotide is not expressed in its native host. In some more particular embodiments, the polynucleotide has at least 90 or 95 percent identity to SEQ ID. NO. 1 or SEQ ID. NO. 3. In some more particular embodiments, the polynucleotide is SEQ ID. NO. 1 or SEQ ID. NO. 3.
Some aspects of the invention include a method of screening. The method includes an isolated polypeptide having at least 90 or 95 percent homology to SEQ ID. NO. 2 or SEQ ID. NO. 4; with at least one exogenous compound and detecting a interaction between said polypeptide and the at least one compound. In some more particular embodiments, the polypeptide has at least 90 or 95 percent identity to SEQ ID. NO. 2 or SEQ ID. NO. 4. In some more particular embodiments, the polypeptide is SEQ ID. NO. 2 or SEQ ID. NO. 4.
SEQ ID NO. 1. Gene Aedes aegypti dopamine receptor AaDOP1 Gene sequence.
SEQ. ID NO. 2. Aedes aegypti Dopamine Receptor AaDOP1 conceptual amino acid sequences.
SEQ. ID NO. 3. AADOP2 gene sequence.
SEQ. ID NO. 4. AADOP2 conceptual amino acid sequence.
SEQ ID NO. 5. Gene Ixodes scapularis dopamine receptor IsDOP1 Gene sequence.
SEQ. ID NO. 6. Ixodes scapularis Dopamine Receptor IsDOP1 conceptual amino acid sequences.
SEQ. ID NO. 7. IsDOP2 gene sequence.
SEQ. ID NO. 8. IsDOP2 conceptual amino acid sequence.
For the purposes of promoting an understanding of the principles of the novel technology, reference will now be made to the preferred embodiments thereof, and specific language will be used to describe the same. It will nevertheless be understood that no limitation of the scope of the novel technology is thereby intended, such alterations, modifications, and further applications of the principles of the novel technology being contemplated as would normally occur to one skilled in the art to which the novel technology relates are within the scope of this disclosure and the claims.
Mosquitoes (Class Insecta; Order Diptera; Family Culicidae) vector multiple neglected tropical diseases (NTDs) affecting human health, including malaria, yellow-fever, dengue and filariasis. Historically, insecticides employed against arthropod disease vectors have reduced the impact of NTDs, but unfortunately, continued disease control is threatened by the widespread development of vector populations that are resistant to insecticidal chemistries [1]. This issue is further complicated by the fact that there has not been a new public health insecticide produced for vector-borne disease control for over 30 years [2]. Recently, philanthropic investment has focused attention toward the development of new drugs to control NTDs in the human population [3]. It is widely recognized that an arsenal of new vector control solutions are required in order to meet this and other public health goals regarding NTDs. Thus, the research community should aggressively pursue the discovery of new mode-of-action chemistries for mosquito control through both traditional phenotypic screening and novel target-based approaches.
Novel insecticide targets potentially exist among the arthropod G-protein coupled receptors (GPCRs). These proteins comprise a large family of membrane-bound molecules that mediate critical biological processes such as neurotransmission, vision, and hormonal regulation, among others [4-5]. GPCRs are extensively targeted for drug development in humans—approximately 40% of prescription pharmaceuticals interact with these receptors [6]—and more recently, Gamo et al. [7] reported multiple GPCR-interacting chemistries as promising anti-malarial leads. Also, the mode-of-action of amitraz, a chemistry registered for tick and insect control, is presumed to have partial agonistic activity at an octopamine sensitive GPCR [8]. More than 100 different GPCRs have been identified in the genomes of multiple insect species, including malaria- and yellow fever-transmitting mosquitoes [9-10]. These studies have provided a basis for the functional characterization of GPCRs and their prioritization as potential subjects for insecticide development.
The biogenic amine-binding GPCRs (rhodopsin-like) are integral components of the central and peripheral nervous systems of eukaryotes and include receptors that bind the neurotransmitters dopamine, histamine, octopamine, serotonin, tyramine, and acetylcholine [11]. The dopamine receptors are classified as either D1- or D2-like [12] based on their differential functional roles. Ligand binding to the D1-like dopamine receptors causes Gαs-mediated stimulation of adenylyl cyclase (AC) production of cAMP. A reciprocal effect is observed following agonist activation of D2-like dopamine receptors, whereby cAMP production by AC is inhibited via Gαi/o proteins. Dopamine and its receptors are essential for complex behavioral mechanisms in arthropods such as locomotion [13-15], arousal [16], and olfactory learning [17-18].
The importance of dopaminergic-related functions has stimulated research to understand these processes in mosquitoes. Dopamine and serotonin have been tied to salivary gland functioning of vectors [19-20] and may have an impact on pathogen acquisition and transmission during blood feeding. Andersen et al. [21] reported that increased levels of dopamine were detected in Aedes aegypti following a blood meal that were implicated in ovarian or egg development, and in newly-emerged adults, presumably as part of the sclerotization process. Much attention has been given to the role of dopamine in the melanization pathway of mosquitoes and other insects, as well as the effect of dopamine on development, pigmentation, reproduction, immune responses to parasites, wound healing, and Wolbachia infection [22-27]. In the mosquito Culex pipiens, dose-dependent increases in cAMP were detected following treatment with dopamine and octopamine in homogenized head tissues, suggesting the presence of Gαs-coupled receptors that are responsive to these biogenic amines [28]. Putative D1-like and D2-like dopamine receptors have been identified in the genomes of the mosquitoes Ae. aegypti [9] and Anopheles gambiae [10], but research investigating their pharmacological properties is lacking. These genomic sequences provide a logical starting point to functionally characterize the receptors, which is needed to improve our comprehension of dopaminergic processes in mosquitoes. Moreover, due to their presumed significance in mosquito neurobiology, these dopamine receptors are attractive candidates to explore as new targets for chemical control.
Described herein is a “genome-to-lead” approach for insecticide discovery that uses perhaps the first reported chemical screen of a G-protein coupled receptor (GPCR) mined from a mosquito genome. A combination of molecular and pharmacological studies was used to functionally characterize two dopamine receptors (AaDOP1 and AaDOP2) from the yellow fever mosquito, Aedes aegypti. Sequence analyses indicated that these receptors are orthologous to arthropod D1-like (Gαs-coupled) receptors, but share less than 55% amino acid identity in conserved domains with mammalian dopamine receptors. Heterologous expression of AaDOP1 and AaDOP2 in HEK293 cells revealed dose-dependent responses to dopamine (EC50: AaDOP1=3.1±1.1 nM; AaDOP2=240±16 nM). Interestingly, only AaDOP1 exhibited sensitivity to epinephrine (EC50=5.8±1.5 nM) and norepinephrine (EC50=760±180 nM), while neither receptor was activated by other biogenic amines tested. Differential responses were observed between these receptors regarding their sensitivity to dopamine agonists and antagonists, level of maximal stimulation, and constitutive activity. Subsequently, a chemical library screen was implemented to discover lead chemistries active at AaDOP2. Fifty-one compounds were identified as “hits”, and follow-up validation assays confirmed the antagonistic effect of selected compounds at AaDOP2. In vitro comparison studies between AaDOP2 and the human D1 dopamine receptor (hD1) revealed markedly different pharmacological profiles and identified amitriptyline and doxepin as AaDOP2-selective compounds. In subsequent Ae. aegypti larval bioassays, significant mortality was observed for amitriptyline (93%) and doxepin (57%), confirming these chemistries as “leads” for insecticide discovery.
This research provides a pipeline useful for prioritization, pharmacological characterization, and expanded chemical screening of additional GPCRs in disease-vector arthropods. The differential molecular and pharmacological properties of the mosquito dopamine receptors highlight the potential for the identification of target-specific chemistries for vector-borne disease management, also reported are dopamine receptor antagonists with in vivo toxicity toward mosquitoes.
Presented herein are the results of a “proof-of-concept” study involving a “genome-to-lead” approach for developing new mode-of-action insecticides for arthropod disease vectors (
Materials and Methods—Aedes aegypti
The gene sequences for the putative dopamine receptors AaegGPRdop1 SEQ I.D. NO. 2 (AAEL003920) and AaegGPRdop2 (AAEL005834) (referred to hereafter as Aadop1 and Aadop2, respectively) in Ae. aegypti [10] were downloaded from VectorBase (http://www.vectorbase.org/index.php) [29]. Sequences of the D1-like dopamine receptors in Drosophila melanogaster were used to identify and compare conserved structural features [30-31].
Gene expression analyses for each receptor were conducted using RNA extracted from the eggs, larvae, pupae, and adult male and female mosquitoes from the Liverpool strain of Ae. aegypti [10]. Total RNA was isolated using TRIzol Reagent (Invitrogen, Carlsbad, Calif.) and then treated with RNase-Free DNase (QIAGEN, Valencia, Calif.). The SuperScript One-Step RT-PCR kit (Invitrogen, Carlsbad, Calif.) was used to amplify receptor mRNA from approximately 150 ng total RNA per reaction using the primers and experimental conditions provided in supporting information, Table 5. RT-PCR amplification products were electrophoresed and compared by size to the DNA HyperLadder I (Bioline USA Inc., Randolph, Mass.). Products were cut from the gel and isolated with the Qiagen Gel Extraction Kit (Qiagen Valencia, Calif.). The cloning procedure was performed using the TOPO TA cloning kit (Invitrogen, Carlsbad, Calif.), according to the manufacturer's instructions. DNA sequencing was conducted at the Purdue University Genomics Core Facility. The resultant DNA sequences were used to predict full-length coding regions that were manually annotated using Artemis software (version 9) [32].
A neighbor joining sequence analysis was performed using the deduced amino acid sequences representing the mosquito dopamine receptor proteins (referred to hereafter as AaDOP1 and AaDOP2), additional representative biogenic amine receptors from the insects D. melanogaster and A. mellifera, and the human D1- and D2-like dopamine receptors. ClustalW 1.83 [33] was used for sequence alignments prior to tree construction in PAUP 4.0b4a [34]. The bootstrap method (100 replicates) was used to provide branch support. Alignments of amino acid sequences for determination of conserved motifs were conducted using Multalin software [35]. Conserved amino acid residues and additional protein features were predicted as described by Meyer et al. [36].
Functional characterization of AaDOP1 and AaDOP2 was conducted by heterologous expression in HEK293 cells (ATCC, Manassas, Va.) [36]. Expression constructs were produced by synthesis (GenScript, Piscataway, N.J.) and included the partial Kozak transcriptional recognition sequence “CACC” added directly upstream of the transcription initiation codon for each gene. Constructs were cloned into pUC57 and then subcloned into the expression vector pcDNA3.1+ (Invitrogen, Carlsbad, Calif.) by GenScript (Piscataway, N.J.). Stable cell lines co-expressing either AaDOP1 or AaDOP2 with a CRELuc reporter construct were developed to permit pharmacological studies in a 384-well format [36-37]. Briefly, cells already stably expressing the CRELuc reporter construct were transfected in a 10 cm dish with 15 μl Lipofectamine-2000 and 3 μg of pcDNA3.1+/Aadop1 or pcDNA3.1+/Aadop2. Clones were maintained as described for the wild-type HEK293 cells [36] with the addition of 2 μg/ml puromycin and 300 μg/ml Geneticin (Sigma-Aldrich, St. Louis, Mo.).
For initial functional analysis, the receptors were transiently expressed in HEK293 cells [36] and analyzed using a competitive binding assay to measure levels of cAMP accumulation [37]. Dose-response curves were generated using cells stably expressing the receptors [36-37]. The compounds used for pharmacological characterization included dopamine hydrochloride, histamine dihydrochloride, 5-hydroxytryptamine hydrochloride (serotonin), (±)-octopamine hydrochloride, and tyramine hydrochloride (Sigma-Aldrich, St. Louis, Mo.) and (−)-epinephrine bitartrate and L (−)-norepinephrine bitartrate (Research Biochemical International, Natick, Mass.). The synthetic dopamine receptor ligands tested included SKF38393 and SKF81297 (Tocris, Ellisville, Mo.), SCH23390 (Tocris, Ellisville, Mo.), and dihydrexidine (DHX). Data was collected from a minimum of three independent replicate experiments with each sample measured in triplicate. Statistical analysis of data was conducted with GraphPad Prism 5 software (GraphPad Software Inc., San Diego, Calif.).
In order to identify novel AaDOP2 receptor antagonists, a Library of Pharmacologically Active Compounds (LOPAC1280) was screened, using HEK-CRELuc-Aadop2 cells. These cells were cultured as described above, expanded, and cryo-preserved, to produce a uniform cell population. Briefly, cells (˜2.5×107) were harvested by non-enzymatic dissociation [0.5 mM EDTA in Ca2+Mg2+ free-phosphate buffered saline (CMF-PBS)] resuspended in cell culture media, and pelleted by centrifugation for 5 min at 100×G. The pellet was resuspended in freezing media (Opti-MEM supplemented with 10% DMSO and 10% FBS) to a concentration of 5×106/ml, frozen step-wise, and held in liquid N2 until use. Cells were rapidly thawed, diluted in Opti-MEM, and 20 μl containing 25,000 cells were plated per well in 384-well plates (Nunc, Fisher Scientific, Pittsburgh, Pa.) using a BiomekFX liquid handling station (Beckman-Coulter, Brea, Calif.). The plates were incubated overnight in a humidified incubator at 37° C. and 5% CO2.
Prior to screen initiation, a “checkerboard” analysis was conducted that included a minimum (300 nM dopamine in combination with 10 μM SCH23390) and maximum (300 nM dopamine) stimulatory condition. The data obtained were analyzed to calculate the Z-factor [38] using a modified equation that accounts for the number of replicates (NIH website: http://www.ncgc.nih.gov/guidance/section2.html#summary-calcs).
All compounds were diluted to appropriate concentrations and suspended in assay buffer (Opti-MEM supplemented with 0.02% ascorbic acid) using a BiomekFX 96-tip head. All LOPAC1280 compounds were screened in quadruplicate at a concentration of 10 μM, including duplicate samples on two separate assay plates in different quadrants to control for plate and automation effects. Each plate contained a dopamine response curve (14 nM-30 μM) and antagonist control wells (10 μM SCH23390 in combination with 300 nM dopamine). Following compound addition, dopamine was added to each test well at a final concentration of 300 nM, and cells were incubated for 2 hr at 37° C. in a humidified incubator. The plates were then equilibrated at 25° C. prior to the addition of Steadylite plus luminescence reagent (PerkinElmer, Shelton, Conn.). Plates were incubated on a shaker at 300 rpm for 5 min, and the luminescence signal was measured using a DTX880 multimode reader (Beckman Coulter, Brea, Calif.) with a 1 sec integration time.
Raw screen data were processed as follows: the average background luminescence (cells in the absence of dopamine or LOPAC1280 compound) was subtracted from the raw data. Values for the positive receptor activation control (300 nM dopamine) were averaged within each assay plate and used to establish a 100% dopamine receptor stimulation level. Similarly, the average response to SCH23390 was calculated within each assay plate to establish a baseline inhibition for antagonist chemistries. The average percent compound effect was calculated for each LOPAC chemistry in comparison to the SCH23390 antagonist control. The minimum criterion for selection of an antagonist “hit” was established as the percent inhibition equivalent to that determined for SCH23390+3 standard deviations.
Subsequent validation assays using both the AaDOP2 and the human D1 dopamine receptor (hD1) (Sunahara et al., 1990) were conducted for select identified chemistries using a competitive binding cAMP accumulation assay. In addition to SCH23390, these included amitriptyline hydrochloride, doxepin hydrochloride, niclosamide, clozapine, (+)-butaclamol hydrochloride, cis-(Z)-flupenthixol dihydrochloride, resveratrol, mianserin hydrochloride (Sigma, St. Louis, Mo.), piceatannol and methiothepin maleate (Tocris, Ellisville, Mo.). The drugs were suspended from dimethyl sulfoxide (DMSO) stocks in Hanks Balanced Salt Solution (HBSS) (HyClone, Logan, Utah) with 0.1% fatty acid free bovine serum albumin (BSA) and 20 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES), and serial dilutions were prepared using a Precision 2000 automated pipetting system (BioTek, Winooski, Vt.). The cAMP accumulation assay was carried out as previously described [36-37] with minor modifications to permit processing of a larger number of samples in a semi-automated fashion. Briefly, AaDOP2- or hD1-expressing cells were harvested using Hank's based non-enzymatic cell dissociation reagent (Invitrogen, Carlsbad, Calif.), triturated in equal parts Dulbecco's modified eagle medium (DMEM) (Invitrogen, Carlsbad, Calif.) dissociation reagent, centrifuged 5 min at 100×G, and resuspended in HBSS supplemented with 0.1% BSA and 20 mM HEPES. Cells were seeded (50,000 cells in 40 μl) in clear 96-well plates and incubated at 37° C. with 5% CO2 for 1 hr. The cAMP accumulation assay was carried out in HBSS supplemented with final concentrations of 0.1% BSA, 20 mM HEPES, 0.5 mM 3-isobutyl-1-methylxanthine (IBMX), and 0.02% ascorbic acid in a final volume of 50 μl. The selected compounds were added to the wells in duplicate, followed by addition of dopamine (final concentration 3 μM for AaDOP2 and 100 nM for hD1). Plates were incubated at room temperature for 1 hr, and the assay was terminated by addition of 25 μl of 9% ice-cold trichloroacetic acid (TCA). Cell lysates were incubated on ice for at least 1 hr prior to quantifying cAMP accumulation as previously described [36-37].
In Vivo Ae. aegypti Bioassays
Single dose-point and dose response in vivo mosquito bioassays were used to assess the toxicity of selected AaDop2 receptor antagonists identified in the chemical screen. Larvae of Ae. aegypti (Liverpool strain) were reared under standard laboratory conditions on a 12 hr day/night cycle at 75% RH and 28° C., and bioassays were conducted at room temperature (22-24° C.). Larvae were transferred from standard rearing trays into six-well tissue culture plates (Corning, Inc. Corning, N.Y.) using a small plastic pipette. Ten L4-stage larvae were included per well, each containing five ml of de-ionized water and an assigned drug concentration. Controls were conducted similarly but lacked a drug treatment. Bioassays employed a double-blind experimental design, and larval mortality was scored 24 hr following administration of drugs. Single dose-point assays were conducted using 400 μM drug and included three biological replicates each consisting of 100 larvae. Dose-response assays were conducted using five doses (400, 200, 100, 50 and 25 μM. Five technical replicates, each including 10 larvae, were performed per dose, and the assay was repeated three times. Statistical analysis included one sample t-tests (single-point assays) and determination of the LC50 and LC90 values (dose-response assays) and were conducted with GraphPad Prism 5 software (GraphPad Software Inc., San Diego, Calif.).
mRNA transcripts for Aadop1 and Aadop2 were detected by RT-PCR from the eggs, larvae, pupae, and adult male and female Ae. aegypti (
The deduced amino acid sequences of AaDOP1 and AaDOP2 were analyzed to identify conserved structural features typically associated with biogenic amine-binding GPCRs (supporting information, Table 6), as well as unique regions that could be potentially exploited for development of mosquito-specific chemistries. Conserved features included sites predicted for ligand binding, protein stability, G-protein coupling, and post-translational modification. Alignments of the full-length AaDOP1 and AaDOP2 amino acid sequences (
In order to study the function of the putative dopamine receptors AaDOP1-IEOID No. 3 and AaDOP2-IEOID No. 4, each receptor was expressed in HEK293 cells. Production of the mosquito receptor transcripts in transiently-transfected cells was first verified using RT-PCR (supporting information,
Subsequently, dose-response curves were generated using cells stably expressing either AaDOP1 or AaDOP2 for seven different biogenic amines (
The AaDOP2 receptor for an antagonist screen of the LOPAC1280 library because of its low constitutive activity and strong dopamine response compared to background (approximately 10-fold) (
The criterion for “hit” detection was established relative to the control antagonist (SCH23390 response+3 standard deviations), such that only those compounds that inhibited the dopamine response by at least 81% were considered (Table 2). Based on this, the screen identified 51 potential antagonists of the AaDOP2 receptor (additional screen results provided in Table 8). These compounds were partitioned into seven different classes based on their known biochemical interactions with mammalian molecular targets that included dopamine receptor antagonists (20), serotonin (6), histamine (2), and acetylcholine receptor ligands (1), biogenic amine uptake inhibitors (9), protein kinase modulators (6), and miscellaneous chemistries such as cell cycle regulators and apoptosis inhibitors (7).
Ten “hit” compounds (amitriptyline hydrochloride, (±)-butaclamol hydrochloride, clozapine, doxepin hydrochloride, cis-(Z)-flupenthixol dihydrochloride, methiothepin maleate, mianserin hydrochloride, niclosamide, piceatannol, and resveratrol), in addition to SCH23390 were selected for screen validation assays. These compounds were tested for their activity in cAMP accumulation assays to control for potential “off-target” effects (i.e. chemistries that affect the CRELuc reporter system). Seven of these compounds were potent antagonists of the AaDOP2 receptor, as shown by the dose-dependent reduction of cAMP accumulation relative to the dopamine-stimulated control (Table 3). Three of the compounds (i.e. niclosamide, piceatannol, and resveratrol) showed no significant antagonistic effects against AaDOP2 in the cAMP accumulation experiments, having IC50 values ≧10 μM.
The active compounds were also tested against hD1 to allow for comparisons of relative potency between species (Table 3). These comparisons clearly indicated a unique pharmacology between AaDOP2 and hD1 with divergent rank order functional potencies that showed no significant correlation (R2<0.15). For example, the prototypical mammalian D1 antagonist, SCH23390 was greater than 3000-fold more selective for hD1 than AaDOP2. In contrast, the data also revealed that two structurally-related tricyclic antidepressants (i.e. amitriptyline and doxepin) had more than 30-fold selectivity for AaDOP2 when compared to hD1. These observations suggest that the significant differences between these receptors could be exploited for the development of AaDOP2-selective compounds.
In Vivo Ae. aegypti Bioassays
The toxicity of the AaDOP2 antagonist screen hits amitriptyline and doxepin was assessed in Ae. aegypti larval bioassays. These chemistries were selected due to their relatively higher potency at AaDOP2 compared to hD1 (Table 3). Single dose point assays at 400 μM effective concentration of drug revealed that amitriptyline (93% average mortality) and doxepin (57% average mortality) each caused significant larval mortality (p<0.05) 24 hours post-treatment relative to the water control (0% mortality), whereas no mortality was observed for SCH23390 during this timeframe (data not shown). Dose response experiments were conducted for amitriptyline, which was prioritized due to the rapid and high mortality effect observed in the single-point assays. The toxicity of amitriptyline was dose-dependent, and the LC50 and LC90 values were determined at 78 μM and 185 μM, respectively (
This work provides the first detailed investigation into the molecular and pharmacological properties of D1-like dopamine receptors, AaDOP1 and AaDOP2, from the mosquito vector of dengue and yellow fever, Ae. aegypti, and the development of a cell-based screen assay to discover antagonists of AaDOP2. Our study employed a novel pipeline utilizing a “genome-to-lead” approach for the discovery of new chemistries for vector control. This research establishes a basis for improving understanding of mosquito dopaminergic processes in vivo and for chemical screening of these and other receptors characterized in arthropod vectors of human disease, such as in the Lyme disease tick, I. scapularis [36,41]. To our knowledge, Lee and Pietrantonio [46] have published the only other study involving the functional characterization of a biogenic amine-binding GPCR in mosquitoes that was focused on a Gαs-coupled serotonin receptor in Ae. aegypti. Furthermore, ligands of only four other cloned GPCRs have been pharmacologically verified in mosquitoes, including those that target an adipokinetic hormone receptor, a corazonin receptor, a crustacean cardioactive peptide receptor [47], and an adipokinetic/corazonin-related peptide receptor in the malaria mosquito, A. gambiae [48].
Typically, insects possess three different dopamine receptors including two D1-like receptors and a single D2-like receptor [42]. Here, RT-PCR data were used to validate the two mosquito D1-like dopamine receptor gene models [10]; this enabled confirmation of intron/exon boundaries and prediction of the complete protein coding regions needed prior to heterologous expression studies. A putative D2-like dopamine receptor gene (AaDOP3) was also identified in Ae. aegypti [10] although this receptor has not yet been functionally characterized. The RT-PCR studies also demonstrated that transcripts for both D1-like dopamine receptor genes were detectable in each developmental stage of Ae. aegypti, suggesting the importance of these receptors throughout the mosquito life cycle. Much progress has been made in determining the life-stage and tissue-specific expression dynamics of the orthologous dopamine receptors in D. melanogaster [14,30-31,39,43,49], A. mellifera [40,42,50-53], and most recently in the Lyme disease tick, I. scapularis [41]. Our research will support future complementary studies needed to localize expression of these dopamine receptors in mosquito tissues to gain further insight toward their neurophysiological roles.
The AaDOP1 and AaDOP2 amino acid sequences were compared and analyzed to identify conserved as well as unique features of the receptors. Several characteristics typically associated with biogenic amine-binding GPCRs were evident, including aspartate residues in TM II as well as TM III that are thought to interact with the amine moieties of catecholamines [54]. The conserved serine residues in TM V and aromatic residues in TM V and VI are also potentially important for ligand interaction [55-56]. In both receptors, the conceptual cytoplasmic region of TM III contained the conserved “DRY” motif associated with G-protein coupling [57-58], and a pair of cysteine residues were located in the extracellular loops I and II that may form a disulfide bond for protein stabilization [57,59-60]. Interestingly, the divergent intracellular loop III was predicted to be almost twice as long in AaDOP2 (115 amino acids) than in AaDOP1 (62 amino acids), but the sizes of the carboxyl tail region were similar between these receptors. This corresponded well with the relative sizes of these features in the fruit fly and honeybee orthologs [42]; however, the significance of these characteristics is yet to be determined in the mosquito. Importantly, the AaDOP1 and AaDOP2 sequences were markedly different from the human D1-like dopamine receptor sequences. Although a modest level of amino acid identity (˜50%) was observed between the TM domains, the N- and C-termini and extracellular and intracellular loop regions were highly divergent (data not shown). These differences suggest that there exists potential for identifying chemistries that are mosquito-specific and, importantly, do not interfere with dopaminergic functioning in humans.
Heterologous expression experiments conducted in HEK293 cells provided experimental evidence that the Ae. aegypti receptors are functional D1-like dopamine receptors. We measured significant increases in cAMP accumulation following dopamine treatment of cells transiently expressing either AaDOP1 or AaDOP2, suggesting that both receptors are coupled to Gαs proteins. This effect was further substantiated in cell lines stably co-expressing either of these receptors and the CRELuc reporter system, as measured by an increase in luciferase activity following dopamine treatment. Future research is needed to determine if these receptors operate through multiple cellular signaling mechanisms, such as was shown for the D. melanogaster dopamine receptor involved with both cAMP and calcium signaling [61].
The stably transformed cell lines were used to compare the pharmacological properties of AaDOP1 and AaDOP2 in response to seven different biogenic amines. For dopamine, we measured EC50 values in the nanomolar range for both AaDOP1 (3.1±1.1 nM) and AaDOP2 (240±16 nM). However, there were differences in the responses of these receptors to the other biogenic amines. AaDOP2 was activated only with dopamine, whereas AaDOP1 was stimulated by dopamine, epinephrine and to a lesser extent, norepinephrine. These results were similar to those reported for the orthologous dopamine receptors in the tick I. scapularis [36,41]. Another difference between AaDOP1 and AaDOP2 was observed regarding constitutive activity. In both transient and stable expression experiments, the AaDOP1 receptor exhibited significant constitutive activity, as determined by the elevated levels of cAMP detected in the absence of a receptor agonist, whereas AaDOP2 did not. Such constitutive activity was also reported for the D1-like dopamine receptors AmDOP1 of A. mellifera [40], CeDOP1 from the nematode Caenorhabditis elegans [62], Isdop1 of I. scapularis [36], and the human D5 receptor [63]. Seifert and Wenzel-Seifert [64] proposed that constitutive activity of a GPCR may enable the maintenance of basal neuronal activity, although evidence is needed to support such activity for AaDOP1 in vivo.
The pharmacological properties of AaDOP1 and AaDOP2 were further explored by testing their responses to synthetic dopamine receptor agonists and antagonists. Both receptors were strongly stimulated by the agonist DHX; however, only AaDOP1 significantly responded to the well characterized D1 agonists SKF81297 and SKF38393. This differential response to the SKF compounds was also observed for the orthologous D1-like dopamine receptors in the tick I. scapularis [36]. Interestingly, neither of the D. melanogaster D1-like dopamine receptors was strongly stimulated by SKF38393 [31,39]. Both AaDOP1 and AaDOP2 were inhibited by the antagonist SCH23390, as were the tick D1-like receptors [36]. This contrasted with the lack of significant inhibition reported by SCH23390 for D-dop1 in the fruit fly [39] and DOP1 of the honeybee [50]. Given the limited number of drugs that have been tested against these receptors, to date, these differential pharmacological responses provide further evidence that it may be possible to discover chemistries that operate specifically at the mosquito dopamine receptors.
Our over-arching goal was to develop a pipeline to identify lead chemistries active at biogenic-amine binding GPCRs in vector arthropods. Broadly speaking, we define a lead chemistry as any molecule, or its analog or derivative, with potential for insecticide development. In our study, this refers to any molecule identified by screening and subsequently confirmed in a variety of “hit-to-lead” assays. The LOPAC1280 library was chosen for our pilot screen because it is enriched with chemistries that influence dopaminergic processes and includes other GPCR-binding ligands. We hypothesized that chemistries that antagonize these dopamine receptors may possess insecticidal properties. Precedent for this notion stems from pest management successes associated with the use of Fipronil and cyclodienes, which antagonize GABA-gated chloride channels and have highly insecticidal properties (Bloomquist, 1996; Casida and Quistad, 1998). This notion was pursued using HEK293 cells stably expressing AaDOP2 because this receptor has a robust response to dopamine and a low constitutive activity, which are properties that aid interpretation of screen data. Our initial screen was directed at the identification of AaDOP2 antagonists; the success of this experiment justifies expanded screening to explore the antagonist chemical “space”, and with assay modification, screens to detect agonists active at this receptor. Moreover, development of the AaDOP1 assay would enable comparative screens against LOPAC1280 chemistries.
Of the 51 hit AaDOP2 antagonists identified in the LOPAC1280 library, 20 (39%) are known antagonists of mammalian dopamine receptors. A majority of these chemistries fall into the benzodiazepine, phenothiazine or thioxanthene classes that in other systems are known to bind other biogenic amine receptors. Included were ligands selective for D1- and D2-like dopamine receptors in mammalian systems, as well as several non-dopamine receptor selective compounds such as (±)-butaclamol, cis-(Z)-flupenthixol, and the atypical antipsychotic, clozapine. These three compounds were tested in a dose-response format for their ability to inhibit dopamine-stimulated cAMP accumulation. The IC50 values demonstrated the following rank order of potency clozapine>cis-flupenthixol>butaclamol. The next largest grouping of identified compounds includes inhibitors of the biogenic amine transporters (9 compounds, 18%). Several serotonin receptor antagonists (6 compounds, 12%) were identified as well. Follow-up dose response studies with selected chemistries from the identified transport inhibitors and serotonin antagonists (i.e. methiothepin, mianserin, amitriptyiline, and doxepin) revealed that these compounds were potent antagonists at the AaDOP2 receptor and were much more potent than the prototypical D1 antagonist, SCH23390 (Table 3). The antagonistic activity of these ligands is not completely surprising; the National Institute of Mental Health's Psychoactive Drug Screening Program (NIMH-PDSP) database reports Ki values for the human D1-like dopamine receptors at 80-900 nM (http://pdsp.med.unc.edu/). However, these observations, combined with the dopamine antagonist screen results, indicate that well studied and clinically used compounds could be used to target invertebrate GPCRs. In fact, a number of the chemistries identified in our screen have been used in humans for decades, suggesting the possibility of “drug repurposing” as insecticides. Further precedent for the concept of insect-specific chemistries can be drawn from the fact that a number of insecticides (e.g., pyrethroids and fipronil) are considerably more selective at invertebrate as opposed to mammalian targets [65]. The screen also identified multiple protein kinase modulators and several agents that regulate germane cellular functions that presumably inhibit the CRE response via non-AaDOP2 mechanisms. Support for this hypothesis was demonstrated in the direct measurement of cAMP accumulation experiments, where resveratrol, pieacetannol, and niclosamide each lacked activity. The remaining three “hit” compound classes included antagonists of either histamine or muscarinic acetylcholine receptors, and this likely reflects the lack of receptor selectivity for these ligands.
The LOPAC1280 library includes several known antagonists of mammalian dopamine receptors that did not qualify as hits in our screen. In part, this can be explained by the fact that we used a highly stringent cut-off to signify antagonistic activity at AaDOP2. Had we reduced the stringency to select for hits with an antagonistic effect equivalent to that of SCH23390+6 standard deviations (69% inhibition), our screen would have returned an additional 13 hit chemistries, including compounds predicted to have a modest antagonistic effect at AaDOP2 and those that are more selective for D2-like dopamine receptors. Considering the substantial divergence between the mosquito and human D1-like dopamine receptor sequences, there is a strong possibility that a subset of the “non-hit” dopamine receptor antagonists are not active at the mosquito receptor. In support of this, the prototypical mammalian D1 antagonist, SCH23390, was greater than 3000-fold more selective for hD1 than AaDOP2. Although our comparison data set is limited to only eight compounds, these experiments suggest a very divergent pharmacology between these human and mosquito dopamine receptors. Thus, our study provides a foundation for subsequent comparative pharmacological analyses of the mosquito and human dopamine receptors.
We observed a strong correlation between in vitro and in vivo data. The AaDOP2 antagonist screen hits, amitriptyline and doxepin, caused significant lethality in the mosquito larval assay. Our finding that these drugs each have a relatively higher potency at the mosquito dopamine receptor than at hD1 has implications for the identification of arthropod-selective chemistries. Drugs with minimal or no impact on the neurological functioning of humans or other vertebrate species are particularly desirable as prospects for insecticide development. Conversely, SCH23390, which is active at AaDOP2 only in the micromolar range and was several fold more selective for hD1 in cAMP assays, did not cause significant larval mortality at 24 hr.
The success of this initial chemical library screen in identifying new mosquitocidal chemical leads justifies the pursuit of an expanded high-throughput screening effort involving thousands or hundreds of thousands of chemistries against mosquito dopamine receptors. Our platform is also amenable for the screening of agonist chemistries active at these mosquito dopamine receptors, as well as for Gαs-coupled biogenic amine targets of other vector arthropods, and also could be modified to screen Gαi/o-coupled receptors [66]. Importantly, the identification of lead AaDOP2 receptor antagonistic chemistries provides a basis for investigating the effect of these or related compounds on mosquito dopaminergic processes in vivo [67]. Follow-up research is needed to determine the precise mechanism(s) of amitriptyline- and doxepin-induced mortality in A. aegypti larvae. Further work is also needed to determine if these chemistries and associated derivatives or analogs identified by chemical screens possess the properties desired of an insecticide (e.g. bioavailability, in vivo potency/toxicity, suitable half-life, lack of effects on non-target organisms, suitability for synthesis and formulation). Molecular modeling of three dimensional GPCR structures and their binding capabilities, as reported for an adipokinetic hormone receptor in A. gambiae [68] and a tyramine receptor in the moth Plodia interpunctella [69], may facilitate in silico chemical screening [70] and ligand-receptor studies that permit the design or refinement of lead molecules active at mosquito GPCRs.
Historically, multiple neuroactive processes in arthropods have been exploited for pest control using insecticides such as chlorinated hydrocarbons, organophosphates, methylcarbamates, pyrethroids, amidines, and phenylpyrazoles. Resistance involving each of these classes (the vast majority of which operate by affecting ion channels and neurotransmitters) has been documented. The development of new mode-of-action insecticides could improve our arsenal against mosquito populations that have developed resistance to existing chemical formulations [1]. We suggest that the two dopamine receptors characterized here, as well as other biogenic amine-binding GPCRs [72-73], represent promising targets for new insecticide research, due to their presumably central roles in insect neurobiology. This “proof-of-concept” study sets the stage for target-specific approaches for vector control. Such efforts, in parallel with activities of organizations such as the Innovative Vector Control Consortium, may help to realize the goal of delivering new insecticides for reduction of vector-borne diseases [2].
Development of a structure-activity relationship (SAR) for AaDOP2 antagonist molecules has been pursued on two main fronts: (1) cell-based SAR studies to identify novel chemistries that antagonize the AaDOP2 target in vitro, and (2) complementary whole-organism in vivo SAR studies to identify those chemistries that cause mortality of key insect pests and have greatest potential for insecticide development.
cAMP Accumulation Confirmation Assays.
Confirmation of screen hits has been carried out using a secondary cellular assay and evaluation of in vitro potency of antagonist molecules (i.e., the inhibitory concentration or IC50 value). An in vitro SAR analysis of hits from the LOPAC1280 screen for antagonists of the AaDOP2 receptor using a Homogenous Time Resolved Fluorescence (HTRF) cAMP accumulation assay (Cisbio Dynamic 2) was completed. These assays were performed as described in Meyer [84] with modifications to miniaturize the assay to 384-well format. We evaluated the potency of more than 40 potential antagonist chemistries at the AaDOP2 target (Table 9). Chemistries were tested in dose-response experiments in the presence of 304 dopamine in AaDOP2 expressing cells in comparison to the prototypical mammalian dopamine antagonist SCH23390 and amitriptyline. The data in Table 9 are representative of three independent experiments performed in duplicate.
To explore more unique chemical space and identify highly novel chemistries that could be pursued as leads for insecticide development, an additional 6,067 compounds in our Cisbio high-throughput in vitro screen were screened. Three libraries were screened in total; these were the NIH Clinical Collection and Clinical Collection 2 (NCC and NCC 2; 727 small molecules with a history of use in human clinical trials), the Spectrum Collection (MicroSource Discovery Systems Inc., Gaylordsville, Conn.; 2,320 compounds; 60% drugs, 25% natural products and 15% other bioactive molecules, including insecticides), and the TimTec Natural Derivatives Library NDL-3000 (TimTec LLC, Newark, Del.; 3,040 compounds; synthetic compounds and synthetically modified pure compounds). Screens were performed as per methods described in Meyer [81] and Ejendal [78] with modifications and appropriate validation procedures to execute the screening assays in 384 well format using the Cisbio HTRF cAMP assay.
In Vivo Aedes aegypti High-Throughput Screen.
As a second step, the in vivo activity of chemistries identified via in vitro studies was evaluated. Initially, chemistries were evaluated in a high-throughput in vivo assay against L3 stage larvae of the yellow fever mosquito, Ae. aegypti. The chemistries were tested in duplicate at a single point dose of 40004 in 24-well plate (BD Bioscience, San Jose, Calif.) assay. Briefly, chemistries were re-suspended in water and added to wells containing five Ae. aegypti L3 larvae to achieve a final concentration of 400 μM per well in 1 ml sterile RO water. The plates were incubated at 22° C. and the assay was scored at 30 minutes and 1, 1.5, 2, 2.5, 3, 24, 48 and 72 hours. The screen was performed as a double-blind experiment and enables a rapid evaluation of the in vivo activity of multiple water-soluble chemistries, and the prioritization of molecules that exhibit rapid and/or high mosquito mortality. The toxicity of 25 AaDOP2 antagonist chemistries using this assay is presented in for three independent screens (n=3) in Tables 8 and 10 and
Ae. aegypti Concentration Response Curves and Time Course Assays.
Chemistries that cause rapid and high mortality of mosquitoes in our high-throughput in vivo screen were subsequently evaluated in concentration response assays to determine LC50 (lethal concentration) and LT50 (lethal time) values as described by Meyer et al [81], with the results presented in Tables 11 and 15. Briefly, chemistries were diluted to appropriate concentration in water and added to the wells of 24-well plates containing 5 L3 stage Ae. aegypti larvae to achieve a final concentration of 400, 200, 100, 50 and 25 μM in 1 ml sterile RO water. Chemistries were tested in quadruplicate at each concentration. The plates are incubated at 22° C. and percent mortality was assessed at 30 minutes and 1, 1.5, 2, 2.5, 3, 24, 48 and 72 hours. Results reported here represent a minimum of three independent CRC experiments (i.e., n=3).
A series of whole organism bioassays was developed to further evaluate the toxicity of chemistries to a range of agricultural, veterinary and public health pests and non-target insects (Tables 12 and 13). These assays are designed to rapidly explore the biological activity spectrum of chemistries as well as to facilitate initial investigations of the mode of action of insecticidal compounds. Briefly, assays were developed for the German cockroach (Blatella germanica, urban pest), termite (Reticulitermes flavipes; urban pest), lone star tick (Amblyomma americana; veterinary and public health pest), soybean aphid (Aphis glycines, agricultural pest) and honeybee (Apis mellifera; beneficial, non-target insect). The assays were designed to test the toxicity and speed of kill of multiple chemistries on contact. Each of these assays was performed as a single point dose experiment in triplicate using either topical application of compound to the ventral thorax (1 μl of a 200 μM compound solution in 1:1 v:v DMSO:EtOH vehicle; 50-80 μg effective dose depending on compound) as is the case for our cockroach and honeybee assays, or 10 second immersion in a solution of the compound diluted in water or DMSO:EtOH to a final concentration of 100 μM. The assays were conducted in 24-well plate format with mortality assessment following incubation at 22° C. as described above. Commercial insecticides were incorporated in these assays to provide a positive control and permit initial comparative analyses. Lead chemistries identified in our in vivo screen have been evaluated in the cockroach and termite assays.
cAMP Accumulation Confirmation Assays.
An in vitro SAR for antagonists of the AaDOP2 receptor using a Homogenous Time Resolved Fluorescence (HTRF) cAMP accumulation CRC assay (Cisbio Dynamic 2) has been completed. The potency (IC50 value) of more than 40 potential antagonist chemistries at the AaDOP2 target (Table 9) was evaluated, identifying fourteen chemistries with high potency (IC50<100 nM), 12 chemistries with moderate potency (IC50 between 100 nM-1 μM), and eight chemistries with low potency (IC50>1 μM). Six antidepressant molecules (e.g., fluoxetine) and the insecticide molecule, Amitraz, showed no inhibition in the assay. From this preliminary SAR, we have identified a total of five general chemical scaffolds (i.e., dibenzocycloheptane derivatives, phenothiazine derivatives, thioxanthene derivatives, butyrophenone derivatives, and diphenyl amine-containing compounds) that provide the basis for an expanded SAR. It is also noteworthy that within several of the chemical scaffolds, are multiple distinct and more specific chemical classes. For example, the dibenzocycloheptane derivatives also contain dibenzazepines and dibenzodiazepines, and the diphenylamines contain three unique diphenyl structures (i.e., diphenylpiperazine, diphenylpiperidine, and diphenylmethoxy derivatives).
The initial SAR analysis in Table 9 provides the following information on ligand requirements: The amine state contributes to the determination of antagonist potency at the AaDOP2 receptor. For example, compounds with tertiary amines (clomipramine, imipramine, amitriptyline, and loxapine) are approximately five to 100 fold more potent than the secondary amine analogs of these compounds (norclomipramine, desipramine, nortiptyline, and amoxapine, respectively). Furthermore, clomipramine is approximately six fold more potent than imipramine, suggesting that aromatic ring substituents can enhance the activity of the identified antagonists. It should also be noted that the data suggest that the central ring is not necessary for antagonist activity, as risperidone, benztropine and amperozide are moderately potent. However, compounds with six or seven-membered central rings are generally the most potent AaDOP2 receptor antagonists, suggesting a relationship between conformational rigidity and antagonist activity. Interestingly, a conformationally rigid molecule, asenapine, is one of the two most potent compounds tested to date. The structural features noted above can be used to assist in guiding future studies.
In total, 93 new hit compounds were identified in our small molecule screens of the AaDOP2 target against the NIH Clinical Collection, Spectrum Collection and TimTec NDL-3000 library (Tables 9 and 14). Specifically, forty-five hit compounds (including 24 new hits) in the NIH Clinical Collections that caused ≧80% inhibition of the AaDOP2 receptor. Seventeen of these hits were previously identified in a screen of the LOPAC1280 library and four additional compounds had been evaluated in follow-up CRC assays due to their structural similarity to LOPAC1280 library hits. 4 antagonist compounds (58 new hits) were identified from the Spectrum Collection that caused ≧80% inhibition of the AaDOP2 target. Thirty-six of these compounds were previously identified in the LOPAC1280 screen. 11 compounds were identified from the TimTec NDL-3000 library that caused ≧60% inhibition of the AaDOP2 target. Within these 11 compounds at least three new scaffolds have been identified, including quinazoline, benzodiazoxide and indole derivatives.
Quantities of a selection of hit chemistries for follow-up in vitro CRC confirmation assays were obtained to confirm screen hits and evaluate the comparative potency (IC50 values) of compounds. Additionally, a screen of 10,000 chemistries from the ChemDiv library (San Diego, Calif.) was completed. An additional 62 hits were identified from this screen and the top 28 of these were obtained for follow-up confirmation assays. The results of the initial 10,000 compound screen were subjected to a chemoinformatics analysis. The results of analysis revealed a number of promising trends enabling the development of a preliminary molecular model of the AaDOP2 receptor. The results also identified three novel putative chemical scaffolds for AaDOP2 (i.e. piperazinylpyrazolopyrimidines, aryldiaminopyrimidines, and hexahydrothienopyridines).
Thus a total of 17,367 compounds have been screened at the AaDOP2 target, and more than 10 chemical scaffolds and 150 unique hits via the “hit-to-lead” process [81] have been pursued.
In Vivo Aedes aegypti SAR Screen.
Twenty-five chemistries were evaluated for toxicity to Ae. aegypti larvae in a high-throughput in vivo assay (Tables 8 and 10;
These 10 chemistries caused higher mortality at 24 hours and were faster acting than amitriptyline. Three chemistries (asenapine, chlorpromazine, and amperozide) were extremely fast acting and caused greater than 70% mortality of the mosquito population within 30 minutes and five chemistries (cis-flupenthixol, chlorprothixene, mianserin, loxapine and methiothepin) caused greater than 70% mortality within three hours. These 10 chemistries were selected for in vivo concentration response analyses and additional insecticide bioassays (see Sections 2.2 and 2.3). Additionally, five chemistries with moderate mosquito toxicity (i.e., 40-70% mortality at 24 hours post exposure) and nine chemistries with limited or no toxicity to mosquito larvae (i.e., 0-40% mortality at 24 hours) were identified. The in vivo SAR data showed a strong correlation with the in vitro SAR results (see Tables 9 and 10), suggesting a possible mechanism for the chemistries to act by disrupting the AaDOP2 target in vivo. However, amperozide, which caused rapid and high mosquito mortality, is a weak inhibitor of AaDOP2 in vitro (570±110 nM), suggesting that other mechanisms may also be involved.
Ae. aegypti Concentration Response Curve and Time Course Assays
The lethal concentration (LC50; Table 11,
Several chemistries (amitriptyline, chlorpromazine and cis-(z)-flupenthixol) that cause moderate toxicity (25-35% mortality) to cockroaches on contact within three to five days post dose (Table 12) have been identified. Using a published feeding assay [80] that involves incorporation of test compounds into an inert bait matrix, no significant oral toxicity of lead chemistries to cockroaches (Table 16) was demonstrated, suggesting that mode of action is associated with absorption via the cuticle. Importantly, four lead chemistries (cis-(z)-flupenthixol, chlorpromazine, amitriptyline and amperozide) were shown to cause rapid toxicity to adult termites within 24 hours, with significant mortality (70-100%) of the test population within five days (Table 13). Although not wishing to be bound by any theory, one possibility is that the smaller biomass and thinner, less melanized cuticle of R. flavipes in comparison to B. germanica significantly increases the bioavailability of test compound and contributes to the higher mortality we observe in termite assays. The insect bioassay results expand the known activity spectrum of our chemistries to three insect orders, namely Diptera (flies and mosquitoes), Blattodea (cockroaches) and Isoptera (termites), thus suggesting significant commercial potential of resultant insecticides
Materials and Methods—Ixodex scapularis.
The intronless genes IscaGPRdop1 and IscaGPRdop2 (referred to henceforth as Isdop1 and Isdop2) correspond to the automated gene models ISCW001496 and ISCW008775 from the I. scapularis genome assembly, respectively, which were downloaded from VectorBase (http://www.vectorbase.org/index.php) (Lawson et al., 2009). These genes were identified with blastn searches (Altschul et al., 1997) using the D. melanogaster dopamine receptor sequences D-Dop1 (Gotzes et al., 1994) and DopR99B (DAMB) (Feng et al., 1996; Han et al., 1996) as queries against the I. scapularis genome assembly. Isdop1 encodes a protein of 425 amino acids and is located on scaffold DS648196 of the genome assembly at positions 133,404e134,681. Isdop2 is 457 amino acids in length and located on scaffold DS812273, spanning from base pair positions 248,624e247,251. Gene annotations were managed with Artemis software (http://www.sanger.ac.uk/resources/software/artemis/) (Rutherford et al., 2000).
A neighbor joining phylogenetic analysis was conducted to determine the relationships of the full-length deduced amino acid sequences for Isdop1 and Isdop2 with multiple biogenic amine receptors from the insects D. melanogaster and A. mellifera, and the human dopamine receptors. Sequence alignments were performed using ClustalW 1.83 (Chema et al., 2003), and tree construction was conducted with PAUP 4.0b4a (Swofford, 2001). Support for branches in the tree was generated with the bootstrap method (100 replicates) in PAUP.
Multalin software (Corpet, 1988) was used to align conceptual amino acid sequences of the dopamine receptors from I. scapularis, D. melanogaster, A. mellifera and Homo sapiens for comparative analyses of key structural components. To provide support for identification of putative TM domains, hydrophobicity plots were generated with ProtScale software (Kyte and Doolittle, 1982) available at the ExPASy Proteomics Server, Swiss Institute of Bioinformatics (http://ca.expasy.org/tools/protscale.html). Kinase-specific protein phosphorylation sites for protein kinase A and protein kinase C were predicted using the NetPhosK 1.0 Server, Technical University of Denmark (http://www.cbs.dtu.dk/services/NetPhosK/). Putative 1-4 N-linked glycosylation sites were identified with the NetNGlyc 1.0 Server, Technical University of Denmark (http://www.cbs.dtu.dk/services/NetNGlyc/) and EnsembleGly (http://turing.cslastate.edu/EnsembleGly/) (Caragea et al., 2007). Putative palmitoylation sites were identified with CSS-Palm 2.0 (http://csspalm.biocuckoo.org/online.php) using a medium threshold (Ren et al., 2008).
Blastn searches (Altschul et al., 1997) were conducted at Gen-Bank and VectorBase to search for expressed sequence tags (ESTs) corresponding to the Isdop1 and Isdop2 sequences. For gene expression, a pooled sample of total RNA was isolated from 10 adult female I. scapularis with TRIzol Reagent (Invitrogen, Carlsbad, Calif.). Total RNA was then treated with RNase-Free DNAse (QIAGEN, Valencia, Calif.) prior to cDNA synthesis. Because both Isdop1 and Isdop2 are single exon genes, it was not possible to design primers that would span intronic regions as a control for contaminating genomic DNA. To circumvent this issue, 30-rapid amplification of complementary DNA ends (30-RACE) experiments were designed according to the directions provided in the GeneRacer_Kit (Invitrogen, Carlsbad, Calif.). Target amplification of synthesized cDNA was conducted using two gene-specific forward primers (nested PCR) and reverse primers (provided in kit) which annealed to the polyadenylated 30 tail. Gene-specific primers used for amplification of the 30 end of Isdop1 included Ixodes_Dop1—30RACE—1F (50-CCTATCGCACCAAGAGAAGCACCATTTG-30) and Ixodes_Dop1—30 RACE—2F (50-GGATGCCGAAGCAACAACAACACTG-30), respectively. Gene-specific primers used for amplification of the 30 end of Isdop2 included Ixodes_Dop2—30RACE—1F (50-TGGATCAACTCCGGCATGAAC CCCATCA-30) and Ixodes_Dop2—30RACE—2F (50-ACCGCCTCCGACGCATCATCAAGGAAGA-30). To amplify cDNA produced in the 30-RACE procedure, a standard PCR protocol was used according to the reagent manufacturers' instructions (Bioline USA Inc., Randolph, Mass.). Parameters used for cDNA amplification in both the initial and nested PCR steps included a denaturation step of 94_C for 2 min followed by a modified touchdown PCR including (i) 5 cycles of denaturation at 94_C for 30 s, primer annealing at 70_C and extension at 70_C for 30 s, (ii) 5 cycles of denaturation at 94_C for 30 s, primer annealing at 65_C and extension at 68_C for 30 s, and (iii) 25 cycles of denaturation at 94_C for 30 s, primer annealing at 60_C and extension at 68_C for 30 s. A final extension period of 10 min was conducted at 68_C. Amplification products were separated on 1% TBE gels and compared by size to DNA Hyper-Ladder I (Bioline USA Inc., Randolph, Mass.), and amplicons of interest were excised and extracted from the gel using the Qiagen Gel Extraction Kit (Qiagen Valencia, Calif.). Purified DNA was cloned with the pCR 2.1-TOPO TA cloning kit using One Shot_Top 10 F′ chemically-competent Escherichia coli cells (Invitrogen, Carlsbad, Calif.). Following transformation, cells were grown in 10 cm petri dishes containing Luria Bertani (LB)-agar media+50 μg/ml kanamycin+10 ml X-Gal/IPTG solution (Bioline USA Inc., Randolph, Mass.) to permit blue/white colony selection. Randomly-selected clones were grown in 5 ml LB media+50 μg/ml kanamycin for 16-18 h at 37_C. Plasmid DNA was isolated using the QIAprep spin miniprep kit (Qiagen, Valencia, Calif.), and DNA sequencing of four independent clones for each cloned amplicon was conducted at the Purdue University Genomics Core Facility. A consensus sequence for each cloned amplicon was produced by first aligning the sequences with Multalin software (Corpet, 1988) and then using a majority-rule criterion at each base pair position to filter out potential sequencing errors or polymorphisms. Consensus sequences were compared to their corresponding regions in the I. scapularis genome sequence assembly for calculations of percent nucleotide identities.
Clones corresponding to the coding regions of Isdop1 and Isdop2 were produced by synthesis (GenScript, Piscataway, N.J.) according to their gene models (see above). The partial Kozak transcriptional recognition sequence “CACC” was synthesized directly upstream of the transcription initiation codon for each gene. Each gene was first cloned into the vector pUC57 and then subcloned into the expression vector pcDNA3.1+ (Invitrogen, Carlsbad, Calif.) by Gen-Script (Piscataway, N.J.).
To functionally characterize Isdop1 and Isdop2, heterologous expression experiments were conducted in HEK 293 cells (ATCC, Manassas, Va.) using transient transfections as follows. HEK 293 cells were cultured in Dulbecco's modified Eagle Medium (DMEM) (Sigma, St. Louis, Mo.) supplemented with 5% bovine calf serum (BCS), 5% fetal clone I serum (Hyclone, Logan, Utah), and 1% AntieAnti (Invitrogen, Carlsbad, Calif.) in a humidified incubator at 37_C with 5% CO2. HEK 293 cells were seeded in 48 well cluster plates (BD Falcon, Franklin Way, N.J.) and transiently transfected with 200 ng of the Isdop1 or Isdop2 construct DNA, or with the vector pcDNA3.1+ alone in the controls, and 0.5 ml Lipofectamine-2000 (Invitrogen, Carlsbad, Calif.) per well. Forty-eight hours post-transfection, the media was decanted and replaced with assay buffer (Earle's balanced salt solution supplemented with 2% BCS, 0.02% ascorbic acid, 15 mM HEPES and 0.5 mM 3-isobutyl-1-methylxanthine) containing each drug (see below) at the concentrations indicated in Table 3 and
The drugs used for pharmacological studies included dopamine hydrochloride, histamine dihydrochloride, 5-hydroxytryptamine hydrochloride, (_) octopamine hydrochloride, and tyramine hydrochloride (SigmaeAldrich, St. Louis, Mo.) and (_)-epinephrine bitartrate and L (_)-norepinephrine bitartrate (Research Biochemical International, Natick, Mass.). To block epinephrine and norepinephrine activation of endogenous b adrenergic receptors in HEK cells, propranolol (Research Biochemical International, Natick, Mass.) was included in experiments at 300 nM (results summarized in Table 19). The synthetic dopamine receptor agonists that were tested included SKF38393 and SKF81297, and the antagonists were SCH23390 and (+)-butaclamol (SigmaeAldrich, St. Louis, Mo.). At least three independent experiments were conducted for each pharmacological assay.
To expand on our initial findings measuring cAMP responses in transiently-expressing cells and facilitate pharmacological studies more amenable to a high-throughput format, we developed stable cell lines co-expressing either Isdop1 or Isdop2 and a cyclic AMP Response Element/Luciferase (CRE-Luc) reporter system, as previously described by Przybyla et al. (2009). To test this system for each I. scapularis receptor, we conducted two different reporter assays consisting of dose-response curves for dopamine, epinephrine, and norepinephrine and an assay involving the antagonist SCH23390 (10 μM) in combination with dopamine (1 μM). Cells previously transfected with the CRE-Luc reporter construct (HEK CRE-Luc cells) were transfected in a 10 cm dish with 15 μl Lipofectamine-2000 and 3 μg of either the Isdop1 or Isdop2 constructs described above. Clones were selected in Geneticin (SigmaeAldrich, St. Louis, Mo.) and screened for receptor function in response to dopamine. The resulting clones were maintained as described for the wild-type HEK 293 cells above, with the addition of puromycin (2 μg/ml) and Geneticin (300 μg/ml). For the receptor reporter assays, approximately 20,000 cells were seeded per well in white clear-bottomed 384-well plates (Nunc, Fisher Scientific, Pittsburgh, Pa.) and incubated overnight. The cell media was decanted and replaced with 50 μl Opti-MEM (Invitrogen, Carlsbad, Calif.) supplemented with 0.02% ascorbic acid containing increasing concentrations of the indicated biogenic amine. After 2 h incubation in a humidified incubator at 37_C with 5% CO2, the assay buffer was decanted, 20 μl PBS and 15 μl Steadylite HTS (PerkinElmer, Waltham, Mass.) were added per well, and the plate was incubated on a shaker for 5 min at room temperature. Luminescence was detected on a VictorLight (PerkinElmer, Shelton, Conn.) as counts per second (CPS). Each condition was carried out in at least triplicate. The GraphPad Prism software (GraphPad Software Inc., San Diego, Calif.) was used for statistical analysis of pharmacology data. Statistical analyses included either a Student's t-test or a one-way ANOVA followed by Dunnett's post-hoc test as indicated in the figure legends. p values <0.05 were considered statistically significant.
Results for Ixodes scapularis
A neighbor joining analysis was conducted to investigate the relationships of the deduced amino acid sequences of each putative I. scapularis dopamine receptor with multiple biogenic amine receptors from D. melanogaster, A. mellifera and the human D1- and D2-like dopamine receptors (
Initially, the seven TM domains of Isdop1 and Isdop2 were identified through comparisons to the TM domain sequences from their putative orthologs in D. melanogaster (Gotzes et al., 1994; Feng et al., 1996) (
The deduced amino acid sequences for Isdop1, SEQ. ID NO. 6, and Isdop2, SEQ. ID NO. 8, were analyzed to identify structural features typically associated with GPCRs and for sites predicted to undergo post-translational modification (Table 18). Each receptor contains a cysteine residue in the first and second extracellular loops, which are believed to form disulfide bonds to maintain GPCR stability (Dixon et al., 1987; Karnik et al., 1988; Fraser, 1989). There are two and four cysteines positioned near the C-terminus of Isdop1 and Isdop2, respectively, which may be suitable sites for post-translational palmitoylation and potentially contribute to the establishment of a fourth intracellular loop (Jin et al., 1999). In the N-terminal extracellular domains, two and four 1e4 N-linked glycosylation sites were predicted in Isdop1 and Isdop2, respectively, which may be needed for localization in the plasma membrane (Karpa et al., 1999). Multiple potential protein kinase A and protein kinase C phosphorylation sites were identified in the intracellular loops and carboxyl termini of each receptor (Namkung and Sibley, 2004). Isdop1 has a short third intracellular loop (47 amino acids), relative to this larger feature in Isdop2 (90 amino acids).
Each I. scapularis sequence was further examined to identify key amino acid residues which are conserved among biogenic amine receptors in the rhodopsin-like receptor subfamily of GPCRs (Table 18). Both receptor sequences contained the signature aspartate (D) residues in TMII and TMIII which are involved with binding the amine groups of catecholamines (Strader et al., 1988). Adjacent to the terminus of TMIII, Isdop2 contained the conserved “DRY” motif believed to function in coupling with G proteins (Dixon et al., 1987; Fraser et al., 1988), whereas Isdop1 included a substitution of phenylalanine (F) for tyrosine (Y) in the third position of this motif. Three conserved serines were present in TMV of each receptor that are reportedly important for hydrogen bonding with catechol hydroxyl groups (Strader et al., 1989; Pollock et al., 1992). Conserved aromatic residues were identified in TMV and TMVI of each receptor which are thought to interact with the catechol aromatic ring during ligand binding (Strader et al., 1995).
We were unable to identify ESTs in GenBank corresponding to Isdop1 or Isdop2, which necessitated gene expression analyses to determine if mRNA transcripts were produced for these receptors in vivo. We were successful in detecting transcripts for both of the single exon genes Isdop1 and Isdop2 using a 30-RACE method, as analyzed by gel electrophoresis. For Isdop1, a single 30-RACE amplification product (1050 bp) was produced and cloned. A consensus sequence was generated that, after the reverse primer sequence was removed, had >99% nucleotide identity (998/1002 bp) to the 30 end of the ISCW001496 gene sequence from the I. scapularis genome assembly. This sequence included a 107 bp region upstream of the final position in the predicted stop codon and then extended 894 bp downstream into the presumed 30 untranslated region (UTR). A putative polyadenylation site (AATAAA) for Isdop1 was identified in the 30-RACE amplified region 862 bp downstream of the final base pair position of the predicted stop codon.
For Isdop2, two prevalent (704 bp and 497 bp) and three less abundant (w1.2 kb, w350 bp, w250 bp) 30-RACE products were produced. For purposes of confirming gene expression, we successfully cloned and sequenced the two most abundant products. The consensus sequences of these amplicons had, after the reverse primer sequences were removed, >99% (649/652 bp) and >99% (448/451 bp) nucleotide identities to the corresponding 30 end of the ISCW008775 gene sequence from the I. scapularis genome assembly. The 50 region of these amplification products corresponded to the same position where the forward primers were designed, 70 bp upstream from the final base pair position of the predicted stop codon. Thereafter, these overlapping sequences differed in length as they extended into the presumed 30-UTR region 380 bp and 580 bp downstream of the final base pair position of the predicted stop codon, respectively. A putative polyadenylation site (AATAAA) was identified in the 30-RACE amplified region 358 bp downstream of the final base pair position in the predicted stop codon.
To investigate the pharmacological characteristics of Isdop1 and Isdop2, each cloned receptor was transiently-expressed in HEK 293 cells, and responses to drug treatments were measured in terms of cAMP accumulation. Dopamine treatment (10 μM) of cells expressing these receptors resulted in a significant increase of intracellular cAMP (
To further characterize the pharmacological profiles of Isdop1 and Isdop2, treatments with receptor agonists and antagonists were conducted (
Clones stably expressing Isdop1 and Isdop2 were developed and used in conjunction with a CRE-Luc reporter assay to facilitate functional activity measurements of these receptors in a 384-well format (
Next, we investigated the effect of co-treatment with the D1-like dopamine receptor antagonist SCH23390 in the CRE-Luc reporter system. Relative to stimulation with dopamine alone, co-treatment with SCH23390 caused a 74% and 69% inhibition of luciferase activity for Isdop1 and Isdop2, respectively (
The biogenic amine receptors represent candidate neuroactive targets for the identification of GPCR-specific chemistries potentially useful for tick control (Lees and Bowman, 2007). This study provided the first pharmacological characterization of two cloned dopamine receptors in the Lyme disease vector, I. scapularis. Previous research indicated that the dopaminergic pathway is a critical component of the salivary secretion mechanisms in ticks, which are stimulated during their bloodfeeding behavior (Sauer et al., 2000; Bowman and Sauer, 2004; Lees and Bowman, 2007). Thus, improvements in our knowledge of these complex processes may open new doors toward understanding the dynamics of blood-borne pathogen acquisition and transmission among species in this medically-important arthropod lineage.
Tick genome sequencing projects provide a means to expedite determinations of gene function (Van Zee et al., 2007). The two D1-like dopamine receptor sequences analyzed here, along with a single uncharacterized D2-like receptor identified in the genome assembly, are believed to compose the primary suite of receptors which regulate dopaminergic processes in I. scapularis. Interestingly, Isdop1 and Isdop2 are single exon genes, which contrasted with that reported for other invertebrate D1-like dopamine receptors that contain introns (Mustard et al., 2005), but corresponded with the intronless human D1-like dopamine receptors (Sunahara et al., 1990, 1991). The chromosomal positions of Isdop1 and Isdop2 are unknown because, to date, chromosomal mapping in I. scapularis includes only a preliminary linkage map (Ullmann et al., 2002) and a fluorescent in situ hybridization-based map of the major tandem repeats in the genome (Meyer et al., 2010). The 30-RACE experiments demonstrated that both Isdop1 and Isdop2 were expressed in adult female ticks, supporting their functionality in vivo. Further analyses of gene expression involving different life stages and specific tissues, as well as supporting investigations of the 50- and 30-untranslated regions, would improve our understanding of tick dopamine receptor biology. During the review process of this manuscript, a relevant study by _Simo et al. (2011) provided evidence that an I. scapularis D1 receptor (a.k.a. Isdop1) was positioned on the luminal surface of cells lining the salivary gland acini and activated by a paracrine signal of dopamine.
The previously characterized D1-like dopamine receptors in the model insects D. melanogaster and A. mellifera were used to predict orthologous relationships in I. scapularis and to identify conserved structural features of the tick receptors. A neighbor joining analysis of the Isdop1 and Isdop2 amino acid sequences separated the tick receptors into two different clades, with each Glade containing an orthologous receptor from the model insects used in the analysis. Each of these clades was separate from the Glade containing the two human D1-like dopamine receptors. Furthermore, this analysis showed that the Glade containing sequences for Isdop2 and its insect orthologs also included the octopamine receptors OAMB (Han et al., 1998) and AmOA1 (Oar) (Grohmann et al., 2003) from D. melanogaster and A. mellifera, respectively, which further distinguished these receptors from the other D1-like dopamine receptors used in the sequence alignment (Mustard et al., 2005). Our notion that Isdop1 and Isdop2 would function as D1-like dopamine receptors was supported by the fact that each of the orthologous receptors in the model insects were previously characterized as functional D1-like dopamine receptors (reviewed by Mustard et al., 2005). Based on alignments of the conserved TM domains, only 47-53% amino acid identities were shared between either of the tick receptors and the human D1-like dopamine receptors, whereas 65-74% identities were observed between the tick receptors and their putative orthologs in the model insects. Together, the sequence analyses provided evidence to suggest that the I. scapularis proteins would likely have different pharmacological characteristics from each other, as well as from the human and orthologous insect D1-like dopamine receptors.
The pharmacological properties of Isdop1 and Isdop2 were compared through heterologous expression experiments in HEK 293 cells. Both Isdop1 and Isdop2 were stimulated with dopamine treatment, as determined by the subsequent accumulation of intracellular cAMP; however, each receptor had a unique pharmacological profile. We detected (
Additional pharmacological differences between Isdop1 and Isdop2 were suggested in studies examining the specificity of these receptors in response to seven different biogenic amines, where Isdop1 was activated by dopamine, epinephrine, and norepinephrine, but Isdop2 only significantly responded to dopamine treatment. The high constitutive activity of Isdop1 made determining the effect of these biogenic amines somewhat challenging, whereas for Isdop2, dopamine was clearly the most robust stimulatory ligand among the biogenic amines tested. The significant response of cells transiently expressing Isdop1 to epinephrine and norepinephrine presumably reflects the lack of absolute receptor specificity to these biogenic amines at the upper ranges of concentrations tested (i.e. 10 μM). Subsequent dose-response studies revealed that the EC50 values were markedly higher for epinephrine and norepinephrine compared to dopamine for both Isdop1 and Isdop2. In terms of determining the suitability of these receptors for future high-throughput screening assays for agonists, it appears that Isdop2 is a better candidate than Isdop1 because of its strong response to stimulation and lack of detectable constitutive activity.
The I. scapularis receptors showed differential activities in their responses to D1-like dopamine receptor agonists, where Isdop1 was significantly more responsive than Isdop2. In spite of the predicted orthology of Isdop1 with the receptors in D. melanogaster and A. mellifera, divergences in their pharmacological profiles were observed with regard to agonist treatments. For example, stimulation of D-Dop1 from D. melanogaster with SKF38393 was 3-fold less than that observed for dopamine treatment, and SKF81927 also had relatively poor stimulatory effects (Gotzes et al., 1994; Sugamori et al., 1995). This also corresponded with that reported for DOP1 of A. mellifera, where SKF38393 was not an effective agonist (Blenau et al., 1998). In contrast to Isdop1, similar trends regarding agonist activities were observed between Isdop2 and its insect orthologs. For example, SKF38393 was only a weak agonist of DopR99B (DAMB) of D. melanogaster (Feng et al., 1996) and also failed to significantly stimulate DOP2 of A. mellifera (Mustard et al., 2003).
Inhibition of the dopamine response with antagonist treatment of cells transiently expressing either Isdop1 or Isdop2 also provided evidence highlighting the distinct pharmacological properties of these receptors and permitted functional comparisons between their orthologs in D. melanogaster and A. mellifera. SCH23390 was an effective antagonist of Isdop1, whereas butaclamol did not appear to have a considerable inhibitory effect. This was in contrast to that reported for the orthologous D-Dop1 receptor in D. melanogaster and the DOP1 receptor in A. mellifera, where receptor activities were significantly reduced with butaclamol, while SCH23390 was less effective (Gotzes et al., 1994; Blenau et al., 1998). For Isdop2, both butaclamol and SCH23390 had antagonistic effects. This corresponded with that reported for DopR99B (DAMB) of D. melanogaster, where modest inhibition was shown for both SCH23390 and butaclamol (Han et al., 1996; Feng et al., 1996). Similarly, each of these antagonists also significantly reduced cAMP levels in cells expressing DOP2 of A. mellifera, relative to the response observed following treatment with dopamine (Mustard et al., 2003). The antagonistic effect of SCH23390 observed in the stable cell lines expressing Isdop1 or Isdop2 was consistent with that in our transient expression experiments.
Further pharmacological studies of the two I. scapularis D1-like dopamine receptors, as well as the uncharacterized D2-like dopamine receptor, are warranted to better understand the basic properties of these receptors. In addition to stimulation of cAMP production by adenylyl cyclase in vivo, it will also be important to determine if the tick dopamine receptors activate Ca2 signaling pathways, as was shown for DopR99B (DAMB) of D. melanogaster (Reale et al., 1997; Feng et al., 1996). Our production of stable cell lines expressing Isdop1 and Isdop2 provided the foundation for ongoing efforts to screen chemical libraries to further characterize the pharmacology of these receptors. The significant inhibition of luciferase activity observed following treatment with the antagonist SCH23390 in cells stably expressing either Isdop1 or Isdop2 demonstrated the utility of the reporter system to readily measure the effect of antagonistic chemistries in a high-throughput format. In addition, inhibition of the constitutive activity of Isdop1 by SCH23390 indicated that this system may also be useful in screening for inverse agonists. Additional research is needed to determine if chemistries can be designed that specifically target these tick receptors but do not influence the dopaminergic processes in other organisms. This line of research in I. scapularis may serve as a model for complementary studies in other ticks for which chemical control is a primary means of population suppression for tick-borne disease management
Receptor activation was detected using a CRELuc luciferase reporter assay in which the luciferase (Luc) reporter gene (pGL3, Promega, Madison, Wis.) is under transcriptional control of five copies of the cAMP response element (CRE). Dopamine-stimulated activation of the receptor results in an accumulation of intracellular Camp leading to an increase in luminescence measured in counts per second (CPS). To study a uniform population of cells expressing the receptors, stable cell lines in HEK293 cells were established and maintained as previously described (Meyer et al., 2011). Cells were trypsinized, resuspended in cell culture media, pelleted by centrifugation for 5 min at 100 g, and resuspended in Opti-MEM media (Invitrogen, Carlsbad, Calif.) before being transferred to white 384-well plates (Nunc, Thermo Fischer Scientific, Rochester, N.Y.) using amultichannel pipette. Cells were plated at a final cell concentration of approximately 35,000 cells in 20 ml per well and incubated overnight in a humidified cell culture incubator (37_C with 5% CO2). CRELuc assays (for
Screening of the LOPAC1280 library (Sigma, St. Louis, Mo.) for Isdop2 receptor antagonists was executed at the Integrated Screening Technologies Laboratory, Discovery Park at Purdue University, West Lafayette, Ind. Approximately 35,000 Isdop2-expressing cells (20 ml/well) were plated in solid white 384-well plates (Nunc, Thermo Fisher Scientific, Rochester, N.Y.) using a BiomekFX instrument (Beckman-Coulter, Brea, Calif.), incubated overnight and the screening assays were carried out as previously described (Meyer et al., 2012) with the exception that fresh cells (in contrast to cryopreserved cells) were used. Briefly, library compounds were added to the cells followed by addition of dopamine, and the assay plates were incubated for 2 h (37_C with 5% CO2). The plates were equilibrated to room temperature, and the reaction was terminated following addition of 10 μl/well Steadylite Plus reagent (PerkinElmer, Waltham, Mass.). The dopamine stimulated response was determined by measuring luminescence on a DTX880 multimode reader (Beckman Coulter, Brea, Calif.) with a 1 s integration time. To evaluate the robustness of the screening protocol, the Z factor (Zhang et al., 1999) was determined for the assay platform as follows. Plates of Isdop2-expressing cells were prepared with alternating wells containing either dopamine (3 μM) or dopamine (3 μM) in combination with the dopamine receptor antagonist SCH23390 (10 μM) in a “checkerboard” pattern. Data from the assay plates were analyzed to calculate the Z0 using a modification of the original equation (Zhang et al., 1999) that considers the number of replicate samples in the screening assay protocol [http://assay.nih.gov/assay/index.php/Section2:Plate_Uniformity_and_Signal_Variability_Assessment]. Library plates were thawed, and aliquots were diluted in assay buffer using a BiomekFX liquid handling station (Beckman-Coulter, Brea, Calif.). Each library compound was screened in quadruplicate, with duplicates of each compound dispensed on separate plates as a control for potential plate-specific effects. Compounds were tested at a final concentration of 10 μM in combination with dopamine (3 μM). Each plate included samples of dopamine (3 μM) as a positive control for receptor activation, and dopamine (3 μM) in combination with SCH23390 (10 μM) as a control for antagonist chemistries (Meyer et al., 2011). A dopamine concentration-response curve (0.14-300 μM DA) was also included for each plate to enable comparison of data between plates. To analyze screen data, the average background luminescence for each plate was subtracted from all values on the corresponding plate. Wells containing dopamine only (no test compound) were used to establish the dopaminestimulated response control value for each plate. The percent inhibition of the dopamine-stimulated response associated with each test compound was calculated by averaging the four replicate values for each compound and normalizing them compared to the response associated with dopamine alone. These values were then compared to the inhibitory effect of the antagonist control (SCH23390) averaged between all 16 plates. Compounds were considered “hits” if their % inhibition value was three standard deviations below the % inhibition value recorded for SCH23390. In the present study the average_S.D. inhibitory value of SCH23390 was calculated as 86±6%. This value minus three standard deviations (i.e. 86%−(3×6%) ¼ 68%) was used to establish the cutoff for identification of hit compounds.
Screen validation analysis was conducted with cells stably expressing Isdop1 or Isdop2 (Meyer et al., 2011) using a cAMP accumulation assay as previously described (Meyer et al., 2012) with minor modifications. Briefly, cryopreserved cells were thawed and resuspended in HBSS-based assay buffer (Hank's Balanced Salt Solution (Hyclone, Logan, Utah) supplemented with 0.1% bovine serum albumin (BSA) (Sigma, St. Louis, Mo.), 20 mM 4-(2-Hydroxyethyl) piperazine-1-ethanesulfonic acid N-(2-Hydroxyethyl)piperazine-N0-(2-ethanesulfonic acid) (HEPES, Sigma, St. Louis, Mo.) and 0.02% ascorbic acid), counted, and seeded (30 ml/well) in clear tissue culture treated 96 well plates (BD Falcon, Franklin Lakes, N.J.). To obtain optimal signals in the cAMP accumulation assay, 50,000 cells (Isdop1) or 100,000 cells (Isdop2) were seeded and allowed to recover for 6 h (Isdop1) or 1 h (Isdop2) in a cell culture incubator before proceeding with the assay. Compounds evaluated in confirmation assays were amitriptyline hydrochloride, ()-butaclamol hydrochloride, cis-(Z)-flupenthixol dihydrochloride, clozapine, doxepin hydrochloride, mianserin hydrochloride (Sigma, St. Louis, Mo.), and methiothepin maleate and R(+)—SCH23390 (R(+)-7-Chloro-8-hydroxy-3-methyl-1-phenyl-2,3,4,5-tetrahydro-1H-3-benzazepine) hydrochloride (Tocris, Ellisville, Mo.). The antagonist hit compounds were prepared as DMSO stocks and then diluted in HBSS assay buffer supplemented with 3-isobutyl-1-methylxanthine (IBMX, final assay concentration 0.5 mM, Sigma, St. Louis, Mo.). Serial dilutions were carried out with a Precision 2000 liquid handling station (BioTek, Winooski, Vt.). Diluted antagonists were added and plates were incubated at room temperature for 10 min followed by addition of dopamine at concentrations corresponding to the approximate EC90 values for Isdop1 (10 nM) and Isdop2 (30 μM) as determined in the cAMP accumulation assay. Assays were carried out in duplicate and data analysis was conducted with GraphPad Prism v.5 software, constraining the top of the curves to 100 for both receptors and the Hill slope to —0.7 for Isdop2.
Two Gas-coupled dopamine receptors (Isdop1 and Isdop2) with distinct molecular and pharmacological profiles have been identified in the Lyme disease tick, I. scapularis (Meyer et al., 2011; _Simo et al., 2011). Due to their potential as new molecular targets for tick control, we aimed to further characterize and identify novel small molecule antagonists of Isdop1 and Isdop2. Our data in HEK293 cells support the notion that both Isdop1 and Isdop2 are Gas coupled, however, our observations do not exclude that these receptors may also signal through other pathways, which were not examined here. The Isdop2 receptor displays a lower constitutive activity and higher fold-stimulation in response to dopamine than Isdop1 (approximately 9-fold vs. 2-fold) when heterologously expressed in HEK293 cells (
Because antagonists of the AaDOP2 receptor showed toxicity to mosquito larvae (Meyer et al., 2012), we reasoned that the tick dopamine receptors may represent a new class of drug targets for acaricide development. To screen for novel antagonists of the Ixodes dopamine receptors, we developed a luciferase reporter-based screening platform using cells that stably express the Isdop2 receptor in conjunction with the CRELuc luciferase reporter construct (Meyer et al., 2011) to measure dopamine receptor activation. The CRELuc reporter assay is based on an indirect measure of cAMP; receptor activation and subsequent increase in cAMP leads to phosphorylation of CREB causing activation of the luciferase reporter gene that is under transcriptional control of the cAMP response element (CRE). The Isdop2 receptor showed a robust response to dopamine activation and greater signal-to background window relative to Isdop1 (FIG. 18), and accordingly Isdop2 was selected for assay development, optimization, and subsequent chemical library screen for receptor antagonists. Assay validation analysis was conducted using the maximal and minimal screening conditions in a checkerboard pattern to assess the robustness and utility of the assay platform and its parameters (e.g. cell number, dopamine and SCH23390 concentrations, and instrumentation). A concentration of 3 μM dopamine (the approximate EC90 value for Isdop2, labeled with an arrow in
To the best of our knowledge, the Isdop2 antagonist screen described herein is the first chemical library screen published for any tick GPCR. The average inhibitory effect of 10 mM SCH23390 (assay antagonist control) across all sixteen library compound plates was calculated as 86±6%. This value minus three standard deviations (i.e. 68%) was used to establish the cutoff for identification of hit compounds. Thus, compounds that resulted in an average inhibition of >68% of the dopamine-stimulated signal were reported as hits. Using these criteria, 85 hit compounds were identified in the Isdop2 antagonist screen, corresponding to a 7% hit rate (Table 20). The demonstrated pharmacological activities of the compounds in the LOPAC1280 small molecule collection, coupled with the composition of this library (>50% of library chemistries target GPCRs, including dopamine receptors) made the LOPAC1280 a suitable choice for screen validation and Isdop2 antagonist discovery. Use of this library also allowed comparisons to a similar screening effort directed at discovery of novel AaDOP2 antagonists (Meyer et al., 2012). Of the 85 hit compounds discovered here, a large fraction includes GPCR ligands with reported activity in mammalian systems. Mammalian dopamine receptor ligands (primarily benzazepines, phenothiazines, and thioxanthenes) constituted the largest class of hit compounds (26 compounds; 31% of hits). Interestingly, this class also included several known dopamine receptor agonists, as well as chemistries selective at Gai coupled dopamine receptors, suggesting a divergent pharmacological profile of Isdop2 compared to mammalian Gas-coupled dopamine receptors. Other major classes of hits included biogenic amine reuptake inhibitors (9 compounds; 11% of hits) and serotonin receptor ligands (9 compounds; 11% of hits), both of which contained multiple tricyclic antidepressant (TCA) derivatives. TCAs are widely used in medicine, and the observation that several TCAs act on the tick dopamine receptor is not entirely surprising because a number of TCAs have reported affinities for mammalian dopamine receptors at modest nanomolar concentrations. Another large group was comprised of protein kinase modulators (14 compounds; 16% of hits). The abundance of structurally diverse kinase modulators among the hits likely reflects off-target effects, rather than specific effects on the Isdop2 receptor, because the luciferase reporter system relies on protein kinase-dependent phosphorylation of CREB. In support of this hypothesis, we previously showed that the kinase inhibitor piceatannol, which was identified as an AaDOP2 antagonist in the luciferase reporter-based screen, did not exhibit measurable activity at the AaDOP2 target in a cAMP accumulation confirmation assay (Meyer et al., 2012).
In order to study the relative potency of a subset of the hit antagonistic compounds, we conducted secondary assays that employed direct measurement of cAMP. Eight compounds belonging to three chemical classes were selected for confirmation assays. These included four compounds from the dopamine receptor ligand class (butaclamol, clozapine, cis-(Z)-flupenthixol, and SCH23390), two serotonin receptor ligands (mianserin and methiothepin) and two biogenic amine reuptake inhibitors (amitriptyline and doxepin). Each of these compounds was also shown to have antagonist properties at the AaDOP2 receptor from the mosquito A. aegypti (Meyer et al., 2012). The effects of antagonist hits on the dopamine-stimulated cAMP accumulation of both Isdop1 and Isdop2 were evaluated in order to compare the pharmacological profiles of the receptors. All eight compounds showed inhibition of the dopamine-stimulated response of the Isdop1 and Isdop2 receptors, but considerable differences in potency and relative selectivity were observed (Table 21). Five of the eight compounds were more potent at the Isdop2 receptor (i.e. exhibited a lower IC50 value) compared to Isdop1, with the exceptions being SCH23390, doxepin, and clozapine, for which the range of IC50 values were overlapping at both receptors. An inherent difference between these two receptors is the high constitutive activity of Isdop1 compared to Isdop2, and in addition to their antagonistic properties, all eight compounds tested here also displayed inverse agonist activity at Isdop1 (
Our data provided an opportunity to compare two chemical library screen data sets between orthologous dopamine receptors from a tick (Table 20) and mosquito (Meyer et al., 2012). We previously reported a chemical screen of the LOPAC1280 library in which 51 candidate antagonists of the A. aegypti AaDOP2 receptor were identified (Meyer et al., 2012). Comparison of the hit chemistries identified in these screens shows an overlap of 46 compounds, and reveals 39 hits (46%) that are unique to
Isdop2, i.e. compounds that did not meet the hit selection criteria for AaDOP2 (Table 20, indicated with superscript d). Of these 39 Isdop2-specific hits, several have intermediate activity at AaDOP2. However, ten compounds (chlorpromazine, clemastine, clomipramine, cortexolone, cyclobenzaprine, cyclosporin A, cyproheptadine, CGP-74514A, PD 98,059, and Z-L-Phe chloromethyl ketone) had no or very minor effects (<25% inhibition) on the dopamine-stimulated activity of AaDOP2 (Meyer et al., 2012). Conversely, five hit compounds (JL-18, ketotifen, maprotiline, mitoxantrone, and blapachone) from the AaDOP2 screen (Meyer et al., 2012), were not revealed as Isdop2 antagonists. Interestingly, detailed examination of the screen data revealed that all of these five chemistries, representing five different chemical classes, showed an intermediate activity at the Isdop2 receptor (50e65% inhibition of the dopaminestimulated response) suggesting that they are also modest antagonists at Isdop2.
We previously reported the effects of these eight compounds (i.e. the ones selected for confirmation assays) in secondary cAMP accumulation confirmation assays at the AaDOP2 mosquito receptor and the human D1 dopamine receptor. Although each compound showed antagonistic effects on AaDOP2 and hD1, we observed important differences in the relative potencies of these compounds (Meyer et al., 2012). Together with the present data, we also noted differences in the responses to the antagonists between the orthologous Isdop2 and AaDOP2 arthropod receptors. For example, the rank order of potency for the antagonists were methiothepin>butaclamol and SCH23390>cis-(Z)-flupenthixol>amitriptyline>doxepin>clozapine compared to amitriptyline and methiothepin>cis-(Z)-flupenthixol>doxepin and clozapine>mianserin>butaclamol>SCH23390, for Isdop2 (Table 21) and AaDOP2 (Meyer et al., 2012), respectively. Notably, amitriptyline and methiothepin were the most potent AaDOP2 antagonists, whereas amitriptyline was the second and third least potent at Isdop1 and Isdop2, respectively. Additionally, whereas SCH23390 was identified as a potent antagonist at both I. scapularis receptors (Table 21) and the hD1, it had very weak antagonistic properties at the AaDOP2 receptor (approximately 3.000-fold less than that for hD1, Meyer et al., 2012). Moreover, whereas five of the eight chemistries were selective for AaDOP2 over hD1, none of the eight compounds were selective for the I. scapularis dopamine receptors over the human dopamine D1 receptor. The Isdop2 and AaDOP2 receptors share molecular characteristics, as shown by the amino acid alignment in
Despite limited sequence similarity between the dopamine receptors of ticks, insects, and vertebrates, our approach uncovered chemistries that exhibit cross reactivity with orthologous receptors in other species, including humans. The issues of receptor selectivity and off-target effects must be rigorously evaluated. This is a primary rationale for incorporation of the orthologous human D1 receptor into our early-phase drug discovery activities (Meyer et al., 2012). Where possible, the inclusion of additional vertebrate and invertebrate receptors could also greatly enhance initial analyses of compound specificity. The LOPAC 1280 library was employed in this study as it is enriched for human GPCR active chemistries and thus well suited for screen validation. Screening the LOPAC1280 library also allowed us to identify three receptor antagonists (i.e. amitriptyline, doxepin, and clozapine) that are FDA approved and have been used in humans for decades, suggesting the potential for drug re-purposing. However, expanded screening of tens and hundreds of thousands of chemistries is required to identify novel chemical scaffolds with in vitro activity against tick dopamine targets. These efforts will explore more unique “chemical space”, thus increasing the opportunity for discovery of molecules with limited off-target effects.
Our development, optimization, and validation of an assay platform for chemical library screening of a tick dopamine receptor may also be adapted to screen other arthropod Gas- and Gaicoupled receptors to discover chemistries that modulate their function. Implementation of this platform and subsequent confirmation assays resulted in identification of several antagonists/inverse agonists of the two D1-like dopamine receptors in I. scapularis. The comparative molecular and pharmacological analyses between these tick receptors with that of an orthologous mosquito dopamine receptor (Meyer et al., 2012) revealed commonalities and distinct pharmacological properties, which may be further exploited in efforts to develop species-specific chemistries for vector control. This research motivates future studies aimed to investigate the in vivo activities of dopamine receptor antagonists as potential leads for acaricide discovery.
HEK293 cells stably expressing both a CRELuc reporter construct and either of the receptors were stimulated with potential agonists. Dose-response curves were plotted and the EC50 values were calculated. Compounds with EC50 values ≧10 μM are considered to lack intrinsic activity at AaDOP2.
†Percent inhibition of receptor response in the presence of test compound relative to the SCH23390 control; *, SCH23390 “antagonist control”; ‡, compound analyzed in cAMP confirmation assay; CDK, cyclin dependent kinase; DAR, dopamine receptor; H, histamine receptor; Lck, lymphocyte-specific protein tyrosine kinase; NOR, norepinephrine; mAChR, muscarinic acetylcholine receptor; Syk, spleen tyrosine kinase; σ, sigma receptor; 5-HT, 5-hydroxytryptamine (serotonin). N.D.=not determined.
H. sapiens D1 dopamine receptor (hD1).
11 ± 1.9 nM
Select chemistries and the assay control (SCH23390) were tested in dose-response experiments in the presence of 3 μM dopamine in AaDOP2- or hD1-expressing cells (
aData reflect RT-PCR results using total RNA samples extracted from adult female A. aegypti.
D. melanogaster
A. mellifera
I. scapularis
H. sapiens
aCalculations of percent identify were based on the aligned TM sequences shown in FIG. 8.
aValues refer to the number of amino acids composing these features
bPresumed to form a disulfide bond for protein stabilization
cPredicted as important for binding the amine moieties of catecholamines
dImplicated in G-protein coupling
ePredicted to form hydrogen bonds with catechol hydroxyl groups
fAromatic residues implicated in ligand interaction
gN.P. = not predicted
a Chemistries were tested in duplicate at a single dose-point of 400 uM.
bRepresents average of seven technical replicates.
17 ± 0.88
a<10% inhibition at 3 μM test compound.
aResults represent three independent experiments;
bControl (water only).
†Experiment ongoing
Blatella germanica in cockroach contact assay.
aCompounds listed in rank order of toxicity;
bResults represent two independent experiments;
cRepresents three independent experiments;
dNegative control (DMSO:EtOH 1:1 v:v).
aCompounds listed in rank order of toxicity;
bNegative control (DMSO:EtOH 1:1 v:v).
cFipronil (500 μM).
Blatella germanica in cockroach contact assay
aResults represent two independent experiments (n = 2)
Table 17 Summary of amino acid similarity between the conserved TM domain regions of Isdop1 and Isdop2 in Ixodes scapularis with the D1-like dopamine receptors in Drosophila melanogaster (D-Dop1; Gotzes et al., 1994; Dop R99B (DAMB); Feng et al., 1996; Han et al. 1996), Apis mellifera (DOP1, DOP2: Mustard et al., 2003), and Homo sapiens (D1: Sunahara et al., 1990; D5; Sunahara et al., 1991).
I. scapularis
D. melanogaster
A. mellifera
H. sapiens
aCalculation based on the number of conserved amino acids among the combined TM domain sequences aligned in FIG. 15.
Table 18 Summary of key features in Isdop1 and Isdop2 based on their deducted amino acid sequences. The number of amino acids composing the N- and C-termini and the intracellular and extracellular loops were determined relative to the predicted transmembrane (TM) domain sequences shown in
aValues refer to the number of amino acids composing these features.
Table 18 Activity of Isdop1 and Isdop2 following stimulation with biogenic amines. Cyclic AMP accumulation was measured in the absence (basal) or presence of the indicated biogenic amines (10 μM). All experiments were conducted in the presence of 300 nM propranolol to block activation of endogenous β adrenergic receptors in HEK cells. The stimulatory response is shown as fold-over-basal, and the data include the mean±standard error of six independent experiments assayed in duplicate, *p<0.05, **p<0.01, ***p<0.001 using one-way ANOVA with Dunnett's post-hoc test comparing each drug treatment condition to basal.
a Calculation based on the number of conserved amino acids among the combined TM domain sequences aligned in FIG. 15.
Table 19 Hit compounds with antagonist properties at Isdop2 identified in the LOPAC1280 chemical library screen.
aPercent inhibition of response to 3 mM dopamine.
bReference compound used to establish hit criteria.
cCompound selected for in vitro confirmation assays.
dUnique Isdop2 hits (i.e. unique to the current screen, not identified as hits in the AaDOP2 screen, Meyer et al., 2012).
Table 18 Estimated IC50 values (in μM) for select antagonists of the Isdop1 and Isdop2 receptors in confirmation assays using a direct measurement of cAMP. The antagonistic properties are represented by average and 95% confidence intervals of IC50 values based on at least three independent experiments.
Abbreviation: DA, Dopamine.
The compounds in Table 19 were identified as candidate agonists of the AaDOP2 target in vitro using the screening methods provided against the LOPAC1280 library. Follow up CRCs of two biogenic amine “hits” show that both epinephrine and norepinephrine behave as partial agonists at the AaDOP2 receptor (
Additional libraries of compounds that can be screened in order to try and identify compounds that have an affinity for gene products that are substantially similar to those encoded by SEQ. ID NOS. 18. These libraries may contain compounds such as: Linear tricyclic compounds, bearing, when possible, hetero atoms in the central ring, including N, O, and S. Compounds with exocyclic bonds, e.g. amitriptyline and doxepin, will be catalytically reduced and the reduced molecules assayed to determine whether that bond is necessary for high activity. One or two members of the series that prove to have high activity will be compared with their secondary and primary amine analogues. For example, nortriptyline is the secondary amine analog of amitriptyline. We might also study the N-methylpiperazine analogues of amitriptyline and doxepin to determine whether an N,N-dimethylamino or N-methylpiperazine is optimum. Still other N,N-dialkyl compounds can also be tested.
Parallel comparisons can also be made of ring-substituted compounds with their non-substituted congeners, e.g., the 8-chloro atom can be removed from clozapine and that molecule (des-chloroclozapine) can be screen to determine if it interacts with these gene products and/or has insecticidal activity.
While the novel technology has been illustrated and described in detail in the figures and foregoing description, the same is to be considered as illustrative and not restrictive in character, it being understood that only the preferred embodiments have been shown and described and that all changes and modifications that come within the spirit of the novel technology are desired to be protected. As well, while the novel technology was illustrated using specific examples, theoretical arguments, accounts, and illustrations, these illustrations and the accompanying discussion should by no means be interpreted as limiting the technology. All patents, patent applications, and references to texts, scientific treatises, publications, and the like referenced in this application are incorporated herein by reference in their entirety.
This application claims the benefit of U.S. Provisional Patent Application No. 61/555,645, filed on Nov. 4, 2011, and incorporated herein by reference in its entirety.
This invention was made with government support under grant contract number 09-CRIS-2-2037/IND011474 awarded by the United States Department of Agriculture/National Institute of Food and Agriculture; and grant contract no. prime award W81XWH-10-1-0085, sub-award no. PU201596 awarded by the Department of Defense/Telemedicine and Advanced Technology Research Center. The United States Government has certain rights in the invention.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/US2012/063585 | 11/5/2012 | WO | 00 | 5/2/2014 |
Number | Date | Country | |
---|---|---|---|
61555645 | Nov 2011 | US |