Merging cultures and disciplines to create a drug discovery ecosystem at Virginia commonwealth university: Medicinal chemistry, structural biology, molecular and behavioral pharmacology and computational chemistry.

The Department of Medicinal Chemistry, together with the Institute for Structural Biology, Drug Discovery and Development, at Virginia Commonwealth University (VCU) has evolved, organically with quite a bit of bootstrapping, into a unique drug discovery ecosystem in response to the environment and culture of the university and the wider research enterprise. Each faculty member that joined the department and/or institute added a layer of expertise, technology and most importantly, innovation, that fertilized numerous collaborations within the University and with outside partners. Despite moderate institutional support with respect to a typical drug discovery enterprise, the VCU drug discovery ecosystem has built and maintained an impressive array of facilities and instrumentation for drug synthesis, drug characterization, biomolecular structural analysis and biophysical analysis, and pharmacological studies. Altogether, this ecosystem has had major impacts on numerous therapeutic areas, such as neurology, psychiatry, drugs of abuse, cancer, sickle cell disease, coagulopathy, inflammation, aging disorders and others. Novel tools and strategies for drug discovery, design and development have been developed at VCU in the last five decades; e.g., fundamental rational structure-activity relationship (SAR)-based drug design, structure-based drug design, orthosteric and allosteric drug design, design of multi-functional agents towards polypharmacy outcomes, principles on designing glycosaminoglycans as drugs, and computational tools and algorithms for quantitative SAR (QSAR) and understanding the roles of water and the hydrophobic effect.

Do 2-(Benzoyl)piperidines Represent a Novel Class of hDAT Reuptake Inhibitors?

2-(Benzoyl)piperidines (analogues of 1a), structural hybrids of the clinically employed ADHD medication methylphenidate (2) and the abused synthetic cathinone pentedrone (3), have been previously reported to act as novel and selective reuptake inhibitors of the human dopamine transporter (hDAT). One of the more potent benzoylpiperidines, as is the case with methylphenidate analogues, is its 3,4-dichloroaryl counterpart. Here, we demonstrate using homology models that these compounds (i.e., benzoylpiperidines and methylphenidate analogues) likely bind in a comparable manner at hDAT. In addition, it is shown here that the 3,4-dichlorobenzoylpiperidine analogue of 1a is more potent than its 3,4-dimethyl counterpart, suggesting that the electronic character of the substituents might play a role in the potency of these hybrids. Furthermore, the 3,4-benz-fused (i.e., naphthyl) benzoylpiperidine analogue acts in the same manner as its corresponding methylphenidate counterpart at hDAT. As with its methylphenidate counterpart, the naphthyl compound also acts, rather uniquely (although with lower potency) relative to other members of the two series, at the human serotonin transporter (hSERT). In conclusion, the benzoylpiperidines represent a novel structural class of hDAT reuptake inhibitors that function in a manner similar to their methylphenidate counterparts.

Binding and functional structure-activity similarities of 4-substituted 2,5-dimethoxyphenyl isopropylamine analogues at 5-HT2A and 5-HT2B serotonin receptors.

Certain 4-substituted analogs of 1-(2,5-dimethoxyphenyl)isopropylamine (2,5-DMA) are psychoactive classical hallucinogens or serotonergic psychedelic agents that function as human 5-HT2A (h5-HT2A) serotonin receptor agonists. Activation of a related receptor population, h5-HT2B receptors, has been demonstrated to result in adverse effects including cardiac valvulopathy. We previously published on the binding of several such agents at the two receptor subtypes. We hypothesized that, due to their structural similarity, the 5-HT2A and 5-HT2B receptor affinities of these agents might be related, and that QSAR studies might aid future studies. For a series of 13 compounds, it is demonstrated here that i) their published rat brain 5-HT2 receptor affinities are significantly correlated with their h5-HT2A (r = 0.942) and h5-HT2B (r = 0.916) affinities, ii) as with r5-HT2 receptor affinity, h5-HT2A affinity is correlated with the lipophilicity of the 4-position substituent (r = 0.798), iii) that eight of the ten compounds examined in functional (Ca+2 mobilization in stable cell lines generated expressing the human 5-HT2B receptor using the Flp-In T-REx system) assays acted as h5-HT2B agonists (4-substituent = H, F, Br, I, OCH2CH3, NO2, nC3H7, tC4H9) and two (n-hexyl and benzyl) as antagonists, iv) h5-HT2B affinity but not action was correlated with the lipophilicity of the 4-position substituent (r = 0.750; n = 10). The findings suggest that h5-HT2B receptor affinity, and its relationship to substituent lipophilicity, might be approximated by rat and h5-HT2A affinity but cannot be used as a predictor of h5-HT2B agonist action of 2,5-DMA analogs. Furthermore, given that certain 2,5-DMA analogs are on the clandestine market, their potential to produce cardiac side effects following persistent or chronic use via activation of h5-HT2B receptors should be considered.

Medicinal chemistry: The key to critical thinking in pharmacotherapy.

INTRODUCTION: Medicinal chemistry is a polarizing subject for pharmacy students where, if not embraced, future pharmacists may be limited in their role as drug experts. An understanding of medicinal chemistry and its structure-activity relationships creates a strong foundation upon which our knowledge of pharmacotherapy is built. PERSPECTIVE: As the field of pharmacy has shifted to an increasingly clinical role, with an emphasis on patient care as a member of the interprofessional team, pharmacy has also seen an increase in postgraduate training, specifically residencies and fellowships. Pharmacy students noting this trend may depreciate medicinal chemistry early in the curriculum and place more focus on therapeutics and clinical rotations. However, forgoing the fundamental understanding of medicinal chemistry may hinder pharmacy students' current breadth and understanding, and the ability to rationalize future developments in their practice. Medicinal chemistry empowers pharmacists with the ability to reason through medications' impact versus simply memorizing their actions. Pharmacists play a unique role as drug experts, with advanced problem-solving and critical thinking skills that set them apart from drug references and search engines. IMPLICATIONS: As the field moves towards pharmacists as a member of the clinical team, the faculty should integrate medicinal chemistry throughout the doctor of pharmacy curricula. Faculty without this ability for a curriculum change should consider integration in their content. The field of pharmacy must take care to not allow clinical knowledge to significantly overshadow the importance of medicinal chemistry or run the risk of saturating the field with underprepared pharmacists.

Non-conserved residues dictate dopamine transporter selectivity for the potent synthetic cathinone and psychostimulant MDPV.

Clandestine chemists are currently exploiting the pyrrolidinophenone scaffold to develop new designer drugs that carry the risk of abuse and overdose. These drugs promote addiction through the rewarding effects of increased dopaminergic neurotransmission. 3,4-Methylenedioxypyrovalerone (MDPV) and its analogs are illicit psychostimulants of this class that are ∼50-fold more potent than cocaine at inhibiting the human dopamine transporter (hDAT). In contrast, MDPV is a weak inhibitor at both the human serotonin transporter (hSERT) and, as it is shown here, the Drosophila melanogaster DAT (dDAT). We studied three conserved residues between hSERT and dDAT that are unique in hDAT (A117, F318, and P323 in dDAT), and one residue that is different in all three transporters (D121 in dDAT). hDAT residues were replaced in the dDAT sequence at these positions using site-directed mutagenesis and stable cell lines were generated expressing these mutant transporters. The potencies of MDPV and two of its analogs were determined using a Ca2+-mobilization assay. In this assay, voltage-gated Ca2+ channels are expressed to sense the membrane electrical depolarization evoked when dopamine is transported through DAT. Each individual mutant slightly improved MDPV's potency, but the combination of all four increased its potency ∼100-fold (2 log units) in inhibiting dDAT activity. Molecular modeling and docking studies were conducted to explore the possible mode of interaction between MDPV and DAT in silico. Two of the studied residues (F318 and P323) are at the entrance of the S1 binding site, whereas the other two (A117 and D121) face the aryl moiety of MDPV when bound to this site. Therefore, these four non-conserved residues can influence MDPV selectivity not only by stabilizing binding, but also by controlling access to its binding site at DAT.

Computational analysis of non-competitive antagonist arylguanidine-α7 nAChR complexes.

meta-Chlorophenylguanidine (1) is a non-competitive α7 nicotinic acetylcholine receptor (nAChR) antagonist. Here we examined the hydrogen bond donor role of the anilinic N1-H on the inhibitory effect of 1 by preparing its N1-CH3 counterpart 2. Analog 2 was found to be at least as potent as 1 as a non-competitive α7 nAChR antagonist in a patch-clamp assay. To establish a structural basis for the mode of interaction of guanidines 1 and 2, we generated 100 homology models of the hα7 nAChR. This was followed by Connolly surface (SYBYL-X2.1) and blind docking (AutoDock 4.1) studies to identify eight possible binding pockets, two of which were supported by empirical data and employed in our docking studies. The optimized model-ligand complexes were analyzed using a Hydropathic INTeractions (HINT) analysis in order to compare and contrast different binding pockets and modes. We identified a potential allosteric binding site and distinct rotameric binding modes for 1 and 2 at α7 nAChRs. These differences in the binding orientations minimized the importance of an anilinic NH function for the antagonist activity at nACh receptors.

N1H- and N1-Substituted Phenylguanidines as α7 Nicotinic Acetylcholine (nACh) Receptor Antagonists: Structure-Activity Relationship Studies.

We previously reported that N-(3-chlorophenyl)guanidine (1) represents a novel α7 nicotinic ACh (nACh) receptor antagonist chemotype. In the present study, a small series of compounds was synthesized with the intent to investigate the structure-activity relationship (SAR). Preliminary data suggested that the N-methyl analog of 1, 2, was several times more potent. Therefore, the chloro group at the aryl 3-position of 1 and its N1-methyl counterpart 2 were replaced with a number of substituents considering the electronic, lipophilic, and steric nature of the substituents. The potencies of the compounds to inhibit acetylcholine (ACh)-induced responses were obtained in Xenopus laevis oocytes expressing human α7 nicotinic ACh receptors (nAChRs) using a two-electrode voltage-clamp assay. We found that the nature of the 3-position substituents had relatively little (i.e., <10-fold) effect on potency, and the presence of an N1-isopropyl substituent was tolerated. Here, we report the first SAR investigation of this novel α7 nAChR antagonist chemotype.

Psychedelic-like Properties of Quipazine and Its Structural Analogues in Mice.

Known classic psychedelic serotonin 2A receptor (5-HT2AR) agonists retain a tryptamine or phenethylamine at their structural core. However, activation of the 5-HT2AR can be elicited by drugs lacking these fundamental scaffolds. Such is the case of the N-substituted piperazine quipazine. Here, we show that quipazine bound to and activated 5-HT2AR as measured by [3H]ketanserin binding displacement, Ca2+ mobilization, and accumulation of the canonical Gq/11 signaling pathway mediator inositol monophosphate (IP1) in vitro and in vivo. Additionally, quipazine induced via 5-HT2AR an expression pattern of immediate early genes (IEG) in the mouse somatosensory cortex consistent with that of classic psychedelics. In the mouse head-twitch response (HTR) model of psychedelic-like action, quipazine produced a lasting effect with high maximal responses during the peak effect that were successfully blocked by the 5-HT2AR antagonist M100907 and absent in 5-HT2AR knockout (KO) mice. The acute effect of quipazine on HTR appeared to be unaffected by serotonin depletion and was independent from 5-HT3R activation. Interestingly, some of these features were shared by its deaza bioisostere 2-NP, but not by other closely related piperazine congeners, suggesting that quipazine might represent a distinct cluster within the family of psychoactive piperazines. Together, our results add to the mounting evidence that quipazine's profile matches that of classic psychedelic 5-HT2AR agonists at cellular signaling and behavioral pharmacology levels.

X-Ray Crystal Structure of a 2-Amino-3,4-dihydroquinazoline 5-HT3 Serotonin Receptor Antagonist and Related Analogs.

Certain 2-amino-3,4-dihydroquinazolines bind at 5-HT3 serotonin receptors and act as antagonists (e.g. 6-chloro) whereas others bind with little to no affinity and lack functional activity (e.g. 8-chloro). The purpose of this investigation was to gain insight as to why this might be the case. X-Ray crystallographic studies revealed that the N-C-N distances in the examined analogs are nearly identical (1.31 - 1.34 Å), suggesting that differences in N-C-N delocalization does not account for differences in affinity/action. Homology modeling hydrophatic interactions (HINT) analysis revealed that the 6-chloro analog formed a greater number, and more favorable, interactions with the receptor, whereas the 8-chloro analog formed fewer, and unfavorable, interactions. The affinity and activity of the 6-chloro quinazoline relative to its 8-chloro counterpart are unrelated to the N-C-N delocalization pattern but might be related to specific (favorable and unfavorable) interactions of quinazoline substituents with certain receptor features as determined by HINT analysis.

Multi-modal antidepressant-like action of 6- and 7-chloro-2-aminodihydroquinazolines in the mouse tail suspension test.

RATIONALE: 2-Amino-6-chloro-3,4-dihydroquinazoline (e.g., A6CDQ) represents a novel putative antidepressant originally thought to act through a 5-HT3 serotonin receptor antagonist mechanism. Here, we investigated this further by examining a positional isomer of A6CDQ (i.e., A7CDQ). MATERIALS AND METHODS: 5-HT3 receptor and transporter activity (uptake-1 and uptake-2) were investigated using a variety of in vitro assays and the in vivo mouse tail suspension test (TST). RESULTS: Although A7CDQ binds at 5-HT3 receptors with low affinity (Ki = 1975 nM) compared to A6CDQ (Ki = 80 nM), it retained 5-HT3 receptor antagonist action (IC50 = 5.77 and 0.26 µM, respectively). In the mouse TST A7CDQ produced antidepressant-like actions (ED50 = 0.09 mg/kg) comparable to that of A6CDQ. In addition, A6CDQ was found to be a 5-HT releasing agent (Km = 2.8 µM) at hSERT and a reuptake inhibitor (IC50 = 1.8 µM) at hNET, whereas A7CDQ was a weak reuptake inhibitor (Km = 43.6 µM) at SERT but a releasing agent (EC50 = 3.3 µM) at hNET. Moreover, A6CDQ and A7CDQ were potent inhibitors of uptake-2 (e.g.; OCT3 IC50 = 3.9 and 5.9 µM, respectively). CONCLUSIONS: A simple shift of a substituent in a common quinazoline scaffold from one position to another (i.e., a chloro group from the 6- to the 7-position) resulted in a common action in the TST but via a somewhat different mechanism. A6CDQ and A7CDQ might represent the first members of a new class of potential antidepressants with a unique multi-modal mechanism of action.

Revised Pharmacophore Model for 5-HT2A Receptor Antagonists Derived from the Atypical Antipsychotic Agent Risperidone.

Pharmacophore models for 5-HT2A receptor antagonists consist of two aromatic/hydrophobic regions at a given distance from a basic amine. We have previously shown that both aromatic/hydrophobic moieties are unnecessary for binding or antagonist action. Here, we deconstructed the 5-HT2A receptor antagonist/serotonin-dopamine antipsychotic agent risperidone into smaller structural segments that were tested for 5-HT2A receptor affinity and function. We show, again, that the entire risperidone structure is unnecessary for retention of affinity or antagonist action. Replacement of the 6-fluoro-3-(4-piperidinyl)-1,2-benz[ d]isoxazole moiety by isosteric tryptamines resulted in retention of affinity and antagonist action. Additionally, 3-(4-piperidinyl)-1,2-benz[ d]isoxazole (10), which represents less than half the structural features of risperidone, retains both affinity and antagonist actions. 5-HT2A receptor homology modeling/docking studies suggest that 10 binds in a manner similar to risperidone and that there is a large cavity to accept various N4-substituted analogues of 10 such as risperidone and related agents. Alterations of this "extended" moiety improve receptor binding and functional potency. We propose a new risperidone-based pharmacophore for 5-HT2A receptor antagonist action.

"Methylene Bridge" to 5-HT3 Receptor Antagonists: Conformationally Constrained Phenylguanidines.

Arylguanidines, depending upon their aromatic substitution pattern, display varying actions at 5-HT3 receptors (e.g., partial agonist, agonist, superagonist). Here, we demonstrate that conformational constraint of these agents as dihydroquinazolines (such as A6CDQ; 1) results in their conversion to 5-HT3 receptor antagonists. We examined the structure-activity relationships of 1. Replacement/removal of any of the guanidinium nitrogen atoms of 1 resulted in decreased affinity. All three nitrogen atoms of 1 are necessary for optimal binding affinity at 5-HT3 receptors. Introduction of substituents as small as an N2-methyl group abolishes affinity. The results are consistent with homology modeling/docking studies and binding data from site-directed mutagenesis studies. Introducing a "methylene bridge" to the arylguanidine structure, regardless of its functional activity, results in a 5-HT3 receptor antagonist.

des-Formylflustrabromine (dFBr): A Structure-Activity Study on Its Ability To Potentiate the Action of Acetylcholine at α4ß2 Nicotinic Acetylcholine Receptors.

The naturally occurring indole alkaloid des-formylflustrabromine (dFBr; 1) is one of the first agents shown to act as a selective positive allosteric modulator (PAM) at α4ß2 nicotinic acetylcholine receptors (nAChRs). We previously deconstructed this agent to determine which of its structural features contribute to its actions and have identified an agent that might serve as the basis for a " working pharmacophore". Here, we elaborate the dFBr (1; EC50 = 0.2 µM) structure to identify how various structural modifications impact its actions. Electrophysiological studies with Xenopus laevis oocytes identified several compounds with dFBr-like potency and one, the 5-bromo analogue of 1 (i.e., 5-bromo dFBr; 25; EC50 = 0.4 µM), with more than twice the efficacy of 1 as a PAM at α4ß2 nAChRs.

A new chemotype inhibitor for the human organic cation transporter 3 (hOCT3).

Human organic cation transporters (OCTs) represent an understudied neurotransmitter uptake mechanism for which no selective agents have yet been identified. Several neurotransmitters (e.g. serotonin, norepinephrine) are low-affinity substrates for these transporters, but possess higher affinity for other transporters (e.g. the serotonin or norepinephrine transporters; SERT and NET, respectively). We have identified a new class of OCT inhibitors with a phenylguanidine structural scaffold. Here, we examine the actions of a series of such compounds and report preliminary structure-activity relationships (SARs) - the first dedicated SAR study of OCT3 action. Initial results showed that the presence of a substituent on the phenyl ring, as well as its position, contributes to the phenylguanidines' inhibitory potency (IC50 values ranging from 2.2 to >450µM) at hOCT3. There is a trend towards enhanced inhibitory potency of phenylguanidines with increased lipophilic character and the size of the substituent at the phenyl 4-position, with the latter reaching a ceiling effect. The first PiPT-based hOCT3 homology models were generated and are in agreement with our biological data.

Reevaluation of fenpropimorph as a σ receptor ligand: Structure-affinity relationship studies at human σ1 receptors.

Fenpropimorph (1) is considered a "super high-affinity" σ1 receptor ligand (Ki=0.005nM for guinea pig σ1 receptors). Here, we examine the binding of 1 and several of its deconstructed analogs at human σ1 (hσ1) receptors. We monitored their subtype selectivity by determining the binding affinity at σ2 receptors. In addition, we validated an existing pharmacophore model at the molecular level by conducting 3D molecular modeling studies, using the crystal structure of hσ1 receptors, and Hydrophatic INTeractions (HINT) analysis. Our structure affinity relationship studies showed that 1 binds with lower affinity at hσ1 receptors (Ki=17.3nM) compared to guinea pig; moreover, we found that none of the fenpropimorph methyl groups is important for its binding at hσ1 receptors, nor is stereochemistry. For example, removal of all methyl groups as seen in 4 resulted in an almost 5-fold higher affinity at hσ1 receptors compared to 1 and 350-fold selectivity versus σ2 receptors. In addition, although the O atom of the morpholine ring does not contribute to affinity at hσ1 receptors (and might even detract from it), it plays role in subtype (σ1 versus σ2 receptor) selectivity.

Synthetic Cathinones: A Brief Overview of Overviews with Applications to the Forensic Sciences.

Catha edulis, the fresh leaves of which (i.e., khat) are used for their central stimulant actions, has been known for many hundreds of years. S(-)Cathinone was identified as a centrally-active khat constituents >30 years ago. Although khat use was a problem long localized to certain Middle Eastern and certain Eastern African nations, 'synthetic cathinones' (synthetic analogs of cathinone) represent a "new" class of abused substances with growing worldwide appeal. To date, about 150 synthetic cathinones have been identified on the clandestine market, and only a dozen or so have been controlled (U.S. Schedule I). Because these agents do not represent a pharmacologically (i.e., behavioral) or mechanistically homogeneous class of agents, synthetic cathinones are being investigated one agent at a time to understand their actions, mechanisms of action, metabolism, toxicity, and abuse potential - the latter to identify possible modes of overdose treatment and for Scheduling purposes. The available agents might represent only the 'tip of the iceberg'; the potential for many more new synthetic cathinones is very real. Investigation of these agents on an agent-by-agent basis is a daunting task. Attempts are being made to understand these agents as a class, by examination of their structure-activity relationships. Here, we provide an overview of review articles that attempts to shed light on these agents as a class, rather than on an agent-by-agent basis. This article is meant to be a reference resource that might expedite the work of those in this field by directing them to where they can find useful information.

Structure-Activity Relationships of Synthetic Cathinones.

Until recently, there was rather little interest in the structure-activity relationships (SARs) of cathinone analogs because so few agents were available and because they represented a relatively minor drug abuse problem. Most of the early SAR was formulated on the basis of behavioral (e.g., locomotor and drug discrimination) studies using rodents. With the emergence on the clandestine market in the last few years of a large number of new cathinone analogs, termed "synthetic cathinones", and the realization that they likely act at dopamine, norepinephrine, and/or serotonin transporters as releasing agents (i.e., as substrates) or reuptake inhibitors (i.e., as transport blockers), it has now become possible to better examine their SAR and even their quantitative SAR (QSAR), in a more effective and systematic manner. An SAR picture is beginning to emerge, and key structural features, such as the nature of the terminal amine, the size of the α-substituent, stereochemistry, and the presence and position of aromatic substituents, are being found to impact action (i.e., as releasing agents or reuptake inhibitors) and transporter selectivity.

Superagonist, Full Agonist, Partial Agonist, and Antagonist Actions of Arylguanidines at 5-Hydroxytryptamine-3 (5-HT3) Subunit A Receptors.

Introduction of minor variations to the substitution pattern of arylguanidine 5-hydroxytryptamine-3 (5-HT3) receptor ligands resulted in a broad spectrum of functionally-active ligands from antagonist to superagonist. For example, (i) introduction of an additional Cl-substituent(s) to our lead full agonist N-(3-chlorophenyl)guanidine (mCPG, 2; efficacy % = 106) yielded superagonists 7-9 (efficacy % = 186, 139, and 129, respectively), (ii) a positional isomer of 2, p-Cl analog 11, displayed partial agonist actions (efficacy % = 12), and (iii) replacing the halogen atom at the meta or para position with an electron donating OCH3 group or a stronger electron withdrawing (i.e., CF3) group resulted in antagonists 13-16. We posit based on combined mutagenesis, crystallographic, and computational analyses that for the 5-HT3 receptor, the arylguanidines that are better able to simultaneously engage the primary and complementary subunits, thus keeping them in close proximity, have greater agonist character while those that are deficient in this ability are antagonists.

Reformulating a Pharmacophore for 5-HT2A Serotonin Receptor Antagonists.

Several pharmacophore models have been proposed for 5-HT2A serotonin receptor antagonists. These typically consist of two aromatic/hydrophobic moieties separated by a given distance from each other, and from a basic amine. Although specified distances might vary, the models are relatively similar in their general construction. Because our preliminary data indicated that two aromatic (hydrophobic) moieties might not be required for such action, we deconstructed the serotonin-dopamine antipsychotic agent risperidone (1) into four smaller structural fragments that were thoroughly examined in 5-HT2A receptor binding and functional (i.e., two-electrode voltage clamp (TEVC) and intracellular calcium release) assays. It was apparent that truncated risperidone analogues behaved as antagonists. In particular, 6-fluoro-3-(1-methylpiperidin-4-yl)benzisoxazole (4) displayed high affinity for 5-HT2A receptors (Ki of ca. 12 nM) relative to risperidone (Ki of ca. 5 nM) and behaved as a potent 5-HT2A serotonin receptor antagonist. These results suggest that multiple aromatic (hydrophobic) moieties are not essential for high-affinity 5-HT2A receptor binding and antagonist activity and that current pharmacophore models for such agents are very much in need of revision.

Receptor-ligand interaction at 5-HT3 serotonin receptors: a cluster approach.

A fundamental understanding of the interaction of ligands with biological receptors is important because many drugs exert their influence via receptors. Using a cluster approach, we have studied the role of structural and electronic parameters on receptor-ligand binding by carrying out density functional theory based calculations. As model systems, we have studied substituted arylguanidines, which activate 5-HT3 receptors in a manner similar to that of serotonin. The geometries of the arylguanidine derivatives were fully optimized to obtain the lowest energy structures. Electronic properties such as binding energies, dipole moments, polarizabilities, and electron affinities, as well as geometric properties, such as molecular volume and dihedral angles were calculated, and their relationship with binding affinity was evaluated. Results obtained were compared to experimental ligand-receptor binding affinity data available. These fundamental studies show that though both electronic and geometric properties of the ligands are important for binding, the electron affinities of the substituent species play a dominant role. Potential new fundamental indices for ligand-receptor affinity are also discussed.