Defining the secondary structural requirements of a cocaine-binding aptamer by a thermodynamic and mutation study

Isothermal titration calorimetry (ITC) was used to measure the binding affinity and thermodynamics of a cocaine-binding aptamer as a function of pH and NaCl concentration. Tightest binding was achieved at a pH value of 7.4 and under conditions of no added NaCl. These data indicate that ionic interactions occur in the ligand binding mechanism. ITC was also used to measure the binding thermodynamics of a variety of sequence variants of the cocaine-binding aptamer that analyzed which regions and nucleotides of the aptamer are important for maintaining high-affinity binding. Individually, each of the three stems can be shortened, resulting in a reduced binding affinity. If all three stems are shortened, no binding occurs. If all three of the stems in the aptamer are lengthened by five base pairs ligand affinity increases. Changes in nucleotide identity at the three-way junction all decrease the affinity of the aptamer to cocaine. The greatest decrease in affinity results from changes that disrupt the GA base pairs and the identity of T19.     

Engineering a structure switching mechanism into a steroid-binding aptamer and hydrodynamic analysis of the ligand binding mechanism

The steroid binding mechanism of a DNA aptamer was studied using isothermal titration calorimetry (ITC), NMR spectroscopy, quasi-elastic light scattering (QELS), and small-angle X-ray spectroscopy (SAXS). Binding affinity determination of a series of steroid-binding aptamers derived from a parent cocaine-binding aptamer demonstrates that substituting a GA base pair with a GC base pair governs the switch in binding specificity from cocaine to the steroid deoxycholic acid (DCA). Binding of DCA to all aptamers is an enthalpically driven process with an unfavorable binding entropy. We engineered into the steroid-binding aptamer a ligand-induced folding mechanism by shortening the terminal stem by two base pairs. NMR methods were used to demonstrate that there is a transition from a state where base pairs are formed in one stem of the free aptamer, to where three stems are formed in the DCA-bound aptamer. The ability to generate a ligand-induced folding mechanism into a DNA aptamer architecture based on the three-way junction of the cocaine-binding aptamer opens the door to obtaining a series of aptamers all with ligand-induced folding mechanisms but triggered by different ligands. Hydrodynamic data from diffusion NMR spectroscopy, QELS, and SAXS show that for the aptamer with the full-length terminal stem there is a small amount of structure compaction with DCA binding. For ligand binding by the short terminal stem aptamer, we propose a binding mechanism where secondary structure forms upon DCA binding starting from a free structure where the aptamer exists in a compact form.     

Label-free detection of enhanced saccharide binding at pH 7.4 to nanoparticulate benzoboroxole based receptor units

Nanoparticles modified with either 6-amino-1-hydroxy-2,1-benzoxaborolane (3-aminobenzoboroxole) or 3-aminophenylboronic acid were prepared by nucleophilic substitution of a styrene-co-DVB-co-vinylbenzylchloride latex (25 nm). Isothermal titration calorimetry (ITC) was used as a label-free detection method for the analysis of the binding between monosaccharides and these two differently derivatized nanoparticle systems at pH 7.4. Because ITC reveals, thermodynamical parameters such as the changes in enthalpy ΔH, free energy ΔG, and entropy ΔS, possible explanations for the higher binding constants can be derived in terms of entropy and enthalpy changes. In case of the modified nanoparticles, the free energy of binding is dominated by the entropy term. This shows that interfacial effects, besides the intrinsic affinity, lead to a higher binding constant compared with the free ligand. The highest binding constant was found for fructose binding to the benzoboroxole modified nanoparticles: Its value of 1150 M(-1) is twice as high as for the free benzoboroxole and five times as high as with phenylboronic acid or 3-aminophenylboronic acid. In contrast to the binding of fructose to free boronic acids, which is an enthalpically driven process, the binding of fructose to the modified nanoparticles is dominated by the positive entropy term.     

Salt bridge exchange binding mechanism between streptavidin and its DNA aptamer – thermodynamics and spectroscopic evidences

Protein‐nucleic acids binding driven by electrostatic interactions typically are characterized by the release of counter ions, and the salt‐inhibited binding association constant (K_a) and the magnitude of exothermic binding enthalpy (ΔH). Here, we report a non‐classical thermodynamics of streptavidin (SA)–aptamer binding in NaCl (140–350 mM) solutions near room temperatures (23–27 °C). By using isothermal titration calorimetry (ITC) and circular dichroism (CD)/fluorescence spectroscopy, we found that the binding was enthalpy driven with a large entropy cost (ΔH ?20.58 kcal mol^?1, TΔS ?10.99 kcal mol^?1, and K_a 1.08 × 10⁷ M^?1 at 140 mM NaCl 25 °C). With the raise of salt concentrations, the ΔH became more exothermic, yet the K_a was almost unchanged (ΔH ?26.29 kcal mol^?1 and K_a 1.50 × 10⁷ M^?1 at 350 mM NaCl 25 °C). The data suggest that no counter Na⁺ was released in the binding. Spectroscopy data suggest that the binding, with a stoichiometry of 2, was accompanied with substantial conformational changes on SA, and the changes were insensitive to the variation of salt concentrations. To account for the non‐classical results, we propose a salt bridge exchange model. The intramolecular binding‐site salt bridge(s) of the free SA and the charged phosphate group of aptamers re‐organize to form the binding complex by forming a new intermolecular salt bridge(s). The salt bridge exchange binding process requires minimum amount of counter ions releasing but dehydration of the contacting surface of SA and the aptamer. The energy required for dehydration is reduced in the case of binding solution with higher salt concentration and account for the higher binding exothermic mainly. Copyright © 2013 John Wiley & Sons, Ltd.     

Evidence for a substrate specific and inhibitable drug efflux system in chloroquine resistant Plasmodium falciparum strains

The mechanism underpinning chloroquine drug resistance in the human malarial parasite Plasmodium falciparum remains controversial. By investigating the kinetics of chloroquine accumulation under varying-trans conditions, we recently presented evidence for a saturable and energy-dependent chloroquine efflux system present in chloroquine resistant P. falciparum strains. Here, we further characterize the putative chloroquine efflux system by investigating its substrate specificity using a broad range of different antimalarial drugs. Our data show that preloading cells with amodiaquine, primaquine, quinacrine, quinine, and quinidine stimulates labeled chloroquine accumulation under varying-trans conditions, while mefloquine, halofantrine, artemisinin, and pyrimethamine do not induce this effect. In the reverse of the varying-trans procedure, we show that preloaded cold chloroquine can stimulate quinine accumulation. On the basis of these findings, we propose that the putative chloroquine efflux system is capable of transporting, in addition to chloroquine, structurally related quinoline and methoxyacridine antimalarial drugs. Verapamil and the calcium/calmodulin antagonist W7 abrogate stimulated chloroquine accumulation and energy-dependent chloroquine extrusion. Our data are consistent with a substrate specific and inhibitible drug efflux system being present in chloroquine resistant P. falciparum strains.     

On the molecular basis of the high affinity binding of basic amino acids to LAOBP,a periplasmic binding protein from Salmonella typhimurium

The rational designing of binding abilities in proteins requires an understanding of the relationship between structure and thermodynamics. However, our knowledge of the molecular origin of high‐affinity binding of ligands to proteins is still limited; such is the case for l ‐lysine–l ‐arginine–l ‐ornithine periplasmic binding protein (LAOBP), a periplasmic binding protein from Salmonella typhimurium that binds to l ‐arginine, l ‐lysine, and l ‐ornithine with nanomolar affinity and to l ‐histidine with micromolar affinity. Structural studies indicate that ligand binding induces a large conformational change in LAOBP. In this work, we studied the thermodynamics of l ‐histidine and l ‐arginine binding to LAOBP by isothermal titration calorimetry. For both ligands, the affinity is enthalpically driven, with a binding ΔCp of ~?300 cal mol^?1 K^?1, most of which arises from the burial of protein nonpolar surfaces that accompanies the conformational change. Osmotic stress measurements revealed that several water molecules become sequestered upon complex formation. In addition, LAOBP prefers positively charged ligands in their side chain. An energetic analysis shows that the protein acquires a thermodynamically equivalent state with both ligands. The 1000‐fold higher affinity of LAOBP for l ‐arginine as compared with l ‐histidine is mainly of enthalpic origin and can be ascribed to the formation of an extra pair of hydrogen bonds. Periplasmic binding proteins have evolved diverse energetic strategies for ligand recognition. STM4351, another arginine binding protein from Salmonella, shows an entropy‐driven micromolar affinity toward l ‐arginine. In contrast, our data show that LAOBP achieves nanomolar affinity for the same ligand through enthalpy optimization. Copyright © 2015 John Wiley & Sons, Ltd.     

Recent trends and some applications of isothermal titration calorimetry in biotechnology

Isothermal titration calorimeters (ITCs) are thermodynamic instruments used for the determination of enthalpy changes in any physical/chemical reaction. This can be applied in various fields of biotechnology. This review explains ITC applications, especially in bioseparation, drug development and cell metabolism. In liquid chromatography, the separation/purification of specific proteins or polypeptides in a mixture is usually achieved by varying the adsorption affinities of the different proteins/polypeptides for the adsorbent under different mobile-phase conditions and temperatures. Using ITC analysis, the binding mechanism of proteins with adsorbent solid material is derived by elucidating enthalpy and entropy changes, which offer valuable guidelines for designing experimental conditions in chromatographic separation. The binding affinity of a drug with its target is studied by deriving binding enthalpy and binding entropy. To improve the binding affinity, suitable lead compounds for a drug can be identified and their affinity tested by ITC. Recently ITC has also been used in studying cell metabolism. The heat produced by animal cells in culture can be used as a primary indicator of the kinetics of cell metabolism, which provides key information for drug bioactivity and operation parameters for process cell culture.     

Blood schizontocidal activity of methylene blue in combination with antimalarials against Plasmodium falciparum

Methylene blue (MB) is the oldest synthetic antimalarial. It is not used anymore as antimalarial but should be reconsidered. For this purpose we have measured its impact on both chloroquine sensitive and resistant Plasmodium strains. We showed that around 5 nM of MB were able to inhibit 50% of the parasite growth in vitro and that late rings and early trophozoites were the most sensitive stages; while early rings, late trophozoites and schizonts were less sensitive. Drug interaction study following fractional inhibitory concentrations (FIC) method showed antagonism with amodiaquine, atovaquone, doxycycline, pyrimethamine; additivity with artemether, chloroquine, mefloquine, primaquine and synergy with quinine. These results confirmed the interest of MB that could be integrated in a new low cost antimalarial combination therapy.     

Plasmodium falciparum malaria: rosettes are disrupted by quinine, artemisinin, mefloquine, primaquine, pyrimethamine, chloroquine and proguanil.

An assay was developed measuring the disruption of rosettes between Plasmodium falciparuminfected (trophozoites) and uninfected erythrocytes by the antimalarial drugs quinine, artemisinin mefloquine, primaquine, pyrimethamine, chloroquine and proguanil. At 4 hr incubation rosettes were disrupted by all the drugs in a dose dependent manner. Artemisinin and quinine were the most effective anti-malarials at disrupting rosettes at their therapeutic concentrations with South African RSA 14, 15, 17 and The Gambian FCR-3 P. falciparum strains. The least effective drugs were proguanil and chloroquine. A combination of artemisinin and mefloquine was more effective than each drug alone. The combinations of pyrimethamine or primaquine, with quinine disrupted more rosettes than quinine alone. Quinine may be an effective drug in the treatment of severe malaria because the drug efficiently reduces the number of rosettes.     

Heat capacity changes associated with guanine quadruplex formation: an isothermal titration calorimetry study

This study addresses the temperature dependence of the enthalpy of formation for several unimolecular quadruplexes in the presence of excess monovalent salt. We examined a series of biologically significant guanine-rich DNA sequences: thrombin binding aptamer (TBA) (d(G(2)T(2)G(2)TGTG(2)T(2)G(2)), PS2.M, a catalytically active aptamer (d(GTG(3)TAG(3)CG(3)T(2)G(2))), and the human telomere repeat (HT) (d(AG(3)(T(2)AG(3))(3))). Using CD spectra and UV melting, we confirmed the presence of quadruplex structures and established the temperature range in which quadruplex conformation is stable. We then performed ITC experiments, adding DNA to a solution containing excess NaCl or KCl. In this approach, only several additions are made, and only the enthalpy of quadruplex formation is measured. This measurement was repeated at different temperatures to determine the temperature dependence of the enthalpy change accompanying quadruplex formation. To control for the effect of nonspecific salt interactions during DNA folding, we repeated the experiment by replacing the quadruplex-forming sequences with a similar but nonfolding sequence. Dilution enthalpies were also subtracted to obtain the final enthalpy value involving only the quadruplex folding process. For all sequences studied, quadruplex formation was exothermic but with an increasing magnitude with increasing temperature. These results are discussed in terms of the change in heat capacity associated with quadruplex formation.     

Structural parameterization of the binding enthalpy of small ligands

A major goal in ligand and drug design is the optimization of the binding affinity of selected lead molecules. However, the binding affinity is defined by the free energy of binding, which, in turn, is determined by the enthalpy and entropy changes. Because the binding enthalpy is the term that predominantly reflects the strength of the interactions of the ligand with its target relative to those with the solvent, it is desirable to develop ways of predicting enthalpy changes from structural considerations. The application of structure/enthalpy correlations derived from protein stability data has yielded inconsistent results when applied to small ligands of pharmaceutical interest (MW < 800). Here we present a first attempt at an empirical parameterization of the binding enthalpy for small ligands in terms of structural information. We find that at least three terms need to be considered: (1) the intrinsic enthalpy change that reflects the nature of the interactions between ligand, target, and solvent; (2) the enthalpy associated with any possible conformational change in the protein or ligand upon binding; and, (3) the enthalpy associated with protonation/deprotonation events, if present. As in the case of protein stability, the intrinsic binding enthalpy scales with changes in solvent accessible surface areas. However, an accurate estimation of the intrinsic binding enthalpy requires explicit consideration of long-lived water molecules at the binding interface. The best statistical structure/enthalpy correlation is obtained when buried water molecules within 5-7 A of the ligand are included in the calculations. For all seven protein systems considered (HIV-1 protease, dihydrodipicolinate reductase, Rnase T1, streptavidin, pp60c-Src SH2 domain, Hsp90 molecular chaperone, and bovine beta-trypsin) the binding enthalpy of 25 small molecular weight peptide and nonpeptide ligands can be accounted for with a standard error of +/-0.3 kcal x mol(-1).     

Impact of protein and ligand impurities on ITC-derived protein–ligand thermodynamics

Background

The thermodynamic characterization of protein–ligand interactions by isothermal titration calorimetry (ITC) is a powerful tool in drug design, giving valuable insight into the interaction driving forces. ITC is thought to require protein and ligand solutions of high quality, meaning both the absence of contaminants as well as accurately determined concentrations.

Methods

Ligands synthesized to deviating purity and protein of different pureness were titrated by ITC. Data curation was attempted also considering information from analytical techniques to correct stoichiometry.

Results and conclusions

We used trypsin and tRNA-guanine transglycosylase (TGT), together with high affinity ligands to investigate the effect of errors in protein concentration as well as the impact of ligand impurities on the apparent thermodynamics. We found that errors in protein concentration did not change the thermodynamic properties obtained significantly. However, most ligand impurities led to pronounced changes in binding enthalpy. If protein binding of the respective impurity is not expected, the actual ligand concentration was corrected for and the thus revised data compared to thermodynamic properties obtained with the respective pure ligand. Even in these cases, we observed differences in binding enthalpy of about 4 kJ ⋅ mol^− 1, which is considered significant.

General significance

Our results indicate that ligand purity is the critical parameter to monitor if accurate thermodynamic data of a protein–ligand complex are to be recorded. Furthermore, artificially changing fitting parameters to obtain a sound interaction stoichiometry in the presence of uncharacterized ligand impurities may lead to thermodynamic parameters significantly deviating from the accurate thermodynamic signature.     

The differences in binding 12‐carbon aliphatic ligands by bovine β‐lactoglobulin isoform A and B studied by isothermal titration calorimetry and X‐ray crystallography

Isoforms A (LGB‐A) and B (LGB‐B) of bovine lactoglobulin, the milk protein, differ in positions 64 (D?G) and 118 (V?A). Interactions of LGB‐A and LGB‐B with sodium dodecyl sulfate (SDS), dodecyltrimethylammonium chloride (DTAC) and lauric acid (LA), 12‐carbon ligands possessing differently charged polar groups, were investigated using isothermal titration calorimetry and X‐ray crystallography, to study the proton linkage phenomenon and to distinguish between effects related to different isoforms and different ligand properties. The determined values of ΔS and ΔH revealed that for all ligands, binding is entropically driven. The contribution from enthalpy change is lower and shows strong dependence on type of buffer that indicates proton release from the protein varying with protein isoform and ligand type and involvement of LA and Asp64 (in isoform A) in this process. The ligand affinities for both isoforms were arranged in the same order, DTAC < LA < SDS, and were systematically lower for variant B. The entropy change of the complexation process was always higher for isoform A, but these values were compensated by changes in enthalpy, resulting in almost identical ΔG for complexes of both isoforms. The determined crystal structures showed that substitution in positions 64 and 118 did not influence the overall structure of LGB complexes. The chemical character of the ligand polar group did not affect the position of its aliphatic chain in protein β‐barrel, indicating a major role of hydrophobic interactions in ligand binding that prevailed even with the repulsion between positively charged DTAC and lysine residues located at binding site entrance. Copyright © 2013 John Wiley & Sons, Ltd.     

Isothermal titration calorimetry for chiral chemistry

One of the most powerful techniques that are currently available to measure thermodynamic parameters such as enthalpy (ΔH), Gibbs free energy (ΔG), entropy changes (ΔS), and binding affinity in chemical reactions is isothermal titration calorimetry (ITC). Recent advances in instrumentation have facilitated the development of ITC as a very essential analytical tool in biology and chemistry. In this article, we will focus on a review of the literature on the application of ITC for the study of chiral systems and chiral interactions. We present studies in which the ITC technique is used to study chiral interactions, for instance in chiral solutions, chiral organometallic complexes, guest‐host chiral binding interactions, and biological macromolecules. Finally, we put strong emphasis on the most recent application of ITC for the study of chirality in nanosystems and at the nanoscale.     

Contribution of the trifluoroacetyl group in the thermodynamics of antigen–antibody binding

We analyzed the binding of the 7C8 antibody to the chloramphenicol phosphonate antigens—one containing a trifluoroacetyl group (CP‐F) and the other containing an acetyl group (CP‐H)—by using isothermal titration calorimetry (ITC). The thermodynamic difference due to the substitution of F by H was evaluated using free energy calculations based on molecular dynamics (MD) simulations. We have previously shown that another antibody, namely, 6D9, binds more weakly to CP‐H than to CP‐F, mainly due to the different hydration free energies of the dissociated state and not due to the unfavorable hydrophobic interactions with the antibody in the bound state. Unlike in the binding of the trifluoroacetyl group with 6D9, in its binding with 7C8, it is exposed to the solvent, as seen in the crystal structure of the complex of 7C8 with CP‐F. The thermodynamic analysis performed in this study showed that the binding affinity of 7C8 for CP‐H is similar to that for CP‐F, but this binding to CP‐H is accompanied with less favorable enthalpy and more favorable entropy changes. The free energy calculations indicated that, upon the substitution of F by H, enthalpy and entropy changes in the associated and dissociated states were decreased, but the magnitude of enthalpy and entropy changes in the dissociated state was larger than that in the associated state. The differences in binding free energy, enthalpy, and entropy changes determined by the free energy calculations for the substitution of F by H are in good agreement with the experimental results. Copyright © 2009 John Wiley & Sons, Ltd.     

Chiral resolution of some antimalarial agents by sub- and supercritical fluid chromatography on an (S)-naphthylurea stationary phase

The behavior of mefloquine, halofantrine, enpiroline, quinine, quinidine, chloroquine and primaquine is studied by subcritical fluid chromatography on a (S)-naphthylurea column (250 mm × 4.6 mm ID) with a subcritical mobile phase composed of carbon dioxide, methanol and triethylamine (flow rate of 3 ml/min). Except for primaquine and chloroquine, each enantiomer was separated at a temperature between 40 and 60°C, and at a pressure below 15 MPa. A 98/2, v/v CO₂/methanol 0.1% triethylamine mixture allowed the separation of halofantrine enantiomers while the enantiomers of the more polar metabolite (N-desbutylhalofantrine) were separated with a 80–20 v/v mixture as used for mefloquine, enpiroline, quinine and quinidine. The influence of temperature, pressure and of the nature of the mobile phase is discussed. © 1993 Wiley-Liss, Inc.     

Thermodynamics of ligand binding to a heterogeneous RNA population in the malachite green aptamer

The malachite green aptamer binds two closely related ligands, malachite green (MG) and tetramethylrosamine (TMR), with nearly equal affinity. The MG ligand consists of three phenyl rings emanating from a central carbon, while TMR has two of the three rings connected by an ether linkage. The binding pockets for MG and TMR in the aptamer, known from high-resolution structures, differ only in the conformation of a few nucleotides. Herein, we applied isothermal titration calorimetry (ITC) to compare the thermodynamics of binding of MG and TMR to the aptamer. Binding heat capacities were obtained from ITC titrations over the temperature range of 15-60 °C. Two temperature regimes were found for MG binding: one from 15 to 45 °C where MG bound with a large negative heat capacity and an apparent stoichiometry (n) of ~0.4 and another from 50 to 60 °C where MG bound with a positive heat capacity and an n of ~1.1. The binding of TMR, on the other hand, revealed only one temperature regime for binding, with a more modest negative heat capacity and an n of ~1.2. The large difference in heat capacity between the two ligands suggests that significantly more conformational rearrangement occurs upon the binding of MG than that of TMR, which is consistent with differences in solvent accessible surface area calculated for available ligand-bound structures. Lastly, we note that the binding stoichiometry of MG was improved not only by an increase in the temperature but also by a decrease in the concentration of Mg(2+) or an increase in the time between ITC injections. These studies suggest that binding of a dynamical ligand to a functional RNA requires the RNA itself to have significant dynamics.     

Determination of substrate binding energies in individual subsites of a family 18 chitinase

Thermodynamic parameters for binding of N-acetylglucosamine (GlcNAc) oligomers to a family 18 chitinase, ChiB of Serratia marcescens, have been determined using isothermal titration calorimetry. Binding studies with oligomers of different lengths showed that binding to subsites −2 and +1 is driven by a favorable enthalpy change, while binding to the two other most important subsites, +2 and +3, is driven by entropy with unfavorable enthalpy. These remarkable unfavorable enthalpy changes are most likely due to favorable enzyme-substrate interactions being offset by unfavorable enthalpic effects of the conformational changes that accompany substrate-binding.     

Inhibition of β-Amylase Activity by Calcium,Magnesium and Zinc Ions Determined by Spectrophotometry and Isothermal Titration Calorimetry

The inhibition effect of metal ions on beta amylase activity was studied. The inhibitor-binding constant (Ki) was determined by spectrophotometric and isothermal titration calorimetric (ITC) methods. The binding of calcium, magnesium and zinc ion as inhibitors at the active site of barley beta amylase was studied at pH=4.8 (sodium acetate 16?mM) and T=300?K. The Ki and enthalpy of binding for calcium (13.4, 13.1?mM and -14.3?kJ/mol), magnesium (18.6, 17.8?mM and -17.7?kJ/mol) and zinc (17.5, 17.7?mM and -20.0?kJ/mol) were found by spectrophotometric and ITC methods respectively.     

Three-dimensional structure-activity relationship modeling of cocaine binding to two monoclonal antibodies by comparative molecular field analysis

Successful immunotherapy of cocaine addiction and overdoses requires cocaine-binding antibodies with specific properties, such as high affinity and selectivity for cocaine. We have determined the affinities of two cocaine-binding murine monoclonal antibodies (mAb: clones 3P1A6 and MM0240PA) for cocaine and its metabolites by [3H]-radioligand binding assays. mAb 3P1A6 (K(d) = 0.22 nM) displayed a 50-fold higher affinity for cocaine than mAb MM0240PA (K(d) = 11 nM) and also had a greater specificity for cocaine. For the systematic exploration of both antibodies' binding specificities, we used a set of approximately 35 cocaine analogues as structural probes by determining their relative binding affinities (RBAs) using an enzyme-linked immunosorbent competition assay. Three-dimensional quantitative structure-activity relationship (3D-QSAR) models on the basis of comparative molecular field analysis (CoMFA) techniques correlated the binding data with structural features of the ligands. The analysis indicated that despite the mAbs' differing specificities for cocaine, the relative contributions of the steric (approximately 80%) and electrostatic (approximately 20%) field interactions to ligand-binding were similar. Generated three-dimensional CoMFA contour plots then located the specific regions about cocaine where the ligand/receptor interactions occurred. While the overall binding patterns of the two mAbs had many features in common, distinct differences were observed about the phenyl ring and the methylester group of cocaine. Furthermore, using previously published data, a 3D-QSAR model was developed for cocaine binding to the dopamine reuptake transporter (DAT) that was compared to the mAb models. Although the relative steric and electrostatic field contributions were similar to those of the mAbs, the DAT cocaine-binding site showed a preference for negatively charged ligands. Besides establishing molecular level insight into the interactions that govern cocaine binding specificity by biopolymers, the three-dimensional images obtained reflect the properties of the mAbs binding pockets and provide the initial information needed for the possible design of novel antibodies with properties optimized for immunotherapy.