Structure Of Ligand Binding Research Articles

Comparison of Several Molecular Docking Programs: Pose Prediction and Virtual Screening Accuracy

Molecular docking programs are widely used modeling tools for predicting ligand binding modes and structure based virtual screening. In this study, six molecular docking programs (DOCK, FlexX, GLIDE, ICM, PhDOCK, and Surflex) were evaluated using metrics intended to assess docking pose and virtual screening accuracy. Cognate ligand docking to 68 diverse, high-resolution X-ray complexes revealed that ICM, GLIDE, and Surflex generated ligand poses close to the X-ray conformation more often than the other docking programs. GLIDE and Surflex also outperformed the other docking programs when used for virtual screening, based on mean ROC AUC and ROC enrichment values obtained for the 40 protein targets in the Directory of Useful Decoys (DUD). Further analysis uncovered general trends in accuracy that are specific for particular protein families. Modifying basic parameters in the software was shown to have a significant effect on docking and virtual screening results, suggesting that expert knowledge is critical for optimizing the accuracy of these methods.

Protein Polarization Is Critical to Stabilizing AF-2 and Helix-2′ Domains in Ligand Binding to PPAR-γ

The peroxisome proliferator-activated receptor (PPAR-gamma) is a ligand-dependent transcription factor that is important in adipocyte differentiation and glucose homeostasis. This paper presents a detailed dynamics study of PPAR-gamma and its binding to the agonist rosiglitazone using both polarized and unpolarized force fields. The numerical result revealed the critical role of protein polarization in stabilizing the activation function-2 (AF-2) in ligand binding to PPAR-gamma and a helix structure (helix-2'). Specifically when nonpolarized force field is used, a critical H-bond in PPAR-gamma binding is broken, which caused AF-2 to adopt random structures. In addition, helix-2' is partially denatured during the MD simulation, due to the breaking of a backbone hydrogen bond. In contrast, when polarized force field is employed in MD simulation, the PPAR-gamma ligand binding structure is stabilized and the local structure of helix-2' remains folded, both being in excellent agreement with experimental observations. The current result demonstrates the importance of electronic polarization of protein in stabilizing hydrogen bonding, which is critical to preserving the native structure of local helices and protein-ligand binding in PPAR-gamma.

The Interplay of Ligand Binding and Quaternary Structure in the Diverse Interactions of Dynein Light Chain LC8

Dynein light chain LC8 is a small, dimeric, and very highly conserved globular protein that is an integral part of the dynein and myosin molecular motors but appears to have a broader role in multiple protein complexes unrelated to molecular motors. LC8 binds to two families of targets: those having a KXTQT sequence fingerprint and those having a GIQVD fingerprint. All known LC8 binding partners containing these fingerprints share a common binding site on LC8 that raises the question of what determines binding specificity. Here, we present the crystal structure of apo-LC8 at 1.7-Å resolution, which, when compared with the crystal structures of several LC8 complexes, gives insight into the mechanism underlying the binding diversity of LC8. Peptide binding is associated with a shift in quaternary structure that expands the hydrophobic binding surface available to the ligand, in addition to changes in tertiary structure and ordering of LC8 around the binding groove. The observed quaternary shift suggests a mechanism by which binding at one of the two identical sites can influence binding at the other. NMR spectra of titrations with peptides from each fingerprint family show evidence of allosteric interaction between the two binding sites, to a differing degree in the two ligand families. Allosteric interaction between the binding sites may be a mechanism to promote simultaneous binding of ligands from the same family, providing a physiological role for the two fingerprints.

A Functional Peroxisome Proliferator-activated Receptor-γ Ligand-binding Domain Is Not Required for Adipogenesis

The nuclear hormone receptor peroxisome proliferator-activated receptor-gamma (PPARgamma) is the central regulator of adipogenesis. Although it is the target for several drugs that function as agonist activators, a high affinity endogenous ligand for this receptor that is involved in regulating adipogenesis has yet to be identified. Here, we investigated the requirement for ligand activation of PPARgamma in fat cell differentiation, taking advantage of a natural mutant of this receptor that does not bind or become activated by any known natural or synthetic ligand. When ectopically expressed in PPARgamma-null fibroblasts, this Q286P allele was able to strongly promote morphological adipogenesis, without any significant difference compared with wild-type PPARgamma. In addition, no significant differences were found in the expression of several adipogenic genes between the wild-type and Q286P mutant alleles. To extend our studies to an in vivo setting, we performed subcutaneous injections of PPARgamma-expressing fibroblasts into nude mice. We found that both wild-type and Q286P mutant-expressing fibroblasts were able to generate fat pads in the mice. These results suggest that the binding and activation of PPARgamma by agonist ligands may not be required for adipogenesis under physiological conditions.

Average and extreme multi-atom Van der Waals interactions: Strong coupling of multi-atom Van der Waals interactions with covalent bonding

BackgroundThe prediction of ligand binding or protein structure requires very accurate force field potentials – even small errors in force field potentials can make a 'wrong' structure (from the billions possible) more stable than the single, 'correct' one. However, despite huge efforts to optimize them, currently-used all-atom force fields are still not able, in a vast majority of cases, even to keep a protein molecule in its native conformation in the course of molecular dynamics simulations or to bring an approximate, homology-based model of protein structure closer to its native conformation.ResultsA strict analysis shows that a specific coupling of multi-atom Van der Waals interactions with covalent bonding can, in extreme cases, increase (or decrease) the interaction energy by about 20–40% at certain angles between the direction of interaction and the covalent bond. It is also shown that on average multi-body effects decrease the total Van der Waals energy in proportion to the square root of the electronic component of dielectric permittivity corresponding to dipole-dipole interactions at small distances, where Van der Waals interactions take place.ConclusionThe study shows that currently-ignored multi-atom Van der Waals interactions can, in certain instances, lead to significant energy effects, comparable to those caused by the replacement of atoms (for instance, C by N) in conventional pairwise Van der Waals interactions.

Raloxifene and ICI182,780 Increase Estrogen Receptor-α Association with a Nuclear Compartment via Overlapping Sets of Hydrophobic Amino Acids in Activation Function 2 Helix 12

The basis for the differential repressive effects of antiestrogens on transactivation by estrogen receptor-alpha (ERalpha) remains incompletely understood. Here, we show that the full antiestrogen ICI182,780 and, to a lesser extent, the selective ER modulator raloxifene (Ral), induce accumulation of exogenous ERalpha in a poorly soluble fraction in transiently transfected HepG2 or stably transfected MDA-MB231 cells and of endogenous receptor in MCF7 cells. ERalpha remained nuclear in HepG2 cells treated with either compound. Replacement of selected hydrophobic residues of ERalpha ligand-binding domain helix 12 (H12) enhanced receptor solubility in the presence of ICI182,780 or Ral. These mutations also increased transcriptional activity with Ral or ICI182,780 on reporter genes or on the endogenous estrogen target gene TFF1 in a manner requiring the integrity of the N-terminal AF-1 domain. The antiestrogen-specific effects of single mutations suggest that they affect receptor function by mechanisms other than a simple decrease in hydrophobicity of H12, possibly due to relief from local steric hindrance between these residues and the antiestrogen side chains. Fluorescence anisotropy experiments indicated an enhanced regional stabilization of mutant ligand-binding domains in the presence of antiestrogens. H12 mutations also prevent the increase in bioluminescence resonance energy transfer between ERalpha monomers induced by Ral or ICI182,780 and increase intranuclear receptor mobility in correlation with transcriptional activity in the presence of these antiestrogens. Our data indicate that ICI182,780 and Ral locally alter the ERalpha ligand binding structure via specific hydrophobic residues of H12 and decrease its transcriptional activity through tighter association with an insoluble nuclear structure.

Ligand–DNA interaction in a nanocage of reverse micelle

We have studied intercalation of ethidium bromide (EB) to genomic DNA encapsulated in a nanospace of an anionic AOT reverse micelle (RM). Circular dichroism (CD) study on the DNA in the RM reveals its condensed form. Here, we have used temporal decay-associated spectra (DAS) and time-resolved area normalized emission spectral (TRANES) techniques to investigate EB-binding to condensed DNA because the interference of emission from unbound EB in the RM makes conventional steady state and picosecond resolved fluorescence spectroscopic techniques challenging. The binding affinity of the ligand EB with the DNA in the RM is found to increase with the size of the RM, reflecting the effect of lessening of DNA condensation on the binding affinity. CD spectra of the DNA in the RM with various sizes indicate the structural change of the condensed DNA with reverse micellar size. DAS and TRANES techniques along with dynamic light scattering studies of the EB-DNA complex in the RM further reveal two kinds of binding modes of the ligand with the condensed DNA even in essentially monodispersed RMs. To investigate the role of RM on the ligand binding and secondary structure of the DNA, we have also studied complexation of EB with two synthetic self-complimentary oligonucleotides of sequences (CGCAAATTTGCG)2 and (CGCGCGCGCGCG)2 in the RM.

Locating an Antagonist in the 5-HT3 Receptor Binding Site Using Modeling and Radioligand Binding

We have used a homology model of the extracellular domain of the 5-HT(3) receptor to dock granisetron, a 5-HT(3) receptor antagonist, into the binding site using AUTODOCK. This yielded 13 alternative energetically favorable models. The models fell into 3 groups. In model type A the aromatic rings of granisetron were between Trp-90 and Phe-226 and its azabicyclic ring was between Trp-183 and Tyr-234, in model type B this orientation was reversed, and in model type C the aromatic rings were between Asp-229 and Ser-200 and the azabicyclic ring was between Phe-226 and Asn-128. Residues located no more than 5 A from the docked granisetron were identified for each model; of 26 residues identified, 8 were found to be common to all models, with 18 others being represented in only a subset of the models. To identify which of the docking models best represents the ligand-receptor complex, we substituted each of these 26 residues with alanine and a residue with similar chemical properties. The mutant receptors were expressed in human embryonic kidney (HEK)293 cells and the affinity of granisetron determined using radioligand binding. Mutation of 2 residues (Trp-183 and Glu-129) ablated binding, whereas mutation of 14 other residues caused changes in the [(3)H]granisetron binding affinity in one or both mutant receptors. The data showed that residues both in and close to the binding pocket can affect antagonist binding and overall were found to best support model B.

PDB-Ligand: a ligand database based on PDB for the automated and customized classification of ligand-binding structures

PDB-Ligand (http://www.idrtech.com/PDB-Ligand/) is a three-dimensional structure database of small molecular ligands that are bound to larger biomolecules deposited in the Protein Data Bank (PDB). It is also a database tool that allows one to browse, classify, superimpose and visualize these structures. As of May 2004, there are about 4870 types of small molecular ligands, experimentally determined as a complex with protein or DNA in the PDB. The proteins that a given ligand binds are often homologous and present the same binding structure to the ligand. However, there are also many instances wherein a given ligand binds to two or more unrelated proteins, or to the same or homologous protein in different binding environments. PDB-Ligand serves as an interactive structural analysis and clustering tool for all the ligand-binding structures in the PDB. PDB-Ligand also provides an easier way to obtain a number of different structure alignments of many related ligand-binding structures based on a simple and flexible ligand clustering method. PDB-Ligand will be a good resource for both a better interpretation of ligand-binding structures and the development of better scoring functions to be used in many drug discovery applications.

Structural plasticity of peanut lectin: an X-ray analysis involving variation in pH, ligand binding and crystal structure.

Until recently, it has only been possible to grow crystals of peanut lectin when complexed with sugar ligands. It is now shown that it is possible to grow peanut lectin crystals at acidic pH in the presence of oligopeptides corresponding to a loop in the lectin molecule. Crystals have also been prepared in the presence of these peptides as well as lactose. Low-pH crystal forms of the lectin-lactose complex similar to those obtained at neutral pH have also been grown. Thus, crystals of peanut lectin grown under different environmental conditions, at two pH values with and without sugar bound to the lectin, are now available. They have been used to explore the plasticity and hydration of the molecule. A detailed comparison between different structures shows that the lectin molecule is sturdy and that the effect of changes in pH, ligand binding and environment on it is small. The region involving the curved front beta-sheet and the loops around the second hydrophobic core is comparatively rigid. The back beta-sheet involved in quaternary association, which exhibits considerable variability, is substantially flexible, as is the sugar-binding region. The numbers of invariant water molecules in the hydration shell are small and they are mainly involved in metal coordination or in stabilizing unusual structural features. Small consistent movements occur in the combining site upon sugar binding, although the site is essentially preformed.

Nuclear Magnetic Resonance Spectroscopy in the Study of Hemoglobin Cooperativity

Bacterial expression systems have greatly advanced NMR-based studies of Hb A and understanding of the relationship between structure and function in Hb. These advances result largely from the ability to isotopically label HbA for high-resolution heteronuclear NMR studies, literally opening up new dimensions in the NMR of hemoglobin, permitting extensive resonance assignments to be made, providing information about dynamics and hydrogen bonding, and yielding precise information regarding quaternary structure in solution. These expression systems also permit site-specific substitutions to be introduced to the Hb tetramer. NMR spectroscopy provides a powerful tool for quickly determining the effects of site-specific substitutions on the structure of the heme pocket, on site-specific distal ligand binding, on subunit interface structure and dynamics, and on the overall quaternary structure in solution.One area where NMR spectroscopy is likely to further advance the understanding of hemoglobin allostery is conformational dynamics. A few sites at the α 1 β 2 subunit interface appear to be quite dynamic over a range of time scales, and these dynamics appear to be coupled to distal ligand binding and quaternary structure. By fully characterizing the range of dynamics of different forms of Hb, NMR spectroscopy provides a picture of hemoglobin allostery that is considerably more detailed than that provided by the static structures obtained from X-ray crystallography. Another area where NMR spectroscopy is likely to advance the understanding of structure and function in hemoglobin is in determination of solution quaternary structure of the T state, using the residual dipolar coupling method. Residual dipolar coupling methods can determine whether deoxy-Hb A in solution conforms well to the crystallographically observed T-quaternary structure, or whether additional structures are also present, as in distally ligated HbCO A. Although hemoglobin has been studied with great intensity by a wide variety of methods over the last half-century, many of the molecular details that control allostery are still unknown. Experimental results on rHbs derived from NMR, equilibrium, and kinetic studies of ligand binding clearly indicate that hemoglobin is a rather flexible protein, and that its conformation can adapt to a variety of perturbations, for example, amino acid substitutions, binding of ligands, and allosteric effectors. These results also indicate that allosteric interactions in Hb A involve multiple pathways of signal transmission.

Low Resolution Structures of the Retinoid X Receptor DNA-binding and Ligand-binding Domains Revealed by Synchrotron X-ray Solution Scattering

Nuclear receptors are ligand-inducible transcription factors that share structurally related DNA-binding (DBD) and ligand-binding (LBD) domains. Biochemical and structural studies have revealed the modular nature of DBD and LBD. Nevertheless, the domains function in concert in vivo. While high-resolution crystal structures of nuclear receptor DBDs and LBDs are available, there are no x-ray structural studies of nuclear receptor proteins containing multiple domains. We report the solution structures of the human retinoid X receptor DBD-LBD (hRXRalphaDeltaAB) region. We obtained ab initio shapes of hRXRalphaDeltaAB dimer and tetramer to 3.3 and 1.7 nm resolutions, respectively, and established the position and orientation of the DBD and LBD by fitting atomic coordinates of hRXRalpha DBD and LBD. The dimer is U-shaped with DBDs spaced at approximately 2 nm in a head to head orientation forming an angle of about 10 degrees with respect to each other and with an extensive interface area provided by the LBD. The tetramer is a more elongated X-shaped molecule formed by two dimers in head to head arrangement in which the DBDs are extended from the structure and spaced at about 6 nm. The close proximity of DBDs in dimers may facilitate homodimer formation on DNA; however, for the homodimer to bind to a DNA element containing two directly repeated half-sites, one of the DBDs would need to rotate with respect to the other element. By contrast, the separation of DBDs in the tetramers may account for their decreased ability to recognize DNA.

New challenges in quantum chemistry: Quests for accurate calculations for large molecular systems.

As quantum chemistry plays a more and more central role in many complicated chemical problems, it has become necessary to obtain accurate results for large molecular systems. Conventional quantum chemistry methods are either too expensive to apply to large systems or too approximate for the results to be reliable, and they fail to satisfy this requirement. A variety of different approaches is being developed with the aim of achieving this goal: local correlation methods; divide-and-conquer methods; linear-scaling density functional methods based on the fast multipole and other approximations; effective potential methods; and hybrid methods. ONIOM (our N-layered integrated molecular orbital plus molecular mechanics method), developed by the authors, is a hybrid method in which a large molecular system is divided into onion-skin-like layers, and different quantum chemistry/molecular mechanics methods are used for different parts of the system; the results are combined to extrapolatively estimate the results of high-level calculation for the real system. Several applications of ONIOM will be discussed.

Aptamer Beacons for the Direct Detection of Proteins

We have designed a new class of molecules, which we term aptamer beacons, for detecting a wide range of ligands. Similar to molecular beacons, aptamer beacons can adopt two or more conformations, one of which allows ligand binding. A fluorescence-quenching pair is used to report changes in conformation induced by ligand binding. An anti-thrombin aptamer was engineered into an aptamer beacon by adding nucleotides to the 5′-end which are complementary to nucleotides at the 3′-end of the aptamer. In the absence of thrombin, the added nucleotides will form a duplex with the 3′-end, forcing the aptamer beacon into a stem-loop structure. In the presence of thrombin, the aptamer beacon forms the ligand-binding structure. This conformational change causes a change in the distance between a fluorophore attached to the 5′-end and a quencher attached to the 3′-end. Aptamer beacon can be a sensitive tool for detecting proteins and other chemical compounds.

Equilibrium Unfolding of Bombyx mori Glycyl-tRNA Synthetase

Unfolding of Bombyx mori glycyl-tRNA synthetase was examined by multiple spectroscopic techniques. Tryptophan fluorescence of wild type enzyme and an N-terminally truncated form (N55) increased at low concentrations of urea or guanidine-HCl followed by a reduction in intensity at intermediate denaturant concentrations; a transition at higher denaturant was detected as decreased fluorescence intensity and a red-shifted emission. Solute quenching of fluorescence indicated that tryptophans become progressively solvent-exposed during unfolding. Wild type enzyme had stronger negative CD bands between 220 and 230 nm than the mutant, indicative of greater alpha-helical content. Urea or guanidine-HCl caused a reduction in ellipticity at 222 nm at low denaturant concentration with the wild type enzyme, a transition that is absent in the mutant; both enzymes exhibited a cooperative transition at higher denaturant concentrations. Both enzymes dissociate to monomers in 1.5 m urea. Unfolding of wild type enzyme is described by a multistate unfolding and a parallel two state unfolding; the two-state component is absent in the mutant. Changes in spectral properties associated with unfolding were largely reversible after dilution to low denaturant. Unfolding of glycyl-tRNA synthetase is complex with a native state, a native-like monomer, partially unfolded states, and the unfolded state.

The Use of Circular Dichroism in the Investigation of Protein Structure and Function

Circular Dichroism (CD) relies on the differential absorption of left and right circularly polarised radiation by chromophores which either possess intrinsic chirality or are placed in chiral environments. Proteins possess a number of chromophores which can give rise to CD signals. In the far UV region (240-180 nm), which corresponds to peptide bond absorption, the CD spectrum can be analysed to give the content of regular secondary structural features such as alpha-helix and beta-sheet. The CD spectrum in the near UV region (320-260 nm) reflects the environments of the aromatic amino acid side chains and thus gives information about the tertiary structure of the protein. Other non-protein chromophores such as flavin and haem moieties can give rise to CD signals which depend on the precise environment of the chromophore concerned. Because of its relatively modest resource demands, CD has been used extensively to give useful information about protein structure, the extent and rate of structural changes and ligand binding. In the protein design field, CD is used to assess the structure and stability of the designed protein fragments. Studies of protein folding make extensive use of CD to examine the folding pathway; the technique has been especially important in characterising molten globule intermediates which may be involved in the folding process. CD is an extremely useful technique for assessing the structural integrity of membrane proteins during extraction and characterisation procedures. The interactions between chromophores can give rise to characteristic CD signals. This is well illustrated by the case of the light harvesting complex from photosynthetic bacteria, where the CD spectra can be analysed to indicate the extent of orbital overlap between the rings of bacteriochlorophyll molecules. It is therefore evident that CD is a versatile technique in structural biology, with an increasingly wide range of applications.

In vitro selection of functional nucleic acids.

In vitro selection allows rare functional RNA or DNA molecules to be isolated from pools of over 10(15) different sequences. This approach has been used to identify RNA and DNA ligands for numerous small molecules, and recent three-dimensional structure solutions have revealed the basis for ligand recognition in several cases. By selecting high-affinity and -specificity nucleic acid ligands for proteins, promising new therapeutic and diagnostic reagents have been identified. Selection experiments have also been carried out to identify ribozymes that catalyze a variety of chemical transformations, including RNA cleavage, ligation, and synthesis, as well as alkylation and acyl-transfer reactions and N-glycosidic and peptide bond formation. The existence of such RNA enzymes supports the notion that ribozymes could have directed a primitive metabolism before the evolution of protein synthesis. New in vitro protein selection techniques should allow for a direct comparison of the frequency of ligand binding and catalytic structures in pools of random sequence polynucleotides versus polypeptides.

Characterization of C3dg Binding to a Recess Formed Between Short Consensus Repeats 1 and 2 of Complement Receptor Type 2 (CR2; CD21)

To allow for a better characterization of the ligand binding structures of human complement receptor type 2 (CR2; CD21), we have established an IgG1 kappa mouse mAb, FE8, that interferes efficiently with binding of C3dg and EBV to CR2. In contrast to mAb OKB7, the only well-characterized mAb with similar specificity, mAb FE8 blocked binding of soluble C3dg or particles carrying multiple copies of surface-bound C3dg to CR2 or induced complete removal of these ligands from the receptor. In vitro EBV infection of B lymphocytes, on the other hand, was abrogated by mAbs FE8 and OKB7 with similar dose-response characteristics. As FE8 was shown to recognize a discontinuous epitope, a series of overlapping peptides derived from SCR1 and -2 and immobilized on cellulose was screened with FE8. The results suggest that up to five discontinuous sequences contributed to the epitope. The sequence 63-EYFNKYS-69, located between the two SCR units, reacted most intensively. Two other sequences, 16-YYSTPI-21 and 105-NGNKSVWCQANN-116, are located between Cys1 and Cys2 of SCR1 and around Cys3 of SCR2, respectively. Based on the solution structure for two factor H SCRs, a three-dimensional model of SCR1 and -2 was generated. The FE8 binding peptide sequences were located in relative proximity to each other, bounding the recess formed between SCR1 and -2. This potential of mAb FE8 is currently unique and may be exploited for interfering with conditions of unwanted recognition of C3dg-coated structures by the immune system.

Ligand binding and covalent structure of an oxygen-binding heme protein from Rhodobacter sphaeroides, a representative of a new structural family of c-type cytochromes.

The amino acid sequence of an oxygen-binding heme protein (SHP) from Rhodobacter sphaeroides has been determined. The cysteines, which bind the single heme group in the 112-residue protein, are located at positions 43 and 46. SHP is similar in size to the large membrane-bound form of the class I cytochrome c5 of Azotobacter vinelandii (116 residues) and in the location of the heme binding site at positions 48 and 51. Two extra cysteines in SHP (residues 89 and 97) are located in positions similar to those of cytochrome c5 (residues 98 and 101) and form a disulfide bridge in both proteins. In total, four regions of alpha-helix are predicted, covering 46% of the protein, which is comparable to that in other small cytochromes. SHP is thus distantly related to small class I c-type cytochromes but is representative of a distinct family by virtue of its high-spin nature, the lack of a strong sixth ligand, and its capacity to bind oxygen. Potentially, the most important characteristic of SHP is its ability to transiently bind oxygen during autoxidation, which occurs with a half-life of 3 min with a 4-fold excess of O2. SHP also binds carbon monoxide, azide, and cyanide. The kinetics of reduction by free flavins indicate that SHP is less reactive than other class I cytochromes c and that the heme is less accessible to solvent. There is localized positive charge (+3) at the site of reduction of SHP, although the overall protein charge is -2. This may account in part for the ability of SHP to bind anions.

Synthesis, ligand binding, and quantitative structure-activity relationship study of 3 beta-(4'-substituted phenyl)-2 beta-heterocyclic tropanes: evidence for an electrostatic interaction at the 2 beta-position.

A set of 3 beta-(4'-substituted phenyl)-2 beta-heterocyclic tropanes was designed, synthesized, and characterized. We discovered that these compounds can function as bioisosteric replacements for the corresponding WIN 35,065-2 analogs which possess a 2 beta-carbomethoxy group. Several of the compounds showed high affinity and selectivity for the dopamine transporter (DAT) relative to the serotonin and norepinephrine transporters. From the structure-activity relationship study, the 3 beta-(4'-chlorophenyl)-2 beta-(3'-phenylisoxazol-5-yl)tropane (5d) emerged as the most potent and selective compound. The binding data for 2 beta-heterocyclic tropanes were found to show a high correlation with molecular electrostatic potential (MEP) minima near one of the heteroatoms in the 2 beta-substituents. In contrast, low correlations were found for other MEP minima near the 2 beta-substituent as well as for calculated log P or substituent volume. These quantitative structure-activity relationship studies are consistent with an electrostatic contribution to the binding potency of these WIN 35,065-2 analogs at the DAT.