The interaction with DNA of unfused aromatic systems containing terminal piperazino substituents. Intercalation and groove-binding

A number of unfused tricyclic aromatic intercalators have shown excellent activity as amplifiers of the anticancer activity of the bleomycins and the 4',6-diphenylpyrimidines, 2a and 2b, with terminal basic functions (4-methylpiperazino groups) have been synthesized to test the structural requirements for amplifier-DNA interactions. The terminal piperazine rings are bulky, have limited flexibility, and are twisted out of the phenyl ring plane in both 2a and 2b. With 2a the pyrimidine is unsubstituted at position 5 and the conformation predicted by molecular mechanics calculations has a 25-30 degrees twist between the phenyl and pyrimidine ring planes. With 2b the 5-position is substituted with a methyl group and this causes a larger twist angle (50-60 degrees) between the phenyl and pyrimidine planes. These conformational variations lead to markedly different DNA interactions for 2a and 2b. Absorption, CD and NMR spectral, viscometric, flow dichroism and kinetics results indicate that 2a binds strongly to DNA by intercalation while 2b binds more weakly in a groove complex. The general structure and conformation of 2a, a slightly twisted, unfused-aromatic system with terminal piperazino groups is more similar to groove-binding agents such as Hoechst 33258 than to intercalators. The fact that 2a forms a strong intercalation complex with DNA is unusual but in agreement with studies on other amplifiers of anticancer drug action. Molecular modeling studies provide a second unusual feature of the 2a intercalation complex. While most well-characterized intercalators bind with their bulky and/or cationic substitutents in the DNA minor groove, the cationic piperazino groups of 2a are too large to bind in the minor groove in an intercalation complex but can form strong interactions with DNA in the major groove. The tricyclic aromatic ring system of 2a stacks well with adjacent base-pairs in the major-groove complex and the piperazino groups have good electrostatic and van der Waals interactions with the DNA backbone.     

The influence of N-dialkyl and other cationic substituents on DNA intercalation and genotoxicity

DNA intercalation by small chemical molecules can result in frameshift mutagenesis and chromosomal breakage. With evidence mounting that broadly diverse structures are capable of intercalating between DNA base pairs, it becomes important to better define those structural features that enhance intercalation strength and those that confer genotoxicity particularly among those intercalators that do not have the classical planar tricyclic fused ring structure. A chemical substituent that is present on many pharmaceutical and other biologically active molecules is the N-dialkyl group. In the present study, we investigate if and how the presence of an aromatic N-dialkyl or other cationic group affects the genotoxicity and DNA intercalation ability of 26 selected acridines, phenothiazines, benzophenones, triphenylethylenes and other classes of molecules. The data were obtained from the literature, from experiments using a cell-based DNA intercalation assay, and from modeling studies using a three-dimensional computational DNA docking program. It is demonstrated that cationic substitution can enhance both genotoxicity and electrostatic interactions within a chemical/DNA intercalation complex.     

Viscometric analysis of the interaction of bisphenanthridinium compounds with closed circular supercoiled and linear DNA.

The interaction with closed circular supercoiled and linear DNA of bisphenanthridinium compounds substituted through both the meta and para positions of the 6-phenyl group, along with appropriate monomer intercalators as controls, has been investigated by viscometric titration. When CPK models for the phenanthridinium rings of the three bis-compounds are oriented in a parallel manner as a model for intercalation, their ring plane to ring plane distances are approximately 7 to 8 A (SR 2430), 11 A (SR 2193), and 15 A (SR 2166). In SR 2430 the two phenanthridines are linked through the para positions of the 6-phenyl group; this chain allows intercalation of the two rings at adjacent binding sites in DNA, but is not long enough to accommodate an excluded site. The viscometric titrations with both superhelical and linear DNA clearly indicate that SR 2430 gives results close to those of the monomer control compounds while SR 2193 and SR 2166 have approximately twice the unwinding angle and DNA length increase on binding to DNA as the monomer compounds. These phenanthridinium compounds, therefore, are capable of bisintercalation only if their linking groups are of sufficient length to allow an excluded binding site between base pairs. This conclusion is supported by DNA thermal denaturation experiments in the presence of these compounds.     

DNA-drug recognition and effects on topoisomerase II-mediated cytotoxicity. A three-mode binding model for ellipticine derivatives.

Cytotoxic effects and topoisomerase II-mediated DNA breaks induced in vitro by ellipticine derivatives were examined in connection with 1H NMR and circular dichroism (CD) studies on molecular structures and interactions of drugs with DNA. The compounds included four 9-hydroxyellipticine and two 7-hydroxyisoellipticine derivatives. Structure-activity relationships indicated that a change in nitrogen atom position in the pyridinic ring greatly affected drug effects both on topoisomerase II action and cytotoxicity to L1210 cells. The four 9-hydroxyellipticine derivatives yielded bell-shaped curves in in vitro topoisomerase II-mediated DNA break assays, whereas the two 7-hydroxyisoellipticine derivatives demonstrated an almost linear increase at the same concentration (0-10 microM). In both cases, the intensity of cleavage was modulated by the position and the degree of methylation on the pyridinic ring, and results were correlated with cytotoxic activity expressed as the in vitro ID50 values for L1210 leukemia cells. 1H NMR experiments performed on free drug molecules in solution revealed that the two protons (alpha and beta) contiguous to the biologically important hydroxy group were sensitive to changes in electron distribution produced by the distant chemical modifications and methylations of the pyridinic ring. A linear relationship was observed between the differences in chemical shifts of alpha and beta protons (delta delta alpha-beta) versus ID50 values. CD experiments indicated that, at weak ionic strength I = 0.02 and at pH 7, drugs interact with the poly[d(A-T)] duplex according to a "three-mode binding model" which is governed by the drug structure and the drug to DNA ratio. The intercalation mode was related to the induction of topoisomerase II-mediated DNA cleavage, while the external binding mode consecutive to intercalation was related to cleavage suppression. These two modes concerned the good intercalators 9-hydroxyellipticines. The third was found for the weak intercalators 7-hydroxyisoellipticines and was characterized by self-stacked molecules bound "outside" DNA, presumably in the minor groove. Ligands either could be intercalated partially or linked at the edge of bases with a small number of molecules filling intercalation sites, for the second alternative. In addition to having different binding modes, 9-hydroxyellipticines were better inducers of DNA distortions than 7-hydroxyisoellipticines. The incidence of the drug binding modes on DNA-topoisomerase II recognition was discussed in connection with the in vitro cytotoxic activity exhibited by the drugs.     

Electrostatic factors in DNA intercalation

The factors that determine the binding of a chromophore between the base pairs in DNA intercalation complexes are dissected. The electrostatic potential in the intercalation plane is calculated using an accurate ab initio based distributed multipole electrostatic model for a range of intercalation sites, involving different sequences of base pairs and relative twist angles. There will be a significant electrostatic contribution to the binding energy for chromophores with a predominantly positive electrostatic potential, but this varies significantly with sequence, and somewhat with twist angle. The usefulness of these potential maps for understanding the binding of intercalators is explored by calculating the electrostatic binding energy for 9-aminoacridine, ethidium, and daunomycin in a variety of model binding sites. The electrostatic forces play a major role in the positioning of an intercalating 9-aminoacridine and a significant stabilizing role in the binding of ethidium in its sterically constrained position, but the intercalation of daunomycin is determined by the side-chain binding. Sequence preferences are likely to be determined by a complex and subtle mixture of effects, with electrostatics being just one component. The electrostatic binding energy is also unlikely to be a major determinant of the twist angle, as its variation with angle is modest for most intercalation sites. Overall, the electrostatic potential maps give guidance on how positively charged chromophores can be chemically adapted by heteroatomic substitution to optimise their binding.     

DAPI (4',6-diamidino-2-phenylindole) binds differently to DNA and RNA: minor-groove binding at AT sites and intercalation at AU sites.

The interaction of DAPI and propidium with RNA (polyA.polyU) and corresponding DNA (polydA.polydT) sequences has been compared by spectroscopic, kinetic, viscometric, Tm, and molecular modeling methods. Spectral changes of propidium are similar on binding to the AT and AU sequences but are significantly different for binding of DAPI. Spectral changes for DAPI with the DNA sequence are consistent with the expected groove-binding mode. All spectral changes for complexes of propidium with RNA and DNA and for DAPI with RNA, however, are consistent with an intercalation binding mode. When complexed with RNA, for example, DAPI aromatic protons signals shift significantly upfield, and the DAPI UV-visible spectrum shows significantly larger changes than when complexed with DNA. Slopes of log kd (dissociation rate constants) versus-log [Na+] plots are similar for complexes of propidium with RNA and DNA and for the DAPI-RNA complex and are in the range expected for an intercalation complex. The slope for the DAPI-DNA complex, however, is much larger and is in the range expected for a groove-binding complex. Association kinetics results also support an intercalation binding mode for the DAPI-RNA complex. The viscosity of polyA.polyU solutions increases significantly on addition of both propidium and DAPI, again in agreement with an intercalation binding mode for both molecules with RNA. Molecular modeling studies completely support the experimental findings and indicate that DAPI forms a very favorable intercalation complex with RNA. DAPI also forms a very stable complex in the minor groove of AT sequences of DNA, but the stabilizing interactions are considerably reduced in the wide, shallow minor groove of RNA. Modeling studies,thus,indicate that DAPI interaction energetics are more favorable for minor-groove binding in AT sequences but are more favorable for interaction in RNA.     

Neocarzinostatin chromophore binds to deoxyribonucleic acid by intercalation

The nonprotein chromophore of neocarzinostatin was found to share many of the characteristics of classical intercalators in its interaction with DNA. Viscosity studies with PM2 DNA indicated that the DNA helix unwinding induced by the chromophore was 0.82 times that of ethidium or 21 degrees. Electric dichroism of the chromophore--DNA complex showed that each bound chromophore molecule lengthened DNA by 3.3 A and that absorbance transitions of the chromophore at 315--385 nm were oriented approximately parallel to DNA bases, as expected for an intercalated aromatic ring. Binding to DNA induced strong hypochromicity and a pronounced red shift in the absorbance spectrum of the chromophore. Spectrophotometric titrations suggested at least two types of chromophore binding sites on DNA; one type of site was saturated at rb = 0.125 chromophore molecule/nucleotide, but binding to additional sites continued to at least rb = 0.3. These physical--chemical studies were performed at pH 4--5 in order to keep the chromophore stable, but chromophore bound to an excess of DNA at pH 7 showed a stable absorbance spectrum identical with that seen at pH 4--5, suggesting that a similar type of binding occurs at neutral pH. Chromophore which had spontaneously degraded in pH 8 buffer did not bind to DNA at all, as judged by absorbance spectroscopy. The degree of protection afforded by DNA against spontaneous chromophore degradation implied a dissociation constant of approximately 5 microM for the DNA--chromophore complex at neutral pH and physiological ionic strength. Supercoiled DNA was nearly twice as effective as relaxed DNA in protecting chromophore from degradation, providing additional evidence for intercalation at neutral pH. Comparison of absorbance, fluorescence, and dichroism spectra suggests that the naphthalene ring system is the intercalating moiety.     

Orientation and linear dichroism characteristics of porphyrin-DNA complexes

The linear dichroism spectra of complexes of tetrakis(N-methyl-4-pyridinio)prophine (H2TMpyP) and its zinc(II) derivative (ZnTMpyP) with DNA oriented in a flow gradient have been investigated. The dichroism of H2TMpyP determined within the Soret band and the Qy band system is consistent with an intercalative conformation in which the plane of the porphyrin ring system is nearly parallel to the planes of the DNA bases. In the case of ZnTMpyP on the other hand, the porphyrin ring system is inclined at angles of 62-67 degrees with respect to the axis of the DNA helix. The pyridyl groups in both cases are characterized by a low degree of orientation with respect to the axis of the helix. In contrast to H2TMpyP which does not significantly affect the degree of alignment of the DNA in the flow gradient, the binding of ZnTMpyP causes a significant decrease (about 50% for a base pair/ZnTMpyP ratio of 20) in the intrinsic dichroism at 260 nm due to the oriented DNA bases; the binding of ZnTMpyP to DNA either gives rise to regions of higher flexibility or causes bends or kinks at the binding sites. Increasing the ionic strength has little influence on the linear dichroism of the ZnTMpyP-DNA complexes, but the number of molecules bound at intercalation sites diminishes in the case of the H2TMpyP-DNA complexes; the accompanying changes in the linear dichroism characteristics suggest that external H2TMpyP complexes are formed at the expense of intercalation complexes.(ABSTRACT TRUNCATED AT 250 WORDS)     

Torsional sensing of small-molecule binding using magnetic tweezers

DNA-binding small molecules are widespread in the cell and heavily used in biological applications. Here, we use magnetic tweezers, which control the force and torque applied to single DNAs, to study three small molecules: ethidium bromide (EtBr), a well-known intercalator; netropsin, a minor-groove binding anti-microbial drug; and topotecan, a clinically used anti-tumor drug. In the low-force limit in which biologically relevant torques can be accessed (<10 pN), we show that ethidium intercalation lengthens DNA ∼1.5-fold and decreases the persistence length, from which we extract binding constants. Using our control of supercoiling, we measure the decrease in DNA twist per intercalation to be 27.3 ± 1° and demonstrate that ethidium binding delays the accumulation of torsional stress in DNA, likely via direct reduction of the torsional modulus and torque-dependent binding. Furthermore, we observe that EtBr stabilizes the DNA duplex in regimes where bare DNA undergoes structural transitions. In contrast, minor groove binding by netropsin affects neither the contour nor persistence length significantly, yet increases the twist per base of DNA. Finally, we show that topotecan binding has consequences similar to those of EtBr, providing evidence for an intercalative binding mode. These insights into the torsional consequences of ligand binding can help elucidate the effects of small-molecule drugs in the cellular environment.     

Different Binding Mode in AT and GC Sequences for Unfused-Aromatic Dications

Abstract

We have previously synthesized a 2,5-diphenylfuranamidine dication (4) and presented evidence that this compound binds to AT sequences in DNA by a minor-groove interaction mode but binds to GC sequences by intercalation (1,2). To probe these sequence-dependent binding modes in more detail, and particularly to obtain additional evidence for the binding mode in GC rich sequences, we have synthesized and studied the DNA complexes of 1–3 which have the furan ring of 4 replaced by 2,6-substituted pyridine (1), pyrimidine (2), or triazine (3) ring systems. The three compounds with a six-membered central ring system bind to AT DNA sequences more weakly than the furan compound, but retain the minor-groove binding mode. The pyridine and pyrimidine derivatives bind to GC sequences of DNA more strongly than the furan, but the triazine derivative binds more weakly. The aromatic proton signals of 1–3, as previously observed with 4 shift upfield by approximately 0.5 ppm or greater on complex formation with polyd(G-C)₂. This and other spectroscopic as well as viscosity and kinetics results indicate that 1–4 bind to GC sites in DNA by intercalation. A nonclassical intercalation model, with the twisted-unfused, aromatic ring system intercalated into an intercalation site of matching structure can explain all of our and the literature results for the GC binding mode of these unfused, aromatic compounds.     

Effect of intercalators on the properties of liquid-crystalline dispersions of nucleic acid-chitosan complex

Molecules of deoxyribonucleic acid and synthetic polydeoxyribonucleotides (NA) in the particles of liquid-crystalline dispersions resulting from interaction with chitosan are accessible to interaction with intercalators. The intercalation is accompanied by alteration in the direction of spatial twist of cholesterics of NA-chitosan complexes. This effect is absent in the case of "classical" cholesterics produced from NA molecules via phase exclusion, i.e., the cholesteric structure of NA-chitosan complex is very "labile" as distinct from "classical" cholesteric NA.     

Computational DNA binding studies of (–)-epigallocatechin-3-gallate

The catechin family of molecules that are present in the leaves of green tea has been under investigation since the antioxidant and anti-inflammatory properties of tea were discovered. Among multiple proposed therapeutic targets of these molecules, the direct interaction with nucleic acids has been proposed and experimentally observed but without clear knowledge about the potential binding modes between these ligands and DNA. One of these catechin structures, (–)-epigallocatechin gallate (EGCG), has three aromatic rings that could interact with double-stranded DNA via terminal base-pair stacking, intercalation, or through groove binding. Using enhanced sampling techniques and molecular dynamics simulations, we have found a stable complex between the EGCG ligand and DNA through intercalation of the trihydroxybenzoate aromatic ring and an ApC step. Moreover, we have calculated the absorption spectra of four possible binding modes and compared these to absorption profiles reported in the literature, and explored the possible DNA sequence preference for the EGCG ligand to bind. Our results suggest that an intercalative mode of interaction through the major groove is possible between the EGCG ligands and DNA with apparently very little DNA sequence selectivity.     

Harnessing DNA intercalation

Numerous small molecules are known to bind to DNA through base pair intercalation. Fluorescent dyes commonly used for nucleic acid staining, such as ethidium, are familiar examples. Biological and physical studies of DNA intercalation have historically been motivated by mutation and drug discovery research. However, this same mode of binding is now being harnessed for the creation of novel molecular assemblies. Recent studies have used DNA scaffolds and intercalators to construct supramolecular assemblies that function as fluorescent 'nanotags' for cell labeling. Other studies have demonstrated how intercalators can be used to promote the formation of otherwise unstable nucleic acid assemblies. These applications illustrate how intercalators can be used to facilitate and expand DNA-based nanotechnology.     

The interaction of ellipticine derivatives with nucleic acids studied by optical and 1H-nmr spectroscopy: Effect of size of the heterocyclic ring system

The DNA interaction of derivatives of ellipticine with heterocyclic ring systems with three, four, or five rings and a dimethylaminoethyl side chain was studied. Optical spectroscopy of drug complexes with calf thymus DNA, poly [(dA-dT) · (dA-dT)], or poly [(dG-dC) · (dG-dC)] showed a 10 nm bathochromic shift of the light absorption bands of the pentacyclic and tetracyclic compounds upon binding to the nucleic acids, which indicates binding by intercalation. For the tricyclic compound a smaller shift of 1–3 nm was observed upon binding to the nucleic acids. Flow linear dichroism studies show that the geometry of all complexes is consistent with intercalation of the ring system, except for the DNA and poly [(dG-dC) · (dG-dC)] complexes of the tricyclic compound, where the average angle between the drug molecular plane and the DNA helix axis was found to be 65°. One-dimensional ¹H-nmr spectroscopy was used to study complexes between d(CGCGATCGCG)₂ and the tricyclic and pentacyclic compounds. The results on the pentacyclic compound show nonselective broadening due to intermediate chemical exchange of most oligonucleotide resonances upon drug binding. The imino proton resonances are in slow chemical exchange, and new resonances with upfield shifts approaching 1 ppm appear upon drug binding, which supports intercalative binding of the pentacyclic compound. The results on the tricyclic compound show more rapid binding kinetics and very selective broadening of resonances. The data suggest that the tricyclic compound is in an equilibrium between intercalation and minor groove binding, with a preference to bind close to the AT base pairs with the side chain residing in the minor groove. © 1994 John Wiley & Sons, Inc.     

Orientation of the bases of single-stranded DNA and polynucleotides in complexes formed with the gene 32 protein of bacteriophage T4. A linear dichroism study

Linear dichroism measurements were performed in the wavelength region 250 to 350 nm on complexes between the single-stranded DNA binding protein of bacteriophage T4 (gp32) and single-stranded DNA and a variety of homopolynucleotides in compressed polyacrylamide gels. The complexes appeared to orient well, giving rise to linear dichroism spectra that showed contributions from both the protein aromatic residues and the bases of the polynucleotides. In most cases the protein contribution appeared to be very similar, and the linear dichroism of the bases could be explained by similar orientations of the bases for most of the complexes. Assuming a similar, regular structure for most of the polynucleotides in complex, only a limited set of combinations of tilt and twist angles can explain the linear dichroism spectra. These values of tilt and twist are close to (-40 degrees, 30 degrees), (-40 degrees, 150 degrees), (40 degrees, -30 degrees) or (40 degrees, -150 degrees), with an uncertainty in both angles of about 15 degrees. Although the linear dichroism results do not allow a choice between these possible orientations, the latter two combinations are not in agreement with earlier circular dichroism calculations. For the complexes formed with poly(rC) and poly(rA), the linear dichroism spectra could not be explained by the same base orientations. In these two cases also the protein contribution to the linear dichroism appeared to be different, indicating that for some aromatic residues the orientations are not the same as those in the other complexes. The different structures of these complexes are possibly related to the relatively low binding affinity of gp32 to poly(rC), and to a lesser extent to poly(rA).     

Sequence-dependent binding of bis-amidine carbazole dications to DNA.

The conventional wisdom argues that DNA intercalators possess a condensed polyaromatic ring whereas DNA minor groove binders generally contain unfused aromatic heterocycles, frequently separated by amide bonds. Recently, this view has been challenged with the discovery of powerful intercalating agents formed by unfused aromatic molecules and groove binders containing a polyaromatic nucleus. Bis-amidinocarbazoles belong to this later category of drugs having a planar chromophore and capable of reading the genetic information accessible within the minor groove of AT-rich sequences [Tanious, F.A., Ding, D., Patrick, D.A., Bailly, C., Tidwell, R.R. & Wilson, W.D. (2000) Biochemistry 39, 12091-12101]. But in addition to the tight binding to AT sites, we show here that bis-amidinocarbazoles can also interact with GC sites. The extent and mode of binding of 2,7 and 3,6 substituted amidinocarbazoles to AT and GC sequences were investigated by complementary biochemical and biophysical methods. Absorption, fluorescence, melting temperature and surface plasmon resonance (SPR) measurements indicate that the position of the two amidine groups on the carbazole ring influences significantly the drug-DNA interaction. SPR and DNase I footprinting data confirm the AT-preference of the compounds and provide useful information on their additional interaction with GC sequences. The 3,6-carbazole binds approximately twice as strongly to the GC-containing hairpin oligomer than the 2,7-regioisomer. The high tendency of the 3,6 compound to intercalate into different types of DNA containing G.C base pairs is shown by electric linear dichroism. This work completes our understanding of the sequence-dependent DNA binding properties of carbazole dications.     

The interaction of multiply charged intercalating heterocycles with DNA

The interaction with DNA of two aromatic nitrogen heterocycles, 1 and 2 , which at pH 6 have two positive charges on their ring systems and two cationic side chains, have been determined. A third similar compound, 3 , with a single side chain and reduced ring charge, was analyzed as a control. Viscometric titrations with sonicated DNA indicated that all three compounds bind to DNA by intercalation. Spectrophotometric binding studies as a function of ionic strength indicated that both 1 and 2 bind to DNA as tetracations at pH 6. These are the first examples of intercalators with two charges directly on the intercalating ring system. Dissociation kinetics experiments as a function of ionic strength confirmed that 1 and 2 bind to DNA as tetracations. Compound 1 has a G · C base-pair binding preference, 2 seems to prefer binding to alternating pyrimidine–purine sequences regardless of the composition, and 3 has no significant binding specificity.     

Comparative molecular dynamics simulations of the potent synthetic classical cannabinoid ligand AMG3 in solution and at binding site of the CB1 and CB2 receptors

The C-1'-dithiolane Delta(8)-tetrahydrocannabinol (Delta(8)-THC) amphiphilic analogue (-)-2-(6a,7,10,10a-tetrahydro-6,6,9-trimethylhydroxy-6H-dibenzo[b,d]pyranyl)-2-hexyl-1,3-dithiolane (AMG3) is considered as one of the most potent synthetic analgesic cannabinoid (CB) ligands. Its structure is characterized by rigid tricyclic and flexible alkyl chain segments. Its conformational properties have not been fully explored. Structure-activity relationship (SAR) studies on classical CBs showed that the alkyl side chain is the most critical structural part for the receptor activation. However, reported low energy conformers of classical CB analogues vary mainly in the conformation of their alkyl side chain segment. Therefore, comparative molecular dynamics (MD) simulations of low energy conformers of AMG3 were performed in order to investigate its structural and dynamical properties in two different systems. System-I includes ligand and amphoteric solvent DMSO, simulating the biological environment and system-II includes ligand at active site of the homology models of CB1 and CB2 receptors in the solvent. The trajectory analysis results are compared for the systems I and II. In system-I, the dihedral angle defined between aromatic ring and dithiolane ring of AMG3 shows more resistance to be transformed into another torsional angle and the dihedral angle adjacent to dithiolane ring belonging in the alkyl chain has flexibility to adopt gauche+/- and trans dihedral angles. The rest of the dihedral angles within the alkyl chain are all trans. These results point out that wrapped conformations are dynamically less favored in solution than linear conformations. Two possible plane angles defined between the rigid and flexible segments are found to be the most favored and adopting values of approximately 90 degrees and approximately 140 degrees. In system-II, these values are approximately 90 degrees and approximately 120 degrees. Conformers of AMG3 at the CB1 receptor favor to establish a cis conformation defined between aromatic and dithiolane ring and a trans conformation in the CB2 receptor. These different orientations of ligand inside the binding pocket of CB1 and CB2 receptors may explain its different binding affinity in the two receptors. The results of this study can be applied to other synthetic classical CB ligands to produce low energy conformations and can be of general use for the molecules possessing flexible alkyl chain(s). In addition, this study can be useful when restraint of the alkyl chain is sought for optimizing drug design.